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Preface 


In  this  thesis  I  have  sought  to  create  a  computer  model  of  the 
MATO  standard  continuously  variable  slope  delta  voice  encoding  system 
with  sufficient  flexability  to^  permit  continued  study  of  the  standard's 
specifications  and  tolerances.  This  investigation  has  started  the  pro¬ 
cess  of  evaluating  the  proposed  NATO  standard,  however,  additional  study 
is  necessary  to  determine  the  standard's  adequacy  to  assure  system 
interoperability . 

I  wish  to  ‘thank  my  thesis  advisor,  Capt.  Kizer,  and  the  members  of 
the  thesis  committee.  It.  Col.  Carpinella  and  Capt  Seward,  for  their 
assistance,  guidance,  and  tolerance  during  the  course  of  this  project. 

Jeffrey  A.  Lersch 
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Abstract 

A  computer  rrodel  of  the  continuously  variable  slope  delta  voice 
encoding  system  specified  in  the  draft  STAf'AG  on  "Analoguc/Digital  Con¬ 
version  of  Speech  Signals  for  Tactical,  Digital,  Area  Communications 
Systems",  dated  June  1978,  is  developed  and  implemented  in  I'OkTR/J'  IV, 
The  model's  performance  is  then  characterized  in  terms  of  idle  channel 
noise,  total  harmonic  distortion,  internodulation  distortion,  signal- 
to-noise  ratio,  and  frequency  response.  For  each  of  these  attributes, 
the  system's  performance  is  presented  graphically  and  compared  to  the 
criteria  established  in  the  draft  standard.  The  model  is  then  exercised 
by  varying  the  system  parameters  to  the  limits  imposed  by  the  standard 
and  the  resulting  performance  compared  to  the  previously  determined 
ideal  system  performance.  The  results  show  that  the  performance 
characteristics  measured  are  most  sensitive  to  the  primary  integrator 
response  and  output  filter  response  v/hen  the  system  parameters  are 
restricted  to  the  range  allov;ed  by  the  draft  MATO  standard. 


INVE!=;TIGATI0N  of  cowtinuously  vakiablk  slopk 
DELTA  (CVSD)  KODULATOR/DEI-'.ODULATOR  COFPATAbILITY 


I,  Introduction 

A  draft  f'ATO  standard  on  the  analog  to  digital  conversion  of 
speech  signals  using  continuously  variable  slope  delta  (CVSD)  modula¬ 
tion  is  presently  being  circulated  among  the  military  services  for  com¬ 
ments,  The  proposed  standard,  "Tlie  Analogue/Digital  Conversion  of 
Speech  Signals  for  Tactical,  Digital,  Area  Communications  Systems," 

June  1978,  seeks  to  assure  trarismission  systems  interoperability  by 
standard! zini'  Uic  system  architecture  and  setting  tolerances  on  l^ey 
system  parameters.  The  draft  standard  (see  Appendix  A)  consists  mainly 
of  end-to-ena  system  performance  criteria,  primarily  signal -to -noise 
ratios  and  amplitude  response  characteristic-;,  i'o  standards  are  specific¬ 
ally  established  for  transnission-end/reception-end  mismatch  performance. 

The  Air  Force  Communications  Command,  AFCC/OA,  has  questioned 
whether  the  limited  number  of  specifications  given  are  adequate  to 
assure  system  performance  when  the  CVSD  encoding  equipment  is  not  per¬ 
fectly  matched  to  tlie  decoding  equipment.  Are  the  tolerances  specified 
sufficiently  narrow  to  assure  no  serious  signal  degradation  when  the 
modulatoi’  and  demodulator  parameters  differ  by  the  maximum  amount 
allowed  by  the  draft  standard?  This  is  the  question  that  this  inves¬ 
tigation  seeKS  to  ansv/er. 

Problem  teir.enT  Determine  the  adverse  effects  on  the  transmitted 
signal  and  their  severity  w’leu  the  CV8D  encoder  and  decoder  parameters 
differ  within  the  limits  allowed  by  the  draft  STANAG  on  "The  Analogue/ 
Digital  Conversion  of  Speech  Signals  for  Tactical,  Digital,  Area  Com¬ 
munications  Systems,"  June  1978. 

Approach  The  approach  of  this  investigation  is  to  perform  a  computer 
simulation  of  the  CVSD  analog  to  digital  conversion  system  then  evaluate 
the  system’s  performance  under  varying  external  and  internal  conditions. 
Initially,  a  basic  mathematical  analysis  of  the  system  components  is 
performed  and  mathematical  models  of  the  CVSD  encoder,  decoder  and  the 


1 


input  and  output  filters  are  developed.  Tliese  models  are  then  translated 
into  computer  subroutines  and  coded  in  FORTRAN,  In  the  next  section,  the 
tests  used  to  characterize  the  model  are  described.  These  tests  consist 
of  the  standard  voice  frequency  measurements  as,  idle  channel  noise,  total 
harmonic  distortion,  intermodulation  distortion,  si^nal-to-noise  ratio, 
and  frequency  response.  The  system  is  first  characterized  with  the  encoder 
and  decoder  parameters  matched.  Each  test  is  performed  at  frequencies  and 
ariR.'litudes  across  the  normal  active  ranfie  of  the  system.  After  system 
performance  under  ideal  conditions  is  established,  the  system  parameters 
are  allov/ed  to  vary  across  the  ranges  allowed  by  the  draft  standard 
and  the  degradation  of  the  transmitted  signal  by  encoder  and  decoder 
parameter  mismatching  evaluated.  Tlie  results  of  the  testing  are  then 
analyzed  to  determine  v/hich  param.eter  mismatches  most  seriously  degrade 
system  performance  and  to  determine  if  the  degradation  is  serious  enoi;gh 
to  prevent  signal  transmission. 
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I I .  Analog /Digital  Conversion  System  Model 

Ihe  basic  analog/ dipital  CVSD  systen  defined  in  blie  draft  si .^n- 
dard  is  sho\m  in  the  block  diagram  in  Figure  1.  It  consists  of  four 
major  components,  the  input  and  output  filters,  the  CVSD  encoder,  and 
the  CVSD  decoder.  Fach  of  these  components  is  discussed  in  the  following 
sections . 

C\'SD  Encooer  Oneration  Tlie  CVSD  encoder  structure  is  shown  in  Figure 
2,  The  bandlimited  signal  from  the  input  lov/-pass  filter  is  applied  to 
one  input  of  the  comparator  and  sampled  at  the  clock  rate,  either  16  or 
32  kb/s.  Each  input  sample  is  compared  to  an  estimate  of  the  signal 
generated  by  the  encoder  feedback  network  from  previous  input  samples. 

In  this  model,  the  comparator  output  is  +1  if  the  input  sample  is  greater 
than  the  signal  estimate  and  -1  if  the  input  sample  is  less  than  the 
estimate.  This  polar  signal  is  converted  to  binary  (+1  =1,  -1  =0)  and 
forms  the  transmitted  data  signal.  To  generate  the  next  signal  estimate, 
the  polar  signal  from  the  comparator  is  routed  to  the  input  of  the  slope 
overload  detector. 

Slope  overload,  as  defined  for  this  system,  is  the  condition  when 
the  last  three  comparator  outputs  are  identical,  either  all  +l's  or 
all  -I's.  This  indicates  that  the  input  signal  amplitude  is  rising  or 
falling,  respectively,  faster  than  the  encoder  can  track  using  the  pre¬ 
sent  step  size.  Other  systems  define  slope  overload  by  different  length 
strings  of  identical  comparator  outputs.  Strings  of  two  or  four  identical 
bits  are  also  commonly  used  to  indicate  slope  overload.  The  last  three 
comparator  outputs  are  stored  in  a  shift  register  within  the  slope  over¬ 
load  detector  and  combinational  logic  circuits  used  to  determine  if  a 
slope  overload  condition  exists.  The  slope  overload  detector  output  controls 
the  position  of  the  switch  shown  in  the  block  diagram.  Under  normal 
conditions,  when  slope  overload  does  not  exist,  the  switch  is  in  the 

V  position.  Upon  occurrence  of  slope  overload,  the  switch  position 
min 

is  changed  and  V  is  applied  to  the  input  of  the  syllabic  filter, 
m&x 

The  syllabic  filter  is  generally  a  simple  single  pole  RC  filter 
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Figure  1,  The  CVSD  Signal  Processing 


whose  output  is  defined  as  tlie  step  size  of  the  CV?!D  encoder.  The  syllabic 

filter  ('onti'ols  t’no  response  of  the  system  to  the  envelope  of  thie  speech 

sipn-'-l  beinp  processed,  ri'oloncol  aopli  ;etion  of  V  .  to  the  input 

min 

of  the  syllabic  filter  causes  the  output  to  decrease  to  a  minimum  non¬ 
zero  vniue  that  approaches  V  .  .  Under  continuous  slope  overload  con- 

min 

ditions,  V  is  continuously  apnlied  to  the  input  of  the  syllabic  filter 

causinp  the  filter  output  to  increase  to  an  average  value  approaching 

V  .  The  magnitude  of  the  syllabic  filter  output  is  used  to  control  the 
max 

amplitude  of  t;;e  output  pulse  of  the  pulse  modula.tor.  Polarity  of  tlie 
pulse  is  contiolled  by  the  last  output  of  the  comparator. 

The  primary  integrator  responds  to  the  square  w- ve  signal  from  tiie 
pulse  modulator  and  its  output  forms  the  signal  estimate  used  by  the 
comparator.  At  the  end  of  each  clock  period  a  new  estimate  is  avail¬ 
able  to  be  used  by  the  comparator  in  generating  the  next  binary  output  ' 
and  the  next  r-ignal  estimate,  line  primary  i  ntevj’ator '  s  response  con¬ 
trols  th€^  maxi:,'.uiT  analog  signal  frequency  that  can  be  processed  ':,'"iX)Ugh  the 
the  eVSD  analog,  to  digital  conversion  syster:. .  In  figure  3,  are  shovm 

sample  waveforrr.s  at  each  stage  of  the  analog  to  dig.it.tl  conversion 
process . 

Encoder  Al:-.or:Phm  The  mathematical  description  of  f.ie  CVSD  encoder  oper¬ 
ation  is  laryciy  a  description  of  its  component  filters,  the  primar.v 
integrfitor  and  the  syllabic  filter.  One  of  the  system  characteristics 
specified  by  the  draft  siaricard  is  the  primary  integrator  response.  The 
impulse  resjixinse ,  in  its  simplist  form,  is  given  as, 

a(t)  =  (1) 

vrtiere 

f  ,  =  the  ool(  frequency  of  the  filter  in  hertz 

cl 

The  pri''’ary  inte"rntor  innut  signal  is  tlio  squ^u'e  vave  output  of  the 
pulse  modulator,  which  for  a  single  pulse  can  be  described  as, 

a(t)  =  0  t  <  0  and  t  >  T 

=  a  0  S  t  s:  T 


\J 


slope 

pulse  syllaDic  overload  ,  primary 

binary  modulator  filter  detector  comparator  integrator  analo; 

output _ output _ output _ output  output  output  input 


sample  rate 


Figure  3.  Sample  Waveforms  at  Various  Points  in  the  Encoder 


where 

T  = 


sample  rate 


The  primar;/  int^\  rator  output  is  determined  by  convolving  the  filter 
impulse  response  with  tne  input  signal. 


c 

x(t)  =  a(t)  *  O'(t)  =  J" a(T)  Q;{t  -  t)  dT 

—  OO 

-  0  ,  for  t  <  0 


r  -P.irr  tn 

[l-e  clJ.forOStST 


27rf 


cl 


-27rf  T-,  -pTTf  (t  -  T) 


(3) 


2wf 


^  [l  -  e  cl  ]  e  ‘  "cl ,  for  t  >  T 
cl 


Since  the  primary  integrator  output  is  of  interest  only  at  the  end  of 
each  clock  period,  when  it  is  used  for  comparison  with  the  input  analog 
signal,  the  continuous  equations  developed  above  can  be  simplified 
as  follows.  For  t  =  nT, 


(1  -  a)  ff"  "  ^  ,  n  =  1,2,  ...  ,11 

1 


(41 


where 


-27rf  ,  T 
a  =  e  cl 


k.  =  2irf 
1  cl 


Using  sunoryxjsi  ti  or. ,  iHg  primary  integrator  output  as  the  result  of  a 
series  of  U  oulses  can  be  described  as, 


—  (1  -  a)  +  (1  -  a) a  +  .  .  .  +  (1  -  a)a' 


n-1 


Eri-n 


(1  -  a)  a”  ,  for  n  -  1,2, 


.  ,  f’  C^) 


n=0 


8 


nils  exproHsion  can  also  be  defined  recurrsivcly ,  depending  only  on  the 
present  input  and  the  last  output.  Tnis  definition  can  be  used  to 
simulate  the  CVSD  encoder  on  a  computer. 


(6) 


The  analysis  of  the  syllabic  filter  output  follows  identically 
that  of  the  primeiry  inte^.rator.  The  impulse  response  of  the  syllabic 
filter  is, 


|3(t) 


r?>rt 

u 


] 


(7) 


\^^lere 

t  =  the  time  constant  of  the  syllabic  filter 
c 

Tlie  recursive  expression  for  the  syllabic  filter  output  is, 

t 

c 

exp  [-  —  ] 
c 

either  V  or  V  .  ,  the  syllabic  filter  input 
max  mm 

Encoder  Computer  Subroutine  Eouations  (6)  and  (8)  are  implemented  in 
the  subroutine  used  to  perform  the  CVSD  encoding  for  this  investigation. 
Figure  3  is  a  flowchart  of  the  subroutine  used  for  encoding  and  Appendix 
B  is  the  FORTRAN  code  used.  All  of  the  system  defining  parameters  are 
transmitted  to  the  subroutine  through  the  calling  statement,  Fno-oding 
is  performed  on  an  array  basis.  The  analog  signal  to  be  analog  to  digital 
converted  is  first  sampled  at  the  clock  rate  and  the  samples  placed  in 
the  input  array,  v;hich  is  of  size  1  x  N,  whore  is  the  n\;;:ibcr  of  samples. 
All  of  the  samples  are  encoded  by  the  subroutine  and  the  binary  data, 
s'  rc  im  placed  in  the  output  array  be  fore  the  subroutine  returns  conti'ol 


v/here 

k. 


Figure  A.  CVSD  Encoder  Subroutine  Flov;chart 


LO 


Figure  4  (continued),  CVSD  Encoder  Subroutine 

Flowchart 
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to  tile  callir.i''  proc.r^m, 

Tv/o  of  the  prirti'iry  fiyssteii  tlofinin^  paraneters,  V?-'AX  and 
must  be  fenerated  by  subroutine  VMAXOPT  (appendix  I)  before  the  encoder 
subroutine  is  called.  It  should  be  noted  tiiat  the  constants  and 
derived  in  equations  (5)  and  (8)  are  not  stieci  fical  ly  included  in 
the  propram  statements  defining  the  primary  integrator  and  syllabic 
filter  responses  but  are  expected  to  be  included  in  the  values  calculated 
for  Vl'AX  and  Vlfin, 

The  variable  DC  in  the  subroutine  is  the  duty  cycle  of  the 
slope  overload  detector.  This  variable  is  not  used  during  the  encoding 
process.  Instead,  it  is  used  only  by  VIIAXOPT  vrtien  the  valuers  of  Vl'AX 
and  VXIH  sre  being  determined, 

CT/^D’  Decoder  Oneration  Tlie  CVSD  decoder  circuit  is  identical  to  the 
encoder  feediback  circuit,  A  block  diagram  of  the  decoder  is  shov/n  in 
figure  5,  Tlie  only  dii i fer-once  between  the  decoder  and  tri-,-  encoder  is 
that  the  decod;er  has  no  comp;irator.  Tlie  binary  signal  fron,  the  encoder 
is  applied  directly  t:o  the  slope  overload  detector  and  the  output  signai 
is  takcii  from  the  primary  integrator.  Tlie  signal  estimate  generated  in 
the  decoder  is  identical  to  that  generated  in  the  encoder,  if  ti.e  parameters 
of  each  unit  are  i  deth  icul .  However,  at  the  decoder  t!ie  signal  estimate 
is  of  interest  at  all  times  and  not  just  at  the  sample  periods,  as  the 
decoder  signal  estimate  is  the  approximation  of  the  analog  signal  trans¬ 
mitted  by  the  CV’^D  encoder.  Figure  6  shova;  sar.ple  waveforms  at  various 
points  within  the  decoder.  The  waveforms  are  identical  to  those  shown 
in  figure  3,  except  that  the  decoder  primary  integrator  output  is  shown 
as  a  continuou  .  c.ignal. 

Decoder  Alforith;-!  Since  the  decoder  circuit  is  identical  to  the  ‘'ni;oder 
circuit  without  the  comparator,  the  mathematical  anrilysis  devem-(:ei  f.'r 
the  encoder  is  also  applicable  to  the  decoder.  One  excei)tion,  hov.'-„'ver, 
is  that  the  simplification  used  to  obtain  eqUcition  (4)  is  not  v  .orally 
applicable  to  the  decoder  since  the  primary  integrator  output  '  Ihe 
decoder  is  required  to  be  continuous.  The  recursive  expression  for  the 
decoder  output  is, 
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3 


I 


sanple  rate 


Figure  6,  Sample  V/aveforns  at  Various  Points  in  the  Decoder 


-  [yC  h-i]  •)]  -"'■'■i'  ^ 


(1  -  e  1  )  (9) 

h 


for  (N-1  )T  <  <  NT 


v.’here 


sample  rate 


=  the  primary  integrator  input  for  (n-l)T  <  t  <  fiT 
=  the  pole  frequencj  of  the  primary  i nt. ■'’rator  in  hertz. 


Analysis  of  ttie  d.:-o.)i\er  syllabic  filter  is  i  e..^ii'r;ical  to  cl.-.t  of  the 
encoder  syllabic  filter  and  eoi;ation  (B)  also  applies  to  the  decoder. 


^N-1^  (1-/3) 


where 


~  t~  '  syllabic  filter  output 


P  =  exp 


t 

c 


V,,  =  either  V  or  V  .  ,  the  syllabic  filter  Input 

N  max  nin 


Decoder  Computer  ""aibrontine  Usin-T  eq’jations  (6)  'u  i  (B)i  the  decoding; 
subroutine  is  ir.ipleinentc-d  as  shown  in  the  flowchart  in  fipur--  7 ,  Equation 
(9)  is  not  used  since  the  straight  line  approximation  provided  by  the 
Tilcomp  plotter  provides  a  sufficiently  accurate  representation  of  the 
decoder  outj  .ut  for-  this  i  nvesti -Tation .  Except  for  the  el  L,„ina  L  ion  of 
th'='  comparison  .step  used  in  tlie  encoder  subrou ti'ie ,  the  decoder  sub¬ 
routine  is  nr-firly  identical  to  the-t  of  the  encoder.  All  comments  applic¬ 
able  to  the  encoder  subroui'inc  are  also  applicable  to  thr-  decoder  fju’o- 
routine.  The  FORTRAh  code  for  the  decoder  subroutine  is  attached  in 
Appendix  C. 
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Figure  7,  CVSD  Decoder  Subroutine  Flowchart 
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Fncoiior  and  rnco-'Vr 


From  thc!  expressions  de vo lo c»' d  in  the 


preceding  sections  rie.scribin,^  tho  CV'-’D  encoder  and  decoder,  iv  l■:an  be 
seen  that  there  are  fonr  paraiTiOters  that  determine  the  characteristics 
of  the  encoder  and  decoder.  They  are,  f  for  tiie  primary  inte'^rator, 
t  for  the  syllabic  filter,  and  V  and  Y  .  v/hoso  value  determines  the 
magnitude  of  i.he  sten  sizes. 

Tlio  draft  standard  specifies  the  value  of  f  ^  explicitly  in  para¬ 
graph  3,2,  V'nen  the  primary  integrator  consists  of  a  single  [.vile  filter, 
the  value  of  f^^  is  required  to  be  between  100  and  300  IIz. ,  Other  poles 
and  zeros  can  be  added  to  the  primary  integrator,  in  accord:  nee  v;ith  the 

draft  standard,  however,  only  f  is  required.  In  this  invest!  f’ation , 

c] 

the  single  iiole  version  of  the  primary  integrator  is  used  in  thc  CYSD 
encoder  and  decoder  models. 

For  the  syllabic  filter,  the  draft  standard  does  not  specify  the 

value  of  t  directly.  Instead,  t  is  specified  in  terms  of  the  decoder 
c  c 

output  signal  when  a  given  input  is  applied  to  tlie  encoder.  V/hen 
the  analog  input  signal  at  000  Hz  is  stepped  from  -42  dOmO  to  0  dBmO, 
the  decoder  output  signal  is  required  to  achieve  901^  of  its  final  value 
within  2  to  4  milliseconds  after  the  output  signal  starts  to  rise.  ("OTF: 
For  this  system,  the  standard  specifies  the  reference  test  level  point 
to  be  -4  dEm .  So,  a  -42  dOnO  is  actually  -46  dBm.  All  measurements  taken 
in  this  investiga L'ion  are  stated  in  dBmO,  unless  explicitly  stated 
otherwise , ) 

The  values  of  V  and  V  .  are  also  not  specified  directly  by  the 
max  nin 

draft  standard,  but  arc  specified  in  terms  of  the  syllabic  filter  output. 

The  syllabic  filter  output,  v/hich  has  previously  been  defined  as  the 
step  size  of  tfie  encoder  and  decoder,  is  required  to  .'s  linear  as  a 
function  oi'  the  slone  overload  detector  duty  cycle.  The  slope  overload 
detector  duty  cycle  is  defined  the  ratio  of  the  number  of  tines 
slope  overload  is  detected  to  the  number  of  samples  in  the  sane  period. 

In  paragraph  3,4  of  the  draft  stand/ird,  the  step  size  is  showi-i  us  varying 
linearly  as  the  duty  cycle  ranges  from  0  to  .5.  The  step  size  ratio,  the 
ratio  of  the  syllabic  filter  output  v/hen  an  800  Hz,  0  dBmO  sigiv  1  is 
applied  to  the  encoder  input,  to  the  syllabic  filter  output  v.’hen  the 
encoder  input  is  grounded  is  required  to  be  34  d3  ±  2  dB.  This  specification 
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in  combination  with  the  specifications  for  f  and  t  determine  the  values 

cl  c 

of  V  and  V 

max  min 

Due  to  the  fact  that  the  parameters  interact  with  each  other,  the 

values  of  t  ,  V  ,  and  V  .  need  to  be  determined  recursively.  A  value  of 
c  max  min 

f  ,  is  chosen  within  the  ranee  given  by  the  standard  and  an  estimate  of 
cl 

t^  chosen  near  its  expected  value.  The  syllabic  filter  determines  the 

system  response  to  the  amplitude  modulation  of  a  voice  signal.  As  the 

highest  frequency  in  the  envelope  of  the  voice  signal  is  generally  about 

100  Hz,  t^  is  estimated  to  be  the  reciprocal  of  this  frequency  or  ,01. 

A  nominal  step  size  ratio  is  given  by  the  draft  standard  to  be  34  dB . 

These  three  narameters  are  used  bo  calculate  the  values  of  V  and  V  ,  , 

max  min 

Figure  8  is  the  flowchart  of  the  program  that  calculates  these  values 

using  the  subroutines  shov/n  in  figures  9  and  10,  then  plots  the  resulting 

syllabic  filter  output  as  a  function  of  slope  overload  detector  duty  cycle. 

Initially,  estirrated  values  of  V  and  V  .  are  used  and  the  slope 
’  max  min 

overload  detector  duty  cycle  and  system  step  size  ratio  calculated  when 

an  800  Hz,  0  dBmO  test  signal  is  input  to  the  CVSD  encoder.  If  the 

calculated  values  are  not  within  the  tolerances  specified,  V  and  V  . 

max  min 

are  adjusted  and  the  calculations  repeated.  This  process  is  continued  un¬ 
til  values  of  V  and  V  .  are  determined  that  produce  a  slope  overload 
max  min 

detector  duty  cycle  of  .5  ±  1%,  and  a  step  size  ratio  within  ,015^  of  the 
input  value , 

After  determining  the  values  of  V  and  V  .  ,  the  entire  CV8D 

max  min 

system  is  tested  to  determine  if  the  rise  time  requirement  is  met  using 

the  parameters  that  have  been  calculated.  The  flowchart  of  the  test  program 

is  shov.Ti  in  figure  11.  To  determine  the  system  rise  time,  a  test  signal 

consisting  of  alternate  series  of  500  samples  of  an  800  Hz,  -42  dBmO  sine 

wave  and  500  samples  of  the  800  Hz  sine  wave  at  0  dBmO,  The  initial  series 

at  -42  dBmO  initialize  the  storage  elements  of  the  slope  overload  detectors 

in  both  the  encoder  and  decoder  and  get  the  system  into  an  initial  steady- 

state  condition.  After  processing  the  test  signal  through  the  system, 

the  output  signal  is  plotted  in  the  vicinity  around  one  of  the  steps  in 

input  signal  power.  The  system  rise  time  is  then  determined  grn pliically , 

Figure  12  shov/s  a  sample  output  from  this  program  for  both  the  16  and  32 

kb/s  sample  rates.  This  test  was  performed  with  f  ,  =  100  Hz,  step  size 

cl 

ratio  =  34  dB,  t  =  ,01  for  the  16  kb/s  sample  rate,  and 
c 
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Figure  9,  CVSD  Gystcn  Parameter  Qilculation 
Giibroutin'-'  Flov/chark  ( 'A'n 
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Figure  10# 


CVSn  System  Parameter  Calculation 
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Fiyure  13 


Input  and  Outp'.it  I  ov.'-PanF.  The  last  of  the  najor  connonents  mkina 

up  tlie  CVSD  nnalop/dipital  conversion  syston  are  the  input  and  output  fil¬ 
ters.  These  filters  are  used  to  limit  both  the  input  and  output  sip,- 
nal  spectrum  to  the  voice  band  frequencies  only.  For  telephonic  communi¬ 
cations,  the  voice  band  is  generally  considered  to  be  those  frequencies  less 
than  3600  !lz,  lor  optimal  system  performance,  these  filters  should  have 
a  very  sharp  cut-off  and  hi;;h  loss  characteristics  in  the  stop  band. 

The  purpose  of  the  input  filter  is  to  limit  the  input  sipnal  spec¬ 
trum  to  prevent  aliasing  due  to  the  sap.plinp  process.  VAien  the  input 
signal  is  sampled,  in  addition  to  the  input  spectrum,  the  output  spectrum 
also  contains  sum  and  diffei'ence  frequency  components  centered  around  the 
sample  frequency.  If  the  input  specti'u;.  v;ere  to  contain  frequencies  very 
much  larger  than  the  desired  spectrum,  aliasing  or  interference  would 
occur  vAien  the  difference  frequencies  fell  into  the  baseband  spectrum. 

For  this  model,  tlie  input  filter  is  considered  to  be  an  ideal  lov;-pass 
filter.  The  test  signal  generator  output  spectrum  is  limited  to  the  voice 
band  frequencies  only,  with  no  components  falling  outside  that  range.  This 
simulates  a  low-pass  filter  with  zero  insertion  loss  in  the  pass  band  and 
infinite  loss  in  the  stop 

The  function  of  the  output  filter  is  also  to  limit  the  signiil  to 
the  voice  band,  hov/ever,  in  this  case,  the  components  outside  the  original 
input  spectrum  are  produced  by  the  non-linearities  of  the  processing 
system.  The  output  filter  smoothfi  t’.ie  si^.nals  and  eliminates  the  harmonic 
components  above  the  voice  band.  This  filter  may  have  the  sane  character — 
istics  as  the  input  filter  or  may  have  a  narrower  pass  band  to  improve 
perforiianc-’ .  The  filter  chosen  for  this  model  is  a  maximally  flat,  linear 
phase  syTr.metrical  finite  impulse  response  (FIR)  filter.  The  mocel  of 
this  filter  was  developed  by  J.F.  Kaiser  of  the  Digital  System  Research 
Department  of  the  Bell  laboratories.  Reference  1  provides  more  complete 
ctocumentation  of  the  filter  model.  There  are  two  parameters  th.-.it  define 
the  response  of  the  filter,  beta  and  gamma.  Figure  Id  shov/s  ths'  response 
of  the  filter  generated  by  this  program  and  wliere  beta  and  gamma 
are  defined.  Beta  is  the  normalized  center  frequency  of  the  transition 
region  and  gamma  is  the  normalized  width  of  the  region,  hortnalization 
is  with  respect  to  the  sample  rate.  This  filter  was  chosen  for  its  flat 


reiuonr-.e  i ’i  th;i;  p  ir.'j  hinci  in  ord/;r  to  -li  n  i  r-.ize  din  turbinico  of  the-  CVfD 
r  ni'o ’,'1'/  ■•.-■r  rr  :';'on:T,  oiive  tbooc  oro  ihe  i.'rjrv.ry  con’.poncr:  tn 

luider  invo:'l.i)',cttioti .  Tho  para:nc  {ors  c'r.on.oti  for  the  filler  are  /3  =  ,1875 
and  y  -  .1  for  the  IG  kb/;:  r  ile  and  P  ^  .I  and  V  ^  .1  for  the  32 

kb/s  sarnile  rate.  Table  11  lists  ttic  filter  coefficients  generated  by  the 
profTCir'  and.  figure  IS  siiov.'s  the  frequency  re3pon;;e  for  both  filters. 


I'iRurc  14  ,  Maximally  l-  lat  !•  IR  filter  h.esvionse 

Cteracteri nti c  and  Parameter  fefinition 

The  center  of  the  transition  region  for  the  16  kb/s  filter  is  3  kHz  and 
3,2  kliz  for  the  32  kb/s  fil1:er.  As  can  be  seen,  the  maximally  flat  char¬ 
acteristic  is  acl'iieved  at  the  expense  of  stop  band  loss,  Hov/ever,  it  will 
be  shown  Iti  the  performance  results  that  the  system  performance  meets 
most  of  the  criteria  soecified  by  the  draft  standard  in  spite  of  the 
poor  filter  performance. 

Filter  Suhrontines  The  filter  pro'-ran  developed  by  J,F.  i'aiser  is 
used  to  generate  the  FIR  filter  coefficients,  however,  it  has  been  mod¬ 
ified  to  be  a  subroutine  that  returns  the  coefficient  values  to  Ih.e  calling 
program  instead  of  printing  them  out.  These  coefficients  are  produced  by 
FI  TRGFJI  then  used  by  subroutine  P'llTEP,  to  actially  filter  the  si<  ual  input 
to  the  filter.  The  maximum  number  of  coefficients  that  can  'no  nroduced 
by  FITRGEN  v/ithout  pirorram  modification  is  20(1.  Subroutine  I  llTFR 
delays  the  output  signal  by  200  sample  period:;  so  that  it  has  at  least 
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TABI.E  li 


FIR  FILTER  C0EFFICIE7ITS 


16  kb/s  sariT'le  rate 


P( 

1) 
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.37*')  .'>2925 
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P< 
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P( 
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P( 
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P( 

1'^) 

r 

■ '  124  22i  27 

e( 

IS) 

= 
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e( 
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= 

.  ■>  ,19‘’7.v  .. 

p( 

IP) 
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..04329561 

p( 

19) 

= 

.  '»  2951926 

p( 

2r) 

= 

•  10  ■52524 

p( 

21) 

= 

01,  31339 

B( 

2?) 

.  3  oq-.’o  1 

P( 

23) 
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-'49  21 

P( 

24) 
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:  0175539 

P( 
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.  13  5?8>  ■: 

P( 
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"  .’}!  146  44 

P( 
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- .  /  '»  21525 

e( 
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-,  -’  )  'O  ■’'5521 
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Figure  15b,  CVSD  System  Model  Output  Filter  Fieeponse 
for  32  kb/s  Sample  Hate 
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ure  16.  Signal  Filterjnc  Subroutine  Flov.'chart  (FILTER) 
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Figure  16,  (continued) 
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200  input  sirnal  samples  can  be  used  by  the  filtorinr,  alrorithm.  Equa¬ 
tion  (10)  shows  the  filtering  expression  itnplementcd  by  the  I'UTKM 
subroutine . 


♦ 


the  nth  output  sample 

the  nth  input  sample 

the  ith  FIH  filter  coefficient 

the  number  of  filter  coefficients 


32 


III. 


Perforrr'tnco  Tests 


A  model  oT  the  continuously  variafcle  slone  d.  1 1  en  iln^  to  dig- 
ito]  conversion  system  is  constructed  from  the  component  models  described 
in  the  previous  sections.  Figure  17  shows  tlir  tent  configuration  sim¬ 
ulated  by  the  computer  model  used  in  this  investigation.  T!'.e  system 
under  test  is  show  in  figure  1.  In  this  simuljition ,  the  test  signal 
generator  is  a  subroutine  that  generates  samples  of  a  sinusoidal  signal 
that  can  be  composed  of  up  to  two  frequency  components  at  individually 
specified  amplitudes.  The  standard  test  signal  used  in  the  performance 
tests  is  an  800  !fe  sine  v;ave  at  -20  dBnO,  unless  otherwise  stated.  As 
previously  indicated,  the  reference  signal  level  is  -4  dBm.  All  power 
measurements  are  made  relative  to  this  level.  The  test  signal  is  generated 
as  an  array  of  5000  samples  for  most  of  the  tests  performed.  This  array 
is  then  processed  through  the  system,  the  output  array  of  each  system 
component  becoming  the  input  array  of  the  next.  The  final  system  output 
signal  is  then  processed  to  determine  the  various  signal  characteristics. 
System  performance  is  measured  in  terms  of  the  commonly  used  voice  frequency 
tests  as,  idle  channel  noise,  total  harmonic  distortion,  internodulation 
distortion,  signal-to -noise  ratio,  and  frequency  response.  These  tests 
are  first  performed  with  the  CVSD  encoder  and  decoder  parameters  matched, 
then  performed  again  with  various  combinations  of  encoder  and  decoder 
parameters  to  show  how  system  performance  degrades  under  mismatched 
conditions . 


Idle  Channel  rfoise  Test  Idle  channel  noise  is  measure  of  the  basic  amount 
of  noise  that  the  processing  system  adds  to  the  output  signal.  The  output 
signal  power  is  measured  vjhile  the  system  input  is  grounded.  Any  non¬ 
zero  pov/er  measured  is  the  idle  channel  noise.  Before  measuring  the  idle 
channel  noise,  however,  the  system  insertion  loss  is  first  set  vo  that  the 
standard  test  signal  experiences  no  change  in  power  after  being  processed 
through  the  system,  idle  channel  noise  is  then  measured  as, 

n=l 
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Figure  17.  Simulated  System  Test  Con 'figuration 


Fij^nre  18.  Idle  Channel  Noise  ProKran  Flo\;chart 


35 


where 


the  output  cij’na]  amplitude 
the  number  of  samples 


y 

M 


n 


Figure  18  is  the  flowchart  of  the  idle  channel  noise  test  used  to  measure 
the  CVSD  system  performance. 


Total  l-!arr,onic  Distortion  Test  Total  harmonic  distortion  is  one  of  the 
measures  of  system  non-linearity.  The  CCITT  procedui'e  for  measuring 
total  harmonic  distortion  is  to  input  a  single  frequency  test  signal 
near  the  center  of  the  system's  pass  band  and  measure  the  magnitude  of 
the  harmonic  components  in  the  output  sj;ectru;.i.  Total  harmonic  dis¬ 
tortion  is  then  calculated  by, 


THD 


■=3  • 


+  E. 


X  lOCTo 


(12) 


where 


E^i  •••  ,Ej,j  =  the  RM8  voltages  of  the  harmonic  signal 

components  in  the  output  spectrum 

E^  =  the  RMS  voltage  of  the  primary  signal  component  in  the 
output  spectrum 

N  =  the  largest  harmonic  within  the  system  pass  band 


Figure  19  is  tiie  flowchart  of  the  computer  program  used  to  calculate 
total  harmonic  distortion.  After  processing  the  single  frequency  sine 
wave  test  signal  throwih  the  CVSD  system,  the  output  signal  spectrum 
is  calculated  using  a  fast  fourier  transform  (FFT)  .  Due  to  the  limitations 
of  the  FFT,  the  standard  test  signal  is  not  used,  instead,  a  1000  Hz 
signal  at  -20  d!3m0  is  used.  The  FFT  procedure  can  only  measure  sig¬ 
nal  components  at  multiples  of  the  minimum  frequency  resolution  './hich  is 
determined  by  the  number  of  samples  in  the  FFT  window.  In  this  case,  the 
v/indow  length  vras  specified  to  be  255  samples,  which  allowed  a  frequency 
resolution  of  62,5  Hz  at  the  16  kb/s  sample  rate  and  125  Hz  at  the  32  kh/s 
sample  rate.  A  1000  ite  test  signal  was  chosen  as  being  both  compatible 
with  the  FFT  and  a  commonly  used  test  signal  in  voice  frequency  measure¬ 
ments  . 


Fjciire  19.  Total  liarnonic  Diatortion  Tea'  ProrTf.:; 
Flowchart  (THD) 
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In‘ or:~'03ul.*.  t  ion  ~!i'.tort1on  Trie  totril  harnonic  di lort  ion  tofjL  often 

does  not  .rive  a  complete  idea  of  the  system  response  non-linearities. 
Intermodulation  distortion  is  another  meanvirc  of  non-linearity  usnd  in 
voice  frequency  system,  'ihe  CCITT  procedure  of  neasurinc  interr.odula  tion 
distortion  is  to  input  a  composite  test  signal  mj/dc  of  of  t'lso  sinus, oi  d'.l 
sipnals  of  equal  amplitude.  Tlie  frequencies  of  the  sisnals  are 

separated  by  an  amount  taat  the  difference  freq-aency  is  -.vitriin  the 
TV'SS  hand  of  the  system.  Internodulation  distortion  ir.  tf.pn  (xils:  ula  sc -i 

fcy. 


whe  re 

E,.„  =  the  KMS  voltap.e  of  M.n  j  ■  :  .■•'••a 

di  f 

in  the  output,  sj.'ect.rum 

=  the  Ri'.S  voltare  <t)f  thf  *  j,-.;;  fre  :  y  ^ 

output  spectrum 

Eg  =  the  R!'‘.S  volta^-e  ot"  th'  ..econd  ‘r-'  -aen  '  ■  ■ 

output  spectrum. 


Fif.ure  ?0  shOv;s  the  flovchart,  of  the  prorram.  i.  t. 
modulation  distortion  for  the  C'J">T  sy.cts-r  ,  Tic  ; 
similar  to  the  total  harmoni  c  di  a  tort  io-.  ;:x.'.  r.-'-  .  > 
tral  components  used  from  the  Ki  "L'  out;)-,.'  -.re  •.t.c  • 
and  the  difference  frequency.  The  test  ;■  i  n.il  .  ■•••  i 
sists  of  7‘i0  iiz  and  1000  Hz  componc-,*::  .a*  ■. 


Signal-to-^'lo  i  se  hatio  ?'e'isurc::r>’i*:  Tiic  a  i ,  ■„!  I -*/i-nt,  i  r,': 

measure  of  how  accurately  the  s.yutem  bcin.c  cn.,r.ic;cr;  :  rc:  .  . 
input  sitinal ,  A  test  sif.nal  is  process  tt.rour.h  tie  s.wtc,-  ., 
inp  output  sipnal  conpxired  to  the  input  signal  after  co:;;cn;a  .■ 
the  system  insertion  loss  and  sit'nal  delay,  sn.'-,  if;  then  cal.  a 
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Ini  t  i  ,-il  i  /-Q 


Fi/^ure  20.  Interrnodulntion  Fi  n  tort  ion  Tort  I'ron.ra 
Flowchart  (  riTItllD ) 


S'! 


ri 


n-1 

X 

ri 

th<' 

ntr:  input  r: 

;i''n,al  ^iample 

'  n 

the 

ntii  outnut 

sirn.tl  sample 

n 

=  tile 

total  numr.fj 

■r  O’'  sam.rlc.'-- 

?'i;'ur-,  r  iio'Vo  th';  riov; ''..ir!:.  frr  the  ri.'nu]  -to-noi  re  ;r'o;;ra"i  Li3>;d  to 
cliarac ter i the  I’VSD  i'er i'or:e.ince  .  rtandird  600  .M;:  test 

sih.'ial  is  ass.'d  to  pcr'''or':,  the  inili  il  systo:  ':;h';:' icU  ri/.'ition  . 


T''r,- qrer cy  i  "s  -  ;o -i-s t  "tX)  ~e!r,oi:-  of  ser forni-ih  frequency 
renronr'o  nenrui’‘:'’:,en t  are  u'.e  1  in  ihi.'  ir:v.:,; s.  i.,<  t,ion  .  "t.e  fii'st  is  flat 
v/eihl'ited  r.onsure'sen. :  v.'hic:.  is  useo  to  de',-,  r:  ine  the  frequency  response  of 
the  entire  CV‘'in  s\''-to-i  sine  this  is  the  ir.etiiod  for  w'aich  the  dr^^ft  stan¬ 
dard  specifies  per fcr-  .nco  criteria.  A  second  r.etiioci  is  the  frequency 
selective  r- oa'-,urrr:,ent  of  the  re3{,x>n;;e  charaeteri  sti  t:s .  This  Method  i; 
used  to  investigate  tiio  frequency  response  of  the  CVSD  encoder  and  decoder 

only. 

Flat  weif'hted  frc'Guency  resionse  meiisurepent  is  performed  by  input- 
tinp  a  sin' le  t'requency  sine  •■•'ave  test  sipnal  at  a  constant  amplitude,  tlien 
rr.c  I uuri riq  thr  ;  y;'t''”-i  oi'.tput  si.-nnl  lower.  Tlie  output  signal  includes 
co.MTonentr,  at  fre quen : it; n  other  than  the  tost  siy.nal  frequency,  however, 
the  t't)V>‘r  of  entire  comnosilo  sir.nnl  is  Measured  without  filterinr. 

The  measured  pain  variations  are  tiien  olo  and  scalrc  m.uch  that  the 

600  Hz  measurement  is  0  63,  The  flov/chart  of  th-’’  rtoqrm^  usi.np  this 
proc;edur(?  is  sh.own  in  fip’ure  62. 

Frequency  selective  mea  sva'cr  ren  t  o  i"  frequeui.y  r-'Si-.'isp  uses  thi' 

:  ■:'e  test  pj-ocedure  ozeept  only  the  ma'',nitude  of  t.he  output  sipnal  com¬ 
ponent  at  the  ter:t  fre  quency  is  neasurt  d.  ’’’I'c  otlier  conro’ ti  is  of  t'ue 
comrt)!  it.”  output  signal  n.ri-  not  inclu.ded  in  !;hln  meu  sure  men  t ,  Tlx- 
measurement;;  are  then  scaled  and  riot  ted  sucii  t!;aL  the  1000  ilz  •noa;;ure- 
m.ent  is  0  dit.  ;''j,uire  23  is  the  f.lov.-c'nart  of  the  proy.ram  to  perform  the 
Crect’-ncy  r-xlecti'/e  mcar-.urr-rr ’nts ,  T!ie  1000  '.<7,  mea^ui'ement  is  used  as  the 
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Fipuro  2?> 


(continued)  Pronram  DGAIfl  Flowchart 


Fif,ure  22,  (continued)  Program  DGAITI  Flov/chart 
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5 


rcforenco  vali"'  since  the  fast  fouricr  transform  used  to  calculate  tiio 
output  si'.nal  spectrum  cannot  measure  the  component  £it  800  liz . 


f 
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IV.  Test  Results 


The  results  of  the  tests  described  in  the  previous  section  are 
presented  here.  Each  test  was  first  performed  v/ith  the  standard  800  Hz 
test  sif.nal,  while  the  CVSD  encoder  and  decoder  parameters  were  matched. 
This  test  characterized  the  ideal  system  performance  v/ith  the  system 
parameters  at  their  nominal  values,  flext,  the  testa  v/cre  performed  allowin 
the  system  parameters  to  vary  across  the  ratires  shown  in  Table  I  and  usin/^ 
test  signals  that  varied  in  frequency  and  pov/er  across  their  normal 
dynamic  ranp,es,  while  still  maintaining  encoder/ decoder  match.  I’inaily, 
the  test  were  repeated  a.qain  v/ith  the  encoder  parameters  held  constant 
at  one  extreme  of  the  permissible  values  and  the  decoder  parameters  allov/ed 
to  ranfe  to  the  opposite  extreme.  Each  test  was  performed  chanping  one 
variable  at  a  time  while  the  other  were  held  at  their  nominal  values. 

Idle  Channel  hoiso  The  results  of  the  idle  channel  noise  tests  are 
shown  in  figure  24.  For  each  sample  rate,  the  idle  channel  noise  perfor¬ 
mance  improves  as  the  step  size  ratio  increases.  This  results  from 
the  decrease  in  minimum  step  size  as  the  step  size  ratio  increases.  The 
output  signal  depends  entirely  on  the  minimum  step  size  v/hen  the  system 
input  is  zero  or  grounded,  Eince  the  minimum  step  size  is  defined  to 
be  non-zero,  the  output  signal  v/ill  alternate  positive  and  negative 
around  zero  attempting  to  approximate  the  zero  input  signal.  The  smaller 
the  deviation  from  zero ,  the  less  the  power  in  the  output  signal  and  the 
better  the  idle  channel  noise  performance ,  System  performance  exceeds 
the  criteria  specified  in  the  draft  standard.  Idle  channel  noise  is 
-88  dHmO  vs,  the  specified  -50  dDmO  at  16  kb/s  sample  speed  and  -97  dBmO 
vs ,  -60  d'dmO  at  32  kb/ s  . 

Encoder/decoder  parameter  mismatch  has  no  effect  on  idle  ci-^annel 
noise.  This  is  a  result  of  the  fact  that  no  matter  v/hat  ti/e  encoder's 
parameters,  the  output  v;iil  always  be  alternatirv  ones  and  zero:;  v/lien  the 
encoder  input  is  grounded.  Therefore,  the  input  signal  at  the  decoder 
will  always  be  the  same  and  the  output  signal  will  only  be  affected  by 
the  decoder  parameters.  The  idle  channel  noise  performance  under  -iiis- 
matched  conditions  will  be  the  same  as  shown  in  figure  24  v/here  the 
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Figure  2/1.  CVSD  Signal  Proce''.sinr'  System  Idle 
Channel  Poise  Performance 
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parameters  are  those  of  the  decoder. 

Total  'harmonic  Pi  ntortion  Tlie  draft  standard  specifies  no  nvaximum 
total  harmonic  distortion  for  the  CVJID  system,  hov.'ever,  it  is  r^eneraliy 
accepted  tiiat  distortion  levels  of  loss  than  2(Xi  'A’iH  not  ustKilly  be 
objectionable  to  the  system  users.  As  can  be  seen  from  the  test  results 
shov/n  in  fui'orc?  fh,  oyster:  ';«r forriance  at  the  IG  kb/s  sample  rate  exceeds 
this  limit  by  4-Pt  .  hystc-m  yier for; vince  •.'/'.'ion  the  sample  rate  is  increased 
to  3?  kb/s  iri’.proves  substantially.  Tlie  total  harmonic  distortion  level 
droas  to  apjjroximatlcy  6h',  Fifiure  26  shows  that  the  distortion  level 
at  both  sample  rates  is  relatively  constant  for  all  input  pov/er  levels 
within  tlie  normal  opei'atin;:  ranee  except  at  the  very  lov/  pov/er  levels. 

V/hen  the  encoder  and  decoder  r-aramoterr,  are  mismatched,  the  total  har¬ 
monic  distortiOTi  perfornanre  ;>'..ov/r;  r.one  dorradation  as  flrures  27,  23, 
and  2*^  irulicate.  The  l-'.rr.or.t  .-:''<)unt  of  deviation  from  tlie  matched 
system  performianco  occur;;  at  the  very  lov;  t.ower  levels  v/here  the  impact 
will  have  the  least  effect.  As  figure’  2'i  shows,  total  haiumonic  distortion 
is  most  sensi  hive  so  mis;natches  of  the  encoder  and  decoder  primary  inte¬ 
grator  pole  fre qviencies .  Syllabic  filter  time  constants  and  step  size 
ratios  have  minimal  isnoact  on  the  system  performance  when  mismatched, 
liowever,  all  h;\ve  '.he  r.ost  impact  at  L!ie  very  low  input  ~x>\K’er  levels, 

Intermodul  tiioe  Pj ‘-.tort ion  Intermodulation  distortion  performance  for 
the  system  model  with  riorrmnal  p/irameter  values  is  shown  in  figure  30. 

As  is  the  case  with  the  total  harmonic  distortion  test,  the  draft  standard 
provides  no  perfromance  criteria.  In  general,  intei'modulation  levels 
of  more  than  4-5"  will  be  objectionable  to  a  system  user.  At  the  16  kb/s 
sample  rate,  the  intermodulation  distortion  measured  ranges  from 
1  to  5%  depending  on  the  syllabic  filter  time  constant  used.  The  dis¬ 
tortion  falls  to  appro'/.imcitely  1/1  when  the  sample  rate  is  increafed  to 
32  kb/s.  Figure  31  shows  the  system  intermodulation  response  as  the 
input  signal  power  is  varied.  System  non-linearities  cause  the  distortion 
levels  to  rise  at  the  very  lov;  signal  levels  and  at  the  high  input  f»wer 
levels.  Across  the  normal  operating  levels  betw'een  -10  dBmO  and  -30 
dBmO,  the  distortion  is  generally  less  than  5".  V/hen  the  encoder  and 
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.01  .015  .02  .025 

I 

Syllabic  Filter  Time  Constant 


Figure  25.  CVSD  Signal  Processinp  System  Total  Harmonic 

Distortion  Performance  with  i.ncnder  and  Ic'concr 
Fhrametcrs  Itetched  (Test  Signal  =  1000  Hz.  -aCdi-.mO) 
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gure  26a.  CVSD  Syste.-n  Totel  Harmonic  Distortion  Performance 
vs,  Input  Signal  Pov;er  with  Encoder  and  Docoder 
Parameters  Hatched  at  16  kb/s  Sample  Pate 
(1000  Hz  Test  Signal) 


Step 

Size 

Ratio 

□ 

-  32 

dB 

o 

-  34 

dB 

A 

-  36 

dB 

ft  ■■ 

S-40.00  -ifl.co  ^.00 


■’  I - 1 - ' - - r - p - ! - T - 

-jo.M  .-;o  ufi.ca  -la.*  -b.do  -..cj  o.ao 

li^pui  FOkLR  (Dtr.O) 


F'igure  26b,  CVSD  System  Total  Harmonic  Distortion  Performance 
vs.  Input  Si.-nal  Pov/er  with  Encoder  and  Decoder 
Parameters  Hatched  at  32  kb/s  Sample  tete 
(1000  Hz  Test  Signal) 
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FiPiire  27  a.  CVST)  Syston  Tot.?i.L  Ilarrionic  Cir>tortion  Perforns.nce 
vs.  Input  Sig  nal,  Povyrr  '-/i  l.h  Finooder  and  Dc;ooder 
Stop  Size  Ratios  .iisi'.a  tciiod  at  16  kb/s  Sanp^.e  Rate 
(1000  iiz  Test  Sij'nal) 


B 


.  CVSD  System  Total  Harmonic  Dir;tortion  Performance 
vs.  Input  Sie.nal  Pov;er-with  Kiu'oder  and  Decoder 
Step  Size  Ratios  :  is:.iatchod  at  32  kb/s  S^ample  Rale 
(1000  Hz  Test  Siyna] ) 


Figure  27  b 


Flijure  28  a  .  CVSD  Systc:^  Total  Har'nionic  Dirstortion  Performaai 
vs,  Input  '11 -’rial  Pov'CT  v;ith  Rncoder  and  Decoder 
Syllabic  Filter  Tir'e  Oofiscants  ."cisrsatched  at 
16  kb/s  Sample  i'iate  (1000  Hz  Test  Sippial) 
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Fif'ure  28  b,  CVSO  Sy:;te:''’  Total  I'ai'r.oMi:'  'a  'tor'.iun  I’'','  u'- 

vs.  Input  Sirrial  rov.’'  r  '.’•'ili  :  u  o  i.-r  -’u.  '  v  j  I  r 
Syllabic  Filter  Time  Piiii,  i  m.:,i  ;u  ;  .,1  -'a  .1  1 1 
32.  kb/s  '-iir-Li.'lo  ;-;.ate  '.in;:.l) 
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Fipaire  29  r.  CVSO  9yste:i  Total  iiarnonic  r’ir.tort-o'i 

v;:,  Inr,ut  Siinnl  .o:'  '.vi  tii  '.’.n 'odfu'  Doc 

Prirnai'i'  lnt''.'pr';"c<''i'  f’olc  "rc-o/ :c'nc ioc,  DJ.firatc 
at  16  kb/a  S.';npie  icace  (lOO'')  iiz  lert  Sipn 
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Fipure  29  Vj  ,  c.VSD  Sya'-R'-i  T'otal  'iarvumi'-  Diritort  ioa  2  .  i -'ruice 
vr; ,  Input  Oj  .-nal  Pea-,'  v;i  1  i  Pneoi-'-r  :.iki  ■  ,  i:;!' 
Prini.'it'y  In  tc/rator  Pole  i’'requ'’''.i'io-:  -  ch-  : 

at  32  'I'.’n/K  'vmple  Ifcitc  (lO'O.'  il.-:  Ti’iit  2j  i,>l) 
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''i.i^ure  ^la  .  CVSD  System  Internodula  tj  on  Pistortion  Perforrran 
vs.  Input  Sinnal  Pov/er  v.'ith  Encoder  and  Decoder 
Parameters  ."atchcd  at  16  kb/s  (750  and  1000  Hz 
Test  Signal ) 
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I'irure  Sib,  CVSD  System  Intermodule  *  ion  ris.tortion  P'^n 
vs.  Input  Sirnal  Power  .  ■  ih  Encoder  an'i  Pc 
Parameters  Piatciied  at  i.:  'r.I/s  S;ir,;'le  Katt- 
(750  and  1000  i!z  Test  Si  -iiai) 


and  decolor  are  r.ii  sniatclie  d ,  the  i  n  teririodviln  ti  on  distortion  perfor¬ 
mance  sliov.-r;  the  same  characteristics  an  the  total  iiariixDnic  distortion. 

The  cten  ni;-,'?  ratio  and  syllabic  filter  time  constant  have  a  minir.al 
im.mct  on  system  performance  as  show  by  the  results  in  fipures  32  and 
33,  V/hon  the  priiiary  intc^'rators  htive  difforinr,  p.ole  frequencies,  tlie 
internodulation  distortion  m.easurenents  show  more  cieviation  from  the 
matched  sysler.  torformance  as  shown  in  figure  3,!. 

dj  eng  1  -  to ise  11a tj o  Figures  3h  to  <13  show  hov/  the  sysi>'i,i  v.  t.-l 
siynal-to-roise  performance  changes  v;ith  variations  is  system  parameters 
and  input  test  signals.  The  SMR  performance  under  matched  conditions 
with  the  standard  test  siynal  shov/s  very  little-  variation  with  differing 
values  of  step  six,e  ratio,  syllabic  filter  time  constant,  and  primary 
integrator  vole  frequency,  Si,qnal-to-noiGe  ratio  vs,  input  frequency 
performance  meets  the  criteria  set  by  the  draft  stnadard  across  most  of 
the  voice  band,  Fncoder/decoder  mismatches  of  steep,  ui;:.-  c  tio  and  syllabic 
filter  time  constant  have  very  little  impact  on  system  p'orformance  ,  A 
mismatch  of  the  primary  integrator  pole  frequencies,  however,  liave  a 
much  larger  effect  on  system  performance.  The  Fpr,  is  dep'raded  below 
the  criteria  set  by  the  draft  standard,  with  the  largest  deviation  fr'oi:' 
th.e  ideal  performance  occurinp  at  the  iov/er  froquonci.?-; ,  At  the  16  kb/s 
sample  rate  the  SriR  in  lov/ered  by  about  5  dB  and  eil;  the  32  k.b/s  speed, 
ahoai;  4  df , 

fip.p.  l-t<j-no  ise  ratio  per  fo  rma  nr  vs,  input  signal  px3Wer  fails 
bo  meet  the  criteria  est'-iblished  in  the  draft  standard.  At  input  levels 
belov/  a nr.roxi irately  -iO  d3m0  the  model's  performance  falls  below  the 
desired  level.  Trends  in  system  performance  as  '.he  result  of  variations 
in  the  system  pare  .eters  c.an  be  observe.''  i..  spite  of  this  poor  :  crf-ii'-rdnce , 
a-.  i  ri.'  i'j'-  A~:  'lO  "3,  fetched  system  (performance  minimal  change 

as  the  system  paracinter-;  ir.-  .  .ried  across  the  ranges  '.rpecified  in  Table 
I,  '.fnen  tlic  encoder  and  decoder  are  not  matchi^d,  SHR,  like  to  1  li  'cronic 
distortion  .'id  ''I'.er.  O'Csla’.  ion  difitor  Lioii ,  .'.vs  lii.Lh  ..  .  Tixi  i 

t.'  "atched  s/.  jierformance  except  at  the  very  low  signal  1.  ’/  Is, 

33';  s'.o'us  the  me;,'.  c::  .'.",e  ."rorn  ideal  system  performance  wiien  the  printrry 
inter, rotor  nol-.  frequencie's  are  mismatatied . 

5fl 
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Fif'-ore  3?.b .  CV3D  Syntem  Intermodulation  Distorti  on  ['crForMance 
vs.  Input  l>i‘’,nal  Power  with  F.ncodor  and  r.co-'.er 
Step  Size  Ratios  ilismatcliod  at  3"  kb/u  S  r,";  !'.-  late 
(750  and  loOO  Hz  Test  Sii-nal) 


59 
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Fipure  33ei  .  CV3D  System  Intermodulation  Distortion  Performance 
vs.  Input  Sjt'nal  Pov/er  v/i  th  Encoder  and  Decoder 
Syllabic  Kilter  Time  Constants  ibsmatched  at  16  kb/s 
Sample  Rate  (750  and  1000  Ha  Test  Sipnal) 


FiRure  33b,  CVSD  System  Intermodulation  Distortion  Ferfor^nance 
vs.  Input  Signal  Pov/er  v/ith  Encoder  and  Decoder 
Syllabic  Kilter  'J’ime  Constants  Kismatchsd  at  32  kb/s 
Sample  Rate  (750  and  1000  Ha  Test  Sirnal) 
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Figure  34a,  CVSD  System  Intermodulation  Distortion  Performance 
vs.  Input  Signal  Power  with  Encoder  and  Decoder 
Primary  Integrator  Pole  i’x'oquencies  .''ismatched 
at  16  kb/s  Sample  Rate  (750  and  1000  !{2  Test  Signal) 


Figure  34b,  CVSD  System  Intermoduln tion  TJistortion  Performance 
vs.  Input  Signal  Power  with  Kneoner  and  Decoder 
Primary  Integrator  Polo  Prcquencios  Picrvi  tclird  at 
32  kb/s  Sample  Rate  (7‘oC  and  1000  T<?3t  Signal) 
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Signal-to-Moise  Ratio  (dB) 


I 


Syllabic  Filter  Time  Constant 


Figure  35  ,  CVSH  Signal  Processing  System  Signal-to-Pois.e 

Performance  with  Fncoder  and  Decoder  Parameuers 
Matched  (Test  Signal  -  300  Me,  -POdP.mO) 


t 


CVSD  System  Si ,'->na].-'';o-i'loi se  Perfornance  vs. 
Frequency  v'itl\  Enoodor  and  Decoder  Parameters 
Itetched  (-20  dPrp.O  tost  ci£:nal )  at  16  kb/s 
Sample  Rate 


CVSD  System  Si.qna  1  - ■:o-nQi no  Per foiiyisnci-  , 
Frequency  with  F.ticod<-;r  and  Decoder  Paranotern 
hatched  (-20  dBmW  test  signal)  at  32  kb/r> 
Sample  Rate 
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Figure  37a.  CVHD  System  Sirnal-to-'loise  PerfoD'p.nce  vs, 
"requency  xith  Encoder  and  Decoder  Step  Size 
Ratio  :;i snatched  at  16  kb/s  Sample  Rate 
(-20  dBmO  test  signal ) 


B 


'igure  37b.  CVSn  System  Signal -to-floise  Perfonijncc  vs. 

hrcqi'.ency  wi^h  I''ncod<'r  and  Decoder  Step  Size 
Ratio  :-iisiiia ;  ehed  at  ‘-2  kb/s  Sarnplc  Fiatc 
(-20  dBmO  test  sign  il  ) 


64 


UNO 


Figure  3g  a.  CV'^JO  '^ivstern  Si  ^’nal-to-'loire  Porforr,  anc'^  vs. 

r’requcncy  Encoder  and  De-;od&r  Priirary 

Integrator  Polo  i'roq  none  ins  "'isriatched  at  16  kb/s 
Sanplo  Rate  (-20  dBr.O  Test  Sirnal) 


Fif'ure  30  b,  CVSD  System  Sirnal-to-'loise  Performance  vs. 

Frequency  with  encoder  and  Decoder  i'ri'-.  a'y 
Inter,rator  Pole  Frequencies  Minna tc!icd  n  t  Id /s 
Sample  P^te  (-20  dBrnO  Test  Sipnal  ) 
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Fipurc  40a,  CVSD  System  SifTnal-to-i'bi  se  F'erfonnance  vs. 

Input  Signal  /;rr;:lltudc  ’.■•i  t,h  nncoder  and  rocoder 
Parameters  hatched  at  IG  ■Go/s  'V'mplc  Rate 
(800  ife  test  i  pnai  ) 


B 


Figure  do  V).  CVr.O  System  Sir,nal-to-r>oir’,o  Performance  vs. 

Input  SLrnal  Amplitude  v/i  th  F.ncodcr  and  Decoder 
Parameters  "a  tehee  at  32  kb/s  Oar.iile  F\a  ke 
(800  llz  test  signal) 


SNR  (OB 


t  rpr; ■ n,- i'  Tho  r:y:;1:c::i  frequency  rc::',X)t..;c-  eh:.' rat' t, i  cr:  arc 
spocjflo!  iti  th'  era  ft  :;eanear'i  v;l.':n  t  C'a;'i.u'er;  by  U.e  flat  v:oj."!ifod  r.ctriO';. 
The  reeultr'  of  t;ie  tertr.  perfor;:ieiJ  \i.sin,q  tbrs  tec/inicue  are  rho'.'.’n  in  fir- 
uro::  44  'txi  '-.I,  Tla  coapufor  n'ociei'r  perforr;ance  coinplier.  with  'l-he  np<'.ci''ie' 
cri''  t  ria  oxet  rt  at  tfie  lop  one  of  the  voice  ban:';.  iiic  riiodei'a  r(;opx3nr;c 
roilcj  off  r-.hrirpjy  at  appi'oxiriately  3  VY.-/. .  V.hile  the  cli'a  f  t  ainrifl.ar'i  alloap 
sor.’.c  I'oil-oi'f,  tlic  rer.ronoe  in  not  .'-lio'A'ed  to  brt'af.  .^.harp'ly  iin<il  G 
Varie.tior.R  in  the  nysten  ;jorameterr.  li.'tve  very  little  effect  on  the  renrxsnne 
charac teri tic s  '.■fnen  the  enco.der  and  decoder  are  ::etchi';d.  Under  t.in- 
p.atched  condition;;,  frec".ency  renponso  p'erforrtjnce  follow;;  the  sar.o 
pattern  fin  that  rn'L.-ibiint-.ed  in  the  previously  desciibe;d  tents.  Tic  step 
size  ratio  and  aylla'oic  filter  tine  cons'k^ani  have  tlnital  inpaci  ,  wliiJe 
the  priifary  integrator  pole  frequency  rauses  increased  deviation  fron 
the  liiatched  syster;  pei'forrv.nce . 

Usinq  tie  frequency  selective  ncasurer'ioti t  technique,  the  rosionse 
characteristics  of  the  CVSh  encoder  and  deco  h-T  connected  bach-to-Dack 
without  the  input  or  output  filter:;  v/ere  ncasi.ireci.  'ilie  rc  5jul.ts  are  nh.o'^.'n 
in  figures  4.8  to  t]  .  These  tests  show'  that  the  response  rolls  off 
at  the  sanvile  freiquency  at  both  the  lb  a>id  34  kb/s  aar'.pie  rfites.  At 
4  kJIz  for  the  it  i:b/s  sanple  rate  and  H  viiz  for  i,he  3.'’  ki;/:-;  samp'le  rat.e 
the  i'rcc'uency  rGSi'onso  cui  ves  breap.  sharpl;';.  Inc^dor  find  necocler  nis- 
natci.  has  rra-ctic;r.ily  no  effect  on  tho  frequency  responr-.c  of  these 
syster.  coriwonents ,  ilie  rruaary  integrator  pole  frecp.iency  shows  sl.iplitly 
ore  effect  <jr',  the  re:;ron;;e  than  tiio  other  parainct.ers .  its  effect  is 
nainlv  at  t:.c  '/ery  Jov/  frequencies  in  the  voice  band. 
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Figure  'J?  a,  CVSP  Systen  a'viii  Variation  vs,  Frecuency  vn  th 

Encoder  and  Decoder  Step  Size  P>atio3  ."isratched 
at  16  kb/s  rarr.ple  Rato  (-20  d>3n0  Test  Signal) 


c:VSD  Systerr.  Sain  Variation  vs  .  Freqvior.cy  v/i 
Encoder  and  Decoder  Step  Siz.e  Ratios  .'i  sz-i tched 
at  32  kb/r;  Stcnple  Kat.o  (-20  dl’-p.O  Test.  ^ii-  ni'I) 


Figure  4‘i  b  . 


frequency  IHZ) 


Figure  47a.  CV'^D  Systcni  Gain  Variation  vs.  Frequency  with 
Encoder  and  Decoder  Primary  Integrator  Pole 
Frequencies  Mismatched  at  16  kb/s  Sample  Rate 
(-20  dBmO  Test  Signal ) 


FREQUENCY  JHZ) 


Figure  47b,  rvSD  System  Gain  Variation  vs.  Frequency  with 
Encoder  and  Docod.-‘r  Fri.’’yir'y  Integrator  Pole 
Frequencies  I'ismatcliecl  at  32  kb/s  Sample  Kate 
(-20  dOmO  Test  Signal) 
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Figure  48  a , 

CVSD  Fn coder /Decoder  nack-to-''’eck  Frequency 
Response  Performnce  v.'i  tb  !Vitched  Para;neters 
at  16  kb/s  Sanple  Ibite  (-21!  dPtriO  'I'est  Signal) 

“^tep  Si7,(’  i.'dio 

n  -  33  d-' 

O  -  34  dH 

A  -  3G  du 


FREQUEKCT  IHZ) 

Firiure  4B  b.  C’VSD  Fncoder/Decorier  Bnck-to-Bacl';  Froquency 
Response  Performnce  v/ith  Ihtchen  Para'iot'Tr, 
at  3?  kb/s  “^'lample  Rar.e  (-?.U  dBniO  Test  Si  -nal) 


FREOUEKCt  (HZ) 

Figure  50  b.  CVSD  Encoder/Decoder  Back-to-Hack  l^rooiiency 

Response  Perfornance  v/ith  SyllaDic  I  1':.  n.e 

Constants  Hisnatched  at  3?  kb/s  Sample 
(-20  dBmO  Test  Sipnal ) 


Figure  51a,  CVSD  Encoder/ Decoder  Eack-to-3ack  Frequency 
Response  Perforrianco  with  Primary  Interrsitor 
Pole  P’requencies  I-iismatched  at  16  kb/s  Sample 
Rate  (-20  dBinO  Test  Signal) 
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Figure  51b.  CVSD  Kncodcr/Decoder  Bacl;-to-Back  Frequency 
Response  Performance  v/ith  I^ir'.ary  Tnueeralo;’ 
Pole  Frequencies  Mismatched  at  32  kb/s  '  impie 
Rate  (-20  dBmO  Test  Signal) 
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V,  Conn]  u;:  i  ons  and  Rccorrnond.ition:; 

The  test  results  show  ttiat  the  computer  model  'iK.ei.s  most  of  tlie 
perforrance  criteria  set  by  the  draft  standard  when  the  system  rxira: ter 
are  mtched  an  '  ^aI  the-ir  nomi*.  il  values.  The  output  filter  Joes  not 
]iave  the  stop  band  loss  characteristics  specified  in  tlie  standard  and 
as  a  result,  t-.e  system  performance  is  marginal,  hi ynal-to-no i se  rh.ios 
fall  belo'w  t'ne  established  criteria  v/lion  the  input  po'wer  is  le;'s  thri 
-10  dBmO.  In  spite  of  this,  the  general  effects  of  variations  in  '•'le 
system  piirameter  valuer:  and  encoder/decoder  mismatches  can  be  observed 
in  the  test  rernilts. 

1.  then  the  encoder  and  decoder  are  r:atciicd,  ci,an;:c;.’  in  the  vr v 
siiie  ratio,  .-.yll'.bic  filter  time  constant,  and  the  primary  intonr.itor 
pole  freciuency  values  within  the  tolerances  allcvred  by  the  draft  standard 
have  a  ne.  pi  i ra ble  effect  on  the  transmitted  Ki.rnal, 

2.  If  tiie  step  si^e  ratio  or  the  syllabic  filter  time  constant 
are  not  the  same  in  both  tlie  encoder  a'id  the  decoder,  the  effect  of  t’ne 
mismatc!)  on  t'io  transmitted  sifpnal  is  neplir.a'ole  except  at  input  r'ower 
levels  less  t’nan  -32  dBmO,  At  these  levels,  the  effects  v/ould  not  be 
noticeable  to  t’ne  system  users, 

3.  System  performance  is  most  sensitive  to  encoder  and  decoder 
primary  integrator  pole  frequency  raisnfiatches ,  All  the  perfor-'ance  tests 
show  a  larper  dovi.atiori  from  the  matciied  sysuem  performance  vdien  the 
primary  int' ■y.rators  are  niisn.atc’r:cd .  Ttiis  tyro  of  p.aro,;:ieter  :iiismatch 
dominate::  misma  tc’nvn  of  tii'A  oliicr  parameters. 

4.  The  I'i'cequoncy  rosncinso  of  ill"  r.ye',.  n,  is  d-  ':or:nin"d  largely 
'ey  the  output  filter  v/hen  the  pafjs  band  is  restricted  to  less  tl'.an 

A  the  sa'^ple  rate.  /.hove  'A  the  sample  rate,  the  resjxjnse  in  deterii'inod 
by  the  CVSI)  encoder  and  decoder  puss  bind.  '.nie  frequency  selective 
m.easure'  if  nt  of  the  encoder/. ’.ecoder  r  unmanse  .-.nnv;:  I  hut  t'ne  :.  y- . is 
incarable  of  :::eetin.v  the  draft  standai'd  ..qain  variation  vs,  frequericy 
criteria  'iven  in  fipure  7a  of  Appendix  A.  T'.ie  cricoder  and  decoder  alo'ie 
have  a  re.  monne  i'.alls  off  sharply  a'oove  d  klh'.  for  t]ie  16  kb/s  samnle 

r.nte,  '.7.h!ie  I'e.e  i^tandard  requires  tiiat  tfie  rcsjxanse  not  fall  off  more 
t’nan  5  d"  unt.Ll  6  p;;v,  i;-  roac'ned. 

b.  T'.e  .;pf  <:  i  f  i  c  u  t  ion.  .ui  toler.inees  '’iven  iri  t];e  .Irc.ft  siandarJ 

?■><  ( 


r'  "  •  r  1  t'j  C'i  iv-rror'v.nt''^'  vc  i •/.!  . 

71.'?  t.r- it; :  'li 1 '■  a  s;  tvil;;  v/i  1  ]  t;  .!':'”-  .M-:.:  a'itt'.r.jtiati.on  line  lo  i'af.'- i,‘ 

f  i  ;;■  !  t'i  1' r;  ii'l.v.'tteM  in.i‘  eneofler  titui  i,''‘crit:lor  .  In  mo:.^  oir,ei  ,  Ih',-  cio  -r’tndi '  i 
io  riiiii  i.'I,  !  .o  v.'f  vi  j  r  ,  , ;  i i i  oi;:  i’  .iiti;  vould  be  necessary  to  rleterrnmc 
if  the  systen  resronso  wo’.ild  coritinue  to  neet  the  draft  atandar  l  critci'ia 
un  do  T'  n  i  s  ntn  i.c  i  te  d  conditions. 

Reconnen  da  ti on  s 

1,  Since  this  model's  per  romance  is  marrinal  under  ioiril 
condi  tio'in  ,  tsnr.wle;'  i ; i : -itn tch  causes  the  ncrfortmnce  to  fall  belov; 

the  criteria  set  in  the  draft  standard.  Testin?,  should  be  repeated  usin,' 
a  filter  that  has  higher  stop  band  loss.  The  additional  testing  siiould 
concetitrate  on  primary  integrator  response  is  mate  lies  between  the 
encoder  and  decoder,  since  the  otiier  paramotors  have  little  effect  on 
the  syten  response. 

2,  System  tolerance  to  bit  errors  in  the  transmission  system 
was  not  tested,  A  mismatch  between  "-.-'o  encoder  and  decoder  may  cau.-.e 
increased  sensitivity  to  transmission  errors,  Testinr,  to  establish  the 
system  response  to  transmission  bit  error  rate  my  be  desireable. 

3,  The  draft  standard  specifications  for  pain  variation  in  the 
rer.ion  botv/een  h  kl!z  and  C  kii;-'.  fo;-  the  IG  hV'i/s  sample  rate  should  be 
modified.  Fer  fort-aa'.*';  si.oui.i  be  aliov.’n-d  l.o  rol'  fiff  as  sl'iaroly  'is; 
possible  J ;  lo  ■/i;-  1 .3  k",: , 

k.  Testing,  wss  p-  r formed  u-in.p  continuous  r.inusoidril  te;  1  ?i,,nalc 
only.  If  pi  lasi-ana  lo  ei  nals  are  expectci.i  to  be  use.i  on  t'nc  sysie;;, 

test!';.  J  'pe  retie^ited  ui-'-jn.'  l.his  tyr  ■  of  r,,l.;p;al  a  vc.riou,.  i'i^i.'.. 

Vf-yir,-  ^  Piece  L'.'.e  .'’.'33  Vorll  .:-  '  ponds  on  the  iiir.h  snnnle  to 

sapoli--  ■  'jrr-’l  P.io'i  c;  it' ;e '' i  c  t  ics  of  !.i;o  Vvi.ee  sipn-'il,  q((asi-.;n;;lo~ 
si.'ricl;.;  cci';  b"  c- .xpe  c  t  e  d  i-O  suffer  nore-  f.e. '.ru  d.  1 1  ion  n'3  t'uc  pc  •  f .  t.  of 
en  ;o  I't/ de.cod'  --  i lc  .a  '  chc s  .  i'O”:-  ri 'V.tr  ic  I  ive  to!  en  mo-.  i  he  tiec- 

e:,.vr.'  to  ’i.;  'tre  n ’.^■•ot.i  i  "  rcrl'orrv.nce  '/i  th  these'  slenals.  In  .a  sPii  lion  , 
nmra  ritundmre  :  '.y  be  neiCe.nitiry,  sucii  as  delay  iis'oi’tion  r.''^cc  j  fi  ca  t  ioni;  . 
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1.1  Analogue/di^iital  conversion  of  telephone  signals  (speech  or 

other  voice-band  signals  shall  be  performed  by  a  delta 
coder/ decoder  using  syllabic  companding  controlled  by  a 
three  bit  logic . 
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2,  Foui'-vfire  to  Foui — wire  Audio  Frequency  Chara'^teristics 

2 .1  Relative  Level  at  Points  A  and  3 

The  relative  levies  at  points  A  and  B  shall  be  -4  dBr. 

2.2  The  absolute  level  is  calculated  by  the  equation  dBm  =  dBr  +  dBmO. 

2.3  Impedance  at  Points  A  and  B 

The  nominal  value  of. the  impedance  at  points  A  and  B  shall  be 
600  ohms , 

2.4  Return  loss  at  Points  A  and  B  against  600  ohms 

The  return  loss  at  points  A  and  B  shall  be  &  16  dB  in  the 
frequency  ranfje  from  300  Hz  to  3400  Hz  against  a  load  resistor 
of  600  ohms  with  an  input  level  of  -20  dBraO. 

2.5  Symmetry  at  Points  A  and  B 

Points  A  and  B  shall  be  balanced  and  not  referred  to  ground, 
i,e.  shall  be  floating. 

3,  Details  of  the  Coder  and  Decoder  Circuits 

3.1  Input  and  Output  Audio-  Filters 

For  frequencies  above  6  kHz,  each  filter  shall  have  an  atten¬ 
uation  of  a  25  dB. 

3.2  Frequency  Response  of  the  Principal  Integrator 

The  ideal  amplitude  frequency  characteristic  between  points  F 
and  G  is  shown  in  Figure  3, 


f,  =  100  Hz  to  300  Hz 

0 

1 

d3/octave  f2  (optional)  W  2700  Hz 

dB 

fV 

f  g  (optional )  RJ  6400  Hz 

k 

\ 

12  dB/octave 

6  dB/octave 

f 

f„ 

f 

1 

2 

3 

f  — 

- 

Fig,  3  -  Ideal  Amplitude  Frequency  Cliaracteristic  of 
the  Principal  Integrator 
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3,3 


3.4 


3.5 


r 


todulation  Level 


A  signal  of  800  Hz  and  0  dBmO,  applied  to  point  A  of  the  coder 
shall  give  a  duty  cycle  (mean  proportion  of  binary  '1'  digits 
at  point  D  each  one  indicating  a  run  of  3  equal  bits  at  point 
C)  of  c ,  =  0,5  at  point  D  of  the  modulation  level  analyzer 
MLA) . 


Compression  and  Expansion 


In  the  Coder  and  decoder  the  quantizing  step  size  q  which 
drives  the  principle  integrator  at  Point  F,  shall  have  an 
essentially  linear  relationship  to  the  duty  cycle  at  point  D 
of  the  MLA  integrator  (see  Figure  4), 


Nominal  Values 


-  Limiting  Values 

(curve  as  nearly 
linear  as  possible 

q  =  Quantizing  Step  Size 
°  at  0  dBmO 


1  +  98  c^ 
50 


Fig,  4  -  Relation  between  MIA  Output  Duty  Cycle  and 
Size  of  Quantizing  Steps 


It  follows  that  the  ratio  of  the  quantizing  step  size  at  point 
F  corresponding  to  a  ducy  cycle  of  c^  =  0,5  at  point  D  of 
the  MLA  integrator  at  the  minimum  step  size  q  shall  be 
34  dB  (provisional  tolerance:  +  2dB),  ° 

Companding  Speed 

The  following  is  valid  for  the  condition  that  C  is  connected  to 
C .  VAien  an  800  Hz  sine  ave  signal  at  point  A  is  suddenly 
changed  from  -42  dBmO  to  0  dBmO  the  output  signal  at  point  B 
shall  reach  90%  of  its  final  value  within  2  mS  to  4  uS. 


NOTE:- 

The  MLA  integrator  circuits  of  the  coder  and  decoder  shall 
have  the  same  characteristics  and  hence  the  same  companding 
speed, 
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Procodure  for  Testing  the  Delta  Decoder 

The  test  bit  sequence  generator  is  connected  to  the  decoder 
input  point  C  (see  Figure  2). 

Testing  is  performed  by  means  of  periodical  test  bit  sequences 
(listed  in  Table  1)  which  result  in  audio  signals  at  800  Hz  at 
the  decoder  output  }X)int  3.  The  800  Hz  levels  at  point  B 
shall  conform  to  the  values  given  in  Table  1 . 

When  the  signal  at  point  C  is  switched  from  the  periodical 
test  bit  sequence  bo  the  periodical  test  bit  sequence  g,  then 
the  output  signal  at  point  B  shall  reach  90%  of  its  final 
value  within  5.5  mS  to  11.5  mS.  I'lhen  the  signal  at  point  C 
is  switched  from  the  periodical  test  bit  sequence  g  to  the 
periodical  test  bit  sequence  a,  then  the  output  signal  at 
point  B  shall  reach  10%  of  the  value  of  the  periodical  test 
bit  sequence  g  v/ithin  4  mS  to  8mS. 

NOTE for  clarification 

For  an  RC  circuit  in  the  MM  integrator  with  time  constants  of 
4  mS  fbr  both  charging  and  discharging,  the  envelope  charac¬ 
teristic  of  the  output  signal  at  point  B  is  shown  in  Figure  5. 

For  the  case  of  switching  the  signal  at  point  C  from  the  seq¬ 
uence  g  to  sequence  a,  the  amplitude  at  the  beginning  of  dischar¬ 
ging  is  at  the  first  moment  after  switching  higher  -  by  a  factor 
of  50  -  than  the  final  value  which  is  reached  asymptotically. 

The  final  value  equals  -42  dBmO,  i.e.  0,00794,  the  amplitude 
at  the  beginning  of  the  discharging  is  hence  0,397  (c^  =  0), 

The  value  of  10/1  is  then  reached  at  5.76  mS. 


Fig,  5  -  Envelope  Cliaracteristic  of  the  Output  Signal 
at  Point  B  (Half  the  Envelope) 
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4 .  Flectrical  Performance  at  Points  A  and  B 

4.1  General 

The  required  values  under  4.2  to  4,8  are  valid  for  the  condi¬ 
tion  that  C  is  connected  to  C’ . 

For  measurement,  the  input  (point  A)  and  the  output  (point  B) 
are  to  be  terminated  with  600  ohms,  and  signals  v/hose  fre¬ 
quencies  are  sub-multiples  of  the  sampling  rate  shall  be 
avoided.  Accordingly,  where  a  nomincil  test  signal  frequency 
of  800  Hz  is  indicated,  the  actual  frequency  shall  be  slightly 
different;  a  preferred  value  is  820  Hz,  but  frequencies  from 
804  to  860  Hz. 

The  measiirements  according  to  Sections  4,2  to  4.5  shall  be 
performed  selectively. 

4.2  Insertion  Loss  between  Points  A  and  3 

The  insertion  loss  between  points  A  and  B  at  800  Hz  with  an 
input  level  of  0  dBmO  shall  be  0  dBnO  -  2  dB,  The  insertion  loss 
contributed  by  the  transmit  and  receive  sides  shall  not  ex¬ 
ceed  one-half  of  the  value. 

4.3  Attenuation  Distortion  v/ith  Frequency 

The  attenuation  distortion  relative  to  800  Hz  measured-*with  an 
input  level  of  -20  dBmO  applied  to  point  A  shall  be  within  the 
limits  of  Figure  6.  The  distortion  contributed  by  the  trans¬ 
mit  side  alone,  measured  at  point  G  of  the  coder,  shall  not 
exceed  the  limits  indicated  by  the  broken  lines  in  Figure  6, 

4 .4  Variation  of  Gain  with  Input  Level 

The  deviation  of  the  output  level  compared  with  the  value  at 
-20  dBmO  shall  not  exceed  the  limits  given  in  Figure  7  for  a 
frequency  of  800  Hz, 

4.5  Idle  Channel  Moise 


Idle  channel  noise  at  16  kbits/s: 

The  idle  channei  noise  at  point  B  shall  not  exceed  -45  dBmOp. 
The  level  of  any  single  frequency,  measured  selectively,  shall 
not  exceed  -50  dBmO  in  the  frequency  range  from  0.3  kHz  to 
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Idle  channel  noise  at  32  kbit/s 

The  idle  channel  noise  at  point  B  shall  not  exceed  -60  dBmOp, 
The  level  of  any  single  frequency,  meaured  selectively,  shall 
not  exceed  -65  dBmO  in  the  frequency  range  from  0.3  kHz  to 
16  kHz. 

4 .6  Variation  of  Quantization  and  Harmonic  Distortion  with  Input 
Level 

The  distortion  shall  be  mesured  unweighted  with  a  sinevra.ve 
test  signal  at  800  Hz.  V/ith  such  a  signal  applied  to  point 
A,  the  ratio  of  signal  to  distortion  power  at  the  output 
point  B  shall  be  a  xjve  the  limits  of  Figure  8. 

4 .7  Variation  of  Quantizing  and  Harmonic  Distortion  with  Frequency 

The  distortion  shall  be  measured  unweighted  with  a  sinewave 
test  signal  of  -20  d.BmO.  With  such  a  test  signal  applied  to 
point  A,  the  ratio  of  signal  to  distortion  power  at  the  output 
point  B  shall  be  above  the  limits  of  Figure  9, 
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Test  Frequency  800  Hz 


f 


Fig.  9b  -  Quantizing  and  Harmonic  Distortion  with  Frequency 
at  a  Digit  Rate  of  32  kbit/s 
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APPErroix  B 

CVSP  mcodin:"  ^.ubrou^in-^. 

«j*<»ooTn«  o«<»cui»^uT,oumrr.n,rf.rct,rc«.FC3.Tc.u»w,v»»iH,DC) 

• . -CWSO  CNCCOeH  4USMUTIHE - - 

TX  susPourrx  cok*pt<  iw^jt  Tire  function  to  tm  ctjTP-jr  iiNpav 

W>TA  STPEWt.  tCTU  INPUT  PNC  CUTWT  19  COX  THSC'JSM  PMPVS. 

iKttttttttttstttt  uPNiPitcs  txtxmxxststtmttttzaittxttntxan 
INPUT  •  AN  AARAV  COHTAIHINO  TX  INPUT  TINE  FUNCTION  8PNPIE9. 

OUTPUT  •  AN  ARRAY  CONTAIHINO  THE  OUTPUT  IlHARV  DATA  STJfEAN. 

N  •  THE  NUNKR  OF  SANPtXS. 

FCI,  FCa.  FC3  >  ROa-OFF  FREOUEMCIES  OF  THE  PRIMRV  INTEORATOd. 

TC  •  THE  TINE  CONSTJRfT  OF  THE  SVtLAIlC  FIITER, 

AlPHA  •  THE  DECAY  RATE  OF  THE  PRIPARV  INTEWATOd. 

lETA  .  THE  DECAY  RATE  OF  THE  SYUAIIC  FILTER. 

EH  •  THE  StCN  CP  TVE  OIFTEREHCE  DETUEEN  THE  CURRENT  INPUT  AND 
THE  CURRENT  ESTIMATE. 

ENt  •  THE  SION  OF  THE  DIFFERENCE  ONE  TX  PERIOD  ACO. 

CW  •  TX  SION  OF  TX  DIFFERENCE  TWO  TIX  PERIODS  AOO. 

UdlN  •  TX  NIHiniM  INPUT  TO  TX  SVtLADlC  FILTER. 

UMAX  •  TX  NAXINUn  INPUT  TO  TX  SYLLADIC  FILTER. 

FS  •  TX  SANPU  RATE. 

OELTAN  •  TX  aRRENT  STEP  SIZE. 

DIF  •  TX  DIFFERENCE  lETVIEEN  TX  CURRENT  IWVT  AND  TX  CURRENT 
CSTtnATE. 

Wl  •  TX  CURRENT  ESTIMATE. 

DC  •  TX  DUTY  CVCLE  OF  TX  SLOPE  OVERLOAD  DCTECTOR  FOR  TX 
CURRENT  INPUT  STRINO. 

Cnx<x<<nn*sstsxzxssxxxzz3»r3(sxxxnunzxxzxxxa»xzxaxxn>zxsxaxu 

C - 9UDR0UTIX  START - 

O—  INITIALIZE  UARIAILES  AND  MRRAYS 

REAL  INPUT(N) 

INTEOCR  ourpur(N) 

DATA  £»!/•>’,  CNS^.  PI/^.I4159ESS3C/ 

SUN  • 

DEITAN  •  UniN 

0——  OALCUIATI  DECAY  RATES  OF  ENCODER  FILTERS 

alpha  •  EXP  <•(«.  t  PI  «  FCI  /  FS)> 

KTA  •  EXP  (-(a.  <  PI  /  TC  FS)) 

C— r  fTART  CXOOIX 

M  M  I  •  t.N 

C - CAWULATt  TX  OUTPUT  OF  TX  CONPARATCR 

PIF»  INPUTd)  -  XN 
CN  «  SlCN(i..01F) 

C - CCXRAX  XW  ESTIAATE 

)«  ■  ALPHA  S  )«  ♦  (I  -  ALPHA)  •  OELTAN  KM 
U  >  WHIN 
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C -  CCICMTE  THE  NEXT  Ol/m/T  OF  T>«  StOfC  OUCJttOAO  OCTECTM 

If  («<(EH  EHl)  .AND.  ENS)  .EO.  1.)  .0«. 

1((<0«  .AND.  END  .AND.  E.n3>  .EQ.  -1.))  U  •  MMX 
If  (V  .EO.  MMX)  SUN  •  6Un  A  1 

C -  6ENERATE  rCXT  STEf  SIZE 

KLTAM  •  lETA  «  OELTAN  A  (1  -  tETA)  t  U 

0— ~  CHIfT  TME  SlOfE  OUERLOAO  SCTECTOA  SHIFT  REOISTEll 

EN3  •  ENl 
ENt  •  EN 
OUTPUTd)  •  EN 

e -  POUW  TO  IINARV  COfWERT 

If  (EH  .EO.  -1)  OUTPl/Td)  •  • 

M  CONTINUE 

0 -  CALCUUTE  S10(>E  OVERLOAD  DCTECTtM  DUTY  CYCLE 

DC  ■  sun  /  H 


RETURN 

END 
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APPEMDIX  C 

CVST)  Decoding  Subroutine 


fUMOUTiNi  Kcooci ( iMUT,oumrr,N.r«.rci,Fci.rc3.TC,umx.uniN,DCi 

- - -CWSO  OCCOOINO  SUIWOTINC - 

TMIi  SUBROUTINE  CtCOOES  TVE  BINRflV  6ATR  STRERN  CONTRIMT®  IN  T« 
INPUT  WRAV  AND  PUTS  T>C  OUTPUT  TINE  FUNCTION  SANPLES  IN  THE  OUT¬ 
PUT  ARRAY. 

CBBSttttttttuMtttsusni  uapiables  tutsnsnnstsxnsnxnxnsnnnx 
C  IrMT  •  AN  ARRAY  CONTAIHINO  THE  INPUT  BINARY  MTA  STRSAN. 

C  OUTPUT  •  AN  ARRAY  CONTAININO  THE  OUTPUT  TIRE  FUNCTION 
C  N  •  TNE  NURBER  OP  SARPUt 
e  PS  •  Ttc  SARPU  RATE 

C  PCI,  pca,  FC3  •  ROLL-OFF  FREQUENCIES  OF  THE  PRIRARV  INTEGRATOR 
C  TC  «  THE  TIRE  CONSTANT  OP  TIC  SYLLABIC  FILTER, 
e  KH  •  nc  CURRENT  OUTPUT  TIRE  SARPLE 
e  ENl  •  THE  SIGN  OF  THE  DIFFERENCE  ONE  TIRE  PERIOD  AGO 

e  oe  •  THE  SION  OP  T>C  DIFFERENCE  TVO  TIRE  PERIODS  AOO. 

e  UIMN  •  THE  NAXIHUR  IlfUT  TO  THE  SYLLABIC  FILTER 

C  URIN  •  the  HINIRUR  imjT  TO  THE  SVIUBIC  FILTER 

C  DELTAN  •  THE  CURRENT  STEP  5I2E 
e  ALPHA  •  THE  DECAY  RATE  OF  THE  PRIRARY  INTEGRATOR 
e  BETA  •  THE  DECAY  RATE  OF  THE  SVLLAJIC  FILTER 
C  DC  •  T»C  SLOPE  OUERLOAO  DETECTOR  DUTY  CVCU. 


C- —  GET  NEXT  IWT  BIT  AND  BIMRV  TO  POLAR  COtWCRT 
CN  «  INPUT! I) 

IF  IlNPUTdl  .EO.  •!  CN  •  -t 

e - GENERATE  NEXT  OUTPUT  TINE  SANPU 

XH  •  ALMA  t  XN  ♦  <1  -  ALPHA!  •  DELTAN  SEN 
U  *  uniN 
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C - CCNCMTC  THE  NEXT  OUmiT  OT  THE  tLon.  omulom  petectcw 

IF  {(<<  EH  .fiNO.  END  .AN9.  EH2)  .E9.  1.)  .0«. 

><((EN  .M<0.  EHtI  .rtxo.  E'3>  .EQ.  -l.)>  U  ■  MMM 
IF  (U  .EO.  UHAX)  SUfl  •  SUH  ♦  1 

C -  GCNEMTE  next  STEF  St2E 

OeiTAN  •  SETS  S  OEITAN  ♦  (I  >  SETA)  S  V 

C - SHIFT  T>C  SLOPE  OtXXLOAO  CrrECTOR  SHIFT  REaism 

ENE  •  ENl 
ENt  «  EN 
OUTPUTd)  •  XN 
M  CONTINUE 

e>—  CALCULATE  THE  SLOPE  OUEffLOAO  DETECTOR  DUTY  CVCU 

DC  •  sun  /  N 

RETURN 

END 
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FIR  Fiiterinp,  P>ubroutine  (FIITI’R) 


sutMurir*  FiLTOKirr,  n.  hp.  ii 

.  . -fl«  FILTER  SUWtOUTIHE— . . 

THIS  SUBROUTIME  FILTE*)?  iNPtJT  TIPS  FL'VCTICH  SP'T’LE  STPIHO  USIMO 
0  FILTER  COtFFICIEhTS  r:>-Z.!rATED  BV  TN  EXTT':‘JRL  FlLTia  Cir  CBRToa 

c  wutine.  tv-.e  fiut.:--  £d  ti:  i  fu'ctim  c-:z  piAcra  in  twe 

C  »<«E  R.TRRV  AS  Th-:  =tjT  F.O  B-TL'N  TO  TrZ  C.-:’.L3  ,3  f. An.  WJC 

C  TO  the  filter  T-CWVK-M.  C:3  SAri>US  NT  TK£  lEOINHINO  NNO  CHB 

0  OF  TKE  SANPLE  STRINQ  Ml  LOST. 

cttuttBttttBMxuuxnxsn  umiailes  sxxxnxxxxxxzxxxtxsxxxxxxttxxxxxx* 

C  XT  •  THE  ARRAV  CONTRTNINO  TVE  INPUT  S6f5>iE  smiNQ  mq  <»rTER  PRO- 
e  CCSSINO,  THE  FILTERRED  SAflPU  STRINO. 

C  N  •  THE  HUraSR  OF  SPfTLES  IN  THE  INPUT  PRRAV 

C  HP  -  THE  NUNBER  OF  FILTER  COEFFICIENTS 

C  ■  •  M*  ARRAY  COfTAININO  T>C  FILTER  CCEFFICIENTS 

CMMtxttmsxmxtxttmttxxtxtxxxxttxtxtTttxtmtxxxmtnxtaxtttxtaxi 

C - SUBROUTINE  START - 

C -  INITIALIZE  UARIABLES  AND  ARRAYS 

OINEKSION  XT(K).  BITT} 


FILTER  THE  IfSVT  SANPLE  STRINO 

00  IM  I  •  1.46M 
K  •  20S  ♦  I 
DO  6S  J  •  1,TS> 

IF  U  .EO.  IJ  SON  •  f<J>  t  XTtCI 

SUN  •  SUN  ♦  B(J)  t  (XT<K  ♦  4  -  1)  ♦  XTIK  -  4  ♦  1)» 

CONTirSiE 

XT(1)  •  SUN 

CONTINUE 

RETURN 

END 


Y 

It 
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r»r»e,r,.* 


APPENDIX  E 


FIR  Filter  Coefficient  Generating  Subroutine 


tUWOUTIfC  riTROCN(BrrA.OM«M.NP.>> 


— fMXtfMLlV  FLAT  FILTCR  PKOMM- 


TMlS  PROGRAH  OUTPUTS  T»«  PIR  FIITER  COXFFICIENTS  CALCUIRTED  IV 
iOlROUTlME  PXFUT.  TXIS  RCUT’rt  F.IT^:R  riTU:i‘*S  TVE  MLCUIATED 
COEFFICIEMTS  TO  TV*  CAtUIVQ  fRCCi/wi  oR  PRIMTS  OUT  TVE  ESSCR  FES- 
SaCES  UHEN  TVE  CCtFFICXENTS  CPT««T  l£  DETC»11N£D  SUE  TO  THE  CHOICE 
OF  IfRVr  PARPnETtRS. 

THIS  SU9R0UTINE  PMB  THE  fWFtPT  AHO  RATPRX  fWIROUTlHES  USED  TO 
CEMERATE  TI->E  PAXir^UV  FtAT  FIB  FIL-TB  IIIE.VTO  A^E  ADAPTEB 

FRon  A  PRccapn  riuELCPEX)  »v  j.  f.  tc.nciB  cf  eell  u^cjatcries. 

THIS  PROCRAd  (.‘AS  ’LIC'-£D  IN  ‘PROCRAdS  FOR  DIGITAL  8IC/mL  PRO- 
CESSIMG*  tv  THE  IEEE  PkESS. 


eitmttttsmuttmu  uariailes 

C  KTA  •  THE  NOAHALIEES  CCKTER  FREQUCiNCV  OF  TIC  TfMNSITIOH  IANS 
C  MfRM  •  THE  N0RIMLI2E0  UIOTH  OF  THE  TRANSITION  lAND 
C  M>  •  THE  NUHKR  OF  FILTER  COEFFICIENTS 
e  t  •  AN  ARRAV  CONTAINIHQ  THE  FILTER  COEFFICIENTS 
C  UNIT  -  THE  LARGEST  NUKKR  OF  FILTER  COEFFICIENTS  ALLOICS 
e  lERR  •  THE  HUntER  OF  THE  ERROR  NESSAGC 
e  A  4  C  <  WORKING  ARRAVS 


— SUSROOTIHE  START- 


It 

«S 

90R7 


INITIALIZE  UARIAILES  AND  ARRAVS 

ftlNENSlON  A<RMM(RM>,e(2SS) 

UNIT  •  S99 

CALL  RXFLATdCTA.  OATIM.  HF.  A,  I.  C,  liniT,  lERR) 
PRINT  RESULTS 

IF  (lERR  .or.  1)  URITE(fi.9SSS>  SETA.  GANNA 
FORNATC  FOR  SETA  •  ‘.FS.S.’  AND  OANNA  >  ‘.FS.I) 

CO  TO  (IS,  £S,  3S,  4S),  lERR 
RETURN 

URZTE(S.8977) 

FORIMTC  lETR  NOT  IN  RANGE  S.  -  .S*> 

STOP 

WRITECS.SOPS) 

FORNATC  GATIW  NOT  IN  RANGE*) 

STOP 

VRITEIS.SCSS) 

FOF^rc  GANNA  TOO  SNALL,  NIH  IS  .SST*) 

STOP 

END 
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OOAAO  A  A  AA  rtf,  ^  r,r. 


APPErroix  F 

Subroutine  nXFTAT  -  Part  of  FIR  Fi'l‘-.gr  Generator 


fUWOUTINC  nxnATdC.  M.  A.  •,  C,  LinXT,  XUOI) 

c — . . sirtnouTiNC  mrtAT - 

THIS  «uBffouTiNe  comntt  the  cocrrxexEMTt  or  a  flAxiMUv  flat  rx« 

lirCAA  PMASE  FILTER. 


tttttttttttnttMntUHI  UARIAOLCO 

K  •  CENTER  OF  THE  TRANSITION  REOICN.  RMIQE  •  A.  TO  .E 
FREQUOICV  IS  NORfMLIZEO  TO  THE  SMVU  RATE. 

OA  •  UIOTH  OF  THE  TRANSITION  RCOION,  WHERE  THE  OUTRUT  ANRIXTUK 
DECREASES  FRON  8SS  TO  6X. 

UNIT  -  THE  NAXIHUn  NUNIER  OF  COEFFICIENTS  IN  TIC  FILTER 

t  •  THE  ARRAY  CONTAININO  TIC  FILTER  COEFFICIENTS 

XERR  •  ERROR  RESSAOES 
mortal  >:£TU:Rf1 
e,  SETA  not  in  range 

3,  GAMMA  NOT  IN  RANGE 

4,  GAMMA  TOO  SMALL,  LESS  THAN  .*4 

C  A  ■  WORKING  ARRAY 

e  e  -  WORKING  ARRAY 

e  K  >  NUKSER  OF  ZEROS  AT  NYQUIST  FREOUCNCY 

e  L  •  NURSER  OF  ZERO  DERIUATIUES  AT  ZERO  FREO 

e  NT  «  FILTER  HALF  ORDER  ■  NP  -  I 


csnnssssxsttwsttTstitsmmmmsxtusmttswitnttstitiiCTmRt 

C - SUSROUTIMC  START - 

e —  initialize  wariaius  and  arrays 

DIMENSION  AILIMIT),  S(IIMIT),  CdlMlT) 
lERR  ■  1 

NP  •  t 

TUOPI  >1.1  ATANIl.Si 

IF  ((BE  .U.  S.)  .OR.  (IE  .OE.  .E)>  00  TO  SO 
m  •  AMINKR.  I  BE.  1.  -  3.  S  BE) 

IF  ((GA  .LC.  •.)  .CR.  (CA  .CC.  BM))  00  TO  SS 
NT  •  1NT(1.  /  (4.  t  GA  B  GAD 
IF  (NT  .OT.  1E8)  CO  TO  ICB 
AC  ■  (1.  T  CCSCrUOFI  S  BE))  /  t. 

«IIN  •  LIMIT 

GAU  RATP«m(AC,  NT,  K,  NP,  QUN) 

N  >  a  1  KP  -  1 
IF  (K  .EQ.  S)  K  >  I 

O— >  OOmTE  HAQNITUX  AT  NP  POINTS 

Cd)  -  1. 
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APPErmix  G 


?=!ubroutine  RATPRX  -  Part  of  FIR  Filter  Generator 


o  €>  o  «  o  o  <#  r 


appendix  h 
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SUI  •  sun  ♦  B€tT» 

M  CONTINUE 

C -  CAtCUWTE  DUTY  CVCU  «ffl)  MIRAGE  STEP  SIZE 

C 

CO(I-l)  •  t.  /  I 
STEP<I-1)  •  SUN  SS*. 

IM  CONTINUE 

C -  CAIXUtATE  T>C  NUERNGC  STEP  UMiJE  FON  N  DUTY  CYCtX  QT  ZERO. 

e 

ceiTN  •  t. 

SUH  ■ 

U  •  UNIN 

DO  EM  J  ■  l.SO*  _ 

DELTA  •  tETA  I  DELTA  ♦  (t.  -  SETA)  X  U 
SUN  •  sun  T  DELTA 
tM  CONTINUE 

CO(SS)  •  S.  „ 

STEPtMl  ■  SUN  /  SM. 

S -  MINT  AND  PLOT  T>C  RESULTS. 


£• 


ISM 


PRINT  ».  •  FOR  THE  FOttOUINO  SVSTEN  PARANETERSI* 

PRINT  t.*  5AWU  RATE  .  '.FS.*  IPS* 

PRINT  t.*  TC  •  '.TC 

PRINT  t,  •  FCl  •  ‘.FCl 

PRINT  I,  *  RATIO  •  ‘.nsTlO 

PRINT  *.*  THE  FOLLOUINO  SVSTEN  PaPANETERS  AREI* 

PRINT  *.*  WAX  •  ■.(.■'■U,*  UIIH  • 

PRINT  «.*  NAXIHUN  STEP  SIZE  •  *.STEPI1).*  HINIHUR  STEP  SIZE  •  *. 
ISTCPISS) 

IP  (IT  .CT.  1)  CO  TO  DM 
CALL  FACTOR! .S) 

CALL  PLOTO.,  2..  -3) 

CALL  SCAl£(CD,19.I.'i«.l) 

CALL  SCAL£(ST£P.6.J.Ji,J) 

CALL  AXlS(a.«,«.«.l>>rEP  SIZE  (V).I3,6.S.M.S.STEP(C1),STEP(53)) 
CALL  AXlSie.O.a.S.K^fwtiTV  CVCLE.'lS.lS.S.t.S.CDCSD.COCSa)) 

CALL  STECTIS..  9.,  IS.,  3) 

CONTINUE 
ICNIW  •  IT  -  1 

CALL  LINECCD.STEP.SS.l.lS.ICHAR) 

CONTINUE 
CALL  PLOTE(N) 

END 
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O  O  O  O  O  O  OO  Q  o  a  vo  o  o  ov 


APPErtDIX  I 


WTAX  Calculatinn  ^^ubroutine 


fuwouTirc  unM(orr(uitAx.uniN.M.rci>Te.MTio) 


C  TMti  iUS90UTlKe  CftLCULATES  TH?  W»U.t  Of  UBAX  AMU  UfllN  RA-’EO  OW 
•'  TM6  IhOUT  S«.'-LE  RATE.  SVUADIC  rifm  TC.  F7: ’‘.HV  INTEC  ATOR 

BOLL-ORR  FREE  -^CV  (RCl).  i.  J  T)-£  ,  ',TIO  6£Tv:..X  TH'  P'Xir.Ufl  STEF 
SIZE  HIM  .1  S*'=  S:ZE  CUT?'JT  C-  T;  E  tvi,'.  JJC  FILT"R.  THIS 
CALCUUATICM  ’5  fZ  -'C-.-rS  AT  A  !.  i.-T-tr'Oi  f'-iZO'  .MV  C*'  ’  3  HZ  A«0 
SICNAt  R'<,=LJT;:tS  C?  0  TFi  CALU^S  CF  t,  J<  A^0  U'lH  A^JE 

CALCULATED  SUCH  1>WT  THE  DUTY  CYCLE  OF  TkS  6LC^  OVERLOAD  DETECTOR 
OUTPUT  IS  .S. 

ttXXttttttttXIXXttXXXXXXXX  UARIAILES  txxxxxxxxxxxxxxxxxxsxxxxtxxxxxxxxx 
UtMX  •  THE  HAXinun  lr«>UT  TO  THE  SYLLAIIC  FILTER 

uniN  •  T>c  NiHiNun  im/T  to  the  syllabic  filter 

FS  •  THE  SARPLC  RATE 

FCt  •  THE  ROLL-OFF  FREOLOCY  OF  THE  PRIRARV  IHTEQRATOK 

TC  •  THE  TIRE  COHSTANT  OF  THE  SYLLABIC  FILTER 

RATIO  •  THE  RATIO  rSTVEEN  THE  nAXlFSJA  STEP  SITE  AND  TIC  HIHWUH 
STEP  SIZE  OUTPUT  OF  THE  SYLLABIC  FILTER 

TS  •  AN  ARRAY  CONTAININO  THE  TEST  SIGNAL  SANPLES 

BINOUr  •  AN  MRAY  CONTAININO  THE  BINARY  OUTPUT  OF  TIC  ENCODER 

FEAKl  .  THE  PEAK  UALUE  OF  THE  TEST  SICNAL  ANPLITUDE 

DIF  •  THE  DIFFERENCE  tSTLEEN  THE  SLOPE  OVERLOAD  DETECTOR  DOTY 
CYCLE  AND  THE  DESIRED  VALUE  OF  .5. 

RAT  •  THE  RATIO  CF  THE  DUTY  CYCLE  OIFFEREHCE  TO  THE  DESIRED  VALUE 

tnnmmstnxnns  subroutikes  used  xxxxxxxxttxxxxxtxxxxtxxxxxtxxn 


SIGNAL  •  TIC  TEST  SICNAL  GENERATOR 

UniNOPT  •  CALCL"LATES  the  value  of  URIN  that  pairs  with  the  CAL¬ 
CULATED  VALUE  OF  VNAX 

C  ENCOOEl  ■  TIC  CVSD  ENCODER 


0 . . -  SUBROUTINE  START - 

0 -  initialize  arrays  and  variables 

DINENSION  TS(409S) 

INTECER  IINOUTUCSSI 

A(0CN9)  •  S;»<T(It.  SSKOSNS  -A.  VIS.)  t  .#•!  S  SM.  >  B  SQRTtl. ) 
PEAK!  •  A(S.) 

VHAN  •  IS. 

e— —  CALCULATE  TEST  SIGNAL  SAMPLES 


C- 


CALL  S10NAL(TS,4S9S.FS,StS..S..PCAKI,S.> 

—  START  MAX  CALCULATION  LOOP 
CONTINUE 

—  CALCULATE  ESTINATEO  ENCODER  PARANCTERS 
CALL  VNINOPTCVNAX.vniN.FS.TC.RATlO) 

—  PROCESS  THE  TEST  SlGTMl 

CALL  CNCODEI(TS,IINOUT.ASOS,FS,FCI,FC8,FC3,TC,Vmx.UNIN.DC » 

—  FIND  TIC  DlFFEPtNCE  lETUEEN  THE  DUTY  CYCLE  USING  THE  ESTIRATED 
MAX  AND  VrtIH  Fuh)  THE  LESIkED  DUTY  CYCLE 

DIF  .  DC  -  .5 
RAT  •  DIF  /  .B 
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VO 


IF  THE  DOTY  CYCLE  IS  WITHIN  lx  OF  THE  DESIRED  UfttUE  KTUHN 
THE  UTMX  AND  WHIN  L>AU£8  TO  THE  CAUINO  PHOCIMH 

IF  (AtSdtAT)  .LE.  .•!>  00  TO  94* 

OTHERWISE  REESTinATE  VmAX  AMD  REFEAT  CALCULATIONS 

UNAX  ■  UNAX  ♦  .S  S  RAT  t  UNAX 
CO  TO  5 
MC  CONTINUE 
RETURN 
END 
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Of»  A  O  on  OA  A  A  A  A  A  AAA 


APPENDIX  J 


Ca]  culalj  ir?  SnrrouMn*^ 


3UIMUTINC  iiniNOPTCUrMX.uniN.FS.TC.IMTlO) 

ktttitttttttitttttttttttttntMtstnutnutxmtxxmnsnuxntnsstsx 


this  SUIHOUTIHE  COLCUtATES  TKS  UAUt  OF  WIN  THAT  PAIFS  UITH 
th€  i^lue  of  twat  i-j  u-m  so  that  tke  ratio  cf  tve  nah- 
inufi  STEP  SIZE  TO  niN:r(.ii  srr?  site  at  t,  -:  cutaut  of  tvsE  svilabic 

FILTER  IS  UIThIN  .OIX  Of  THE  VALUE  EAtCIFIED. 

ut«t«««tx*uu«*»xx«(  VARIABLES  ttxxxxttttxssmtntttmssxtxttta 
VTMX  •  THE  nAXlflUn  IHfVr  VALUE  OF  1>C  SYLLABIC  FILTER 
UniH  •  THE  niNinUH  value  input  to  tic  SYLLABIC  FILTER 
FS  ■  THE  SAHPLE  RATE 

TC  •  THE  TIRE  COMSTAMT  OF  THE  SYLLABIC  FILTER 

HAXSTEP  •  THE  mXinUR  STEP  SIZE  AT  THE  OUTPUT  OF  THE  SYLLABIC  FIL¬ 
TER 

niNSTEP  •  THE  RINinun  STEP  SIZE  AT  THE  OUTPUT  OF  THE  SYLLABIC  FIL¬ 
TER 

BETA  •  THE  DECAY  RATE  OF  THE  SYLLABIC  FILTER 
SUN  «  TIC  RUHNINO  sun  OF  STEP  SIZES 

RATIO  •  THE  fESIRED  RATIO  SErLTEN  THE  PAXimn  STEP  SIZE  AfO  THE 

niNinun  step  size  at  the  output  cf  the  syllabic  filter  ih  db 

e  R  •  THE  U0L7AQE  RATIO  EOUIVALENT  OF  RATIO 
C  DELTA  •  THE  CIWREHT  STEP  SIZE 

CttusxmtttxxxttxttxxmtxxmxxzmttxmtxxxmtxxxtttUMXtxxxtttxa 

e— -  INITIALIZE  VARIABUS  AHD  ARRAYS 

REAL  NAXSTEP,  PINSTEP 
DATA  PI/3.HlS":S3o/ 

R  •  IB.  SX  (RATIO  /  £4.1 

—  CALCULATE  SYLLABIC  FILTER  DECAY  RATE 

BETA  •  (-(3.  S  PI  /  TC  H  FS)) 

e -  ESTINATE  INITIAL  VALUE  OF  MIIN 

umn  ■  MMX  /  IM. 

e  START  CALCULATION  LOOP 

B  CONTINUE 

e -  XHIHAL  RMMINQ  SUN  AND  STEP  SIZE 


SUN  ■  S. 
BELTS  •  B. 


0— . 


IB 


CALCULATE  AUERACE  NAXINUN  STEP  flIZC 


DO  l«  I  •  I.5M 
IDT  •  1  3 

IT  •  IDT  (  3 
V  •  VNIN 

IF  (IT  .£0.  1)  V  •  VNAX 
DELTA  •  BETA  «  DELTA  ♦  (t 
SUN  •  SUN  ♦  DELTA 

continue 

NAXSTEP  •  SUN  H  SBt. 


BETA)  B  V 


C 


REINITIALIZE  RUNNIHO  SMI 
DELTA  •  f. 

SUN  •  B. 


AND  STEP  BXZt 
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1 


C' 


CALCUWTt  CURROfT  ESTWTt  OT  THE  HINIHUH  rTER  «I2£ 


IS 

c- 


c - 

c 

999 


DO  15  I  •  1.5:« 

DEITR  -  8£T»  t  DE1.TR  ♦  (1.  -  tETR)  t  WITH 

sun  •  SUM  ♦  delta 

CONTINUE 

HINSTEP  •  sun  /  5M. 

FINS  THE  DIFFERENCE  KTUEEN  THE  ESTIMATE  AND  T>C  SPECIFIES  RATIO 

TNIN  •  NAXSTEP  /  » 

OIF  •  THIN  -  niNSTEP 

RAT  •  ASS  (OIF)  /  TNIH  S  IM. 

IF  THE  DIFFERENCE  IS  LESS  THAN  .AIX  RETURN  UNIN 

IF  (RAT  .LE.  .tl)  00  TO  999 

IF  THE  DIFFERENCE  IS  CREATER.  THEN  REESTINATE  W1IN  AND  REPEAT 
CALCULATIONS 

UNIN  •  UNIN  «  .S  S  OIF 

CO  TO  S 

CONTINUE 

RETURN 

END 


I 


( 

! 


1 
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APPPJNDIX  K 


CVSP  't'T.I.o-'i  Si,!'-' 


mOCRWi  «n.SttlN«rr.OUTWT.TAPC€*OUT«JT,^LOT) 

C  this  PROCPAB  plots  TuE  OITTPUT  SIGNAL  OP  TNP  CUSO  TPAK5HISSI0N 
C  SVSTEfl  L-HtH  THE  I-'OLIT  S:r.“AL  IS  A  Z79  HERTZ  SINE  UAvE 

C  that  OAPIES  15  ft-»L:T'..:E  rPf'fl  -■'3  Z'-'i  TO  8  I ' -0.  TEST 

C  SIGNAL  GE'tEPATOR  ALTEEha-'-lv  SCO  fA^LES  AT  -aR  O'!"# 

C  AND  0  Dr“0  IN  CRIER  TO  IT  -N-STRATr  y„-:  sVjTE'''C  STEP  PESPONSE 
c  Characteristics,  th:  sv-j-in  unc;.r  test  consists  cp  -ve  input 

C  PILTEH.THE  CU3D  ENCCIEH  AND  DECOIER.  AND  Th£  OUTPUT  FILTER. 

UARIAILCS  XtttttXXtttXtttttttXXXXXtXXXXttXtXXX 

C  TSIN  •  AN  ARRAV  CONTAINING  THE  INPUT  TINE  SERIES  SAfVLES. 

C  TSOUT  •  AN  ARRAV  CONTAINING  FIRST  THE  DECCDER  OUTPUT  TI"E  SERIES 
C  SAHPIES,  THEN  THE  OUTPUT  TINE  SERIES  SANPLES  OF  THC  FIR  FILTER. 

C  •  •  AN  ARRAV  CONTAINING  THE  FILTER  COEFFICIENTS. 

C  TINE  •  AN  ARRAV  CONTAINING  THE  TINE  THAT  THE  FIRST  20*  SANPLES 
C  ARE  TAKEN  SO  THAT  THEV  BAY  8E  PLOTTED. 

C  IINOUT  •  AN  ARRAY  CONTAINING  THE  BINARY  OUTPUT  OF  TIC  CUSD  ENCODER 

C  AflPl  •  THE  ANPLITUSE  OF  THE  TEST  SIGNAL  IN  DBA*. 

C  FS  •  THE  SAMPLE  RATE. 

C  FCl,  FCa.  FC3  •  THE  ROLL-OFF  FREQUENCIES  OF  THE  PRINARV  INTECRA- 

C  TORS. 

C  TC  •  THE  TINE  CONSTANT  OF  THE  SYLLABIC  FILTERS. 

C  UMAX  t  UNIN  •  THE  NAXINUN  AND  NININUN  INPUTS  TO  THE  SYLLAIIC  FIl- 

C  TER. 

C  tCTA  •  THE  N0RBALI2ED  CENTER  OF  THE  TRANSITION  BAND  OF  THE  LOW 
e  PASS  FILTER. 

C  GANNA  •  THE  NORBALISEO  WIDTH  OF  THE  ROLL-OFF  REOICN  CF  THE  OUTPUT 
C  FILTER.  THE  KEGICN  IS  THE  FREOOENCV  BAND  BETWEEN  THE  95B  AND 

C  SX  OUTPUT  ABPLITUDES. 

C  PEAKl  •  Tie  NAXINUN  ANPLITUDE  OF  THE  TEST  SIGNAL  IN  WOLTS. 

C  NP  •  THE  NUNBER  OF  FILTER  COEFFICIENTS. 

C  DC  •  THE  duty  CYCLE  OF  THE  SLOPE  OVERLOAD  DETECTOR. 

c  Ratio  •  the  ratio  between  the  naxinun  step  size  and  tic  nininun 

C  STEP  SIZE  IN  DB 

Cxxttxttttttttttttittxxxxttxxxtixxxxxxxxxtxxxxxxtxxttxxxxxxxxxxxxxxxtxxx 
Cxxtxtxxtttxttsxxtxtx  SUSROUTIHES  USED  tXZSlX»t(XIUS»X»SZa»(l»XX*» 

C  FLTRGCN  •  THE  SUBROUTINE  THAT  GENERATES  THE  OUTPUT  FILTER  COEFFI- 
C  CIENTS. 

C  FACTOR,  PLOT,  AXIS,  SCAU.  RECT,  LINE,  PLOTE  •  CALCONP  PIOTTINO 
C  ROUTINES 

C  SIGNAL  •  THE  TEST  SIGNAL  GENERATOR.  PRODUCES  SANPLES  OF  SINUSOI- 
e  DAL  WAVES  WITH  AT  HOST  TWO  FREQUENCY  COBPOtCNTS, 

C  CHCOOEl  •  THE  CVSO  ENCODER  SUBROUTINE  WITH  A  SINGLE  ROLL-OFF 
C  FREQUENCY  IN  THE  PRIMARY  INTEGRATOR. 

C  OCCODEl  •  THE  CVSO  DECODING  SUBROUTINE  WITH  A  SINGLE  ROLL-OFF 
C  FREQUENCY  IN  THE  PRINARV  INTEGRATOR. 

C  FILTER  •  THE  SUBROUTINE  THAT  FILTERS  THE  INPUT  TINE  SERIES  SANPUS 
C  USING  THE  FILTER  COEFFICIENTS  GENERATED  BY  FLTRGEN. 

C  PLTIBE  •  A  LINEAR  PLOTTING  ROUTINE  TO  PLOT  SIGNAL  ANPLITUDE  US. 

C  TINE. 

C  UHAXOPT  •  GENERATES  WAX  AND  WIN  FOR  THE  CUSD  ENCODER  AND  DECODER 
C  SUBRUTIHES 

CXSSXXXxtstxtxxxxxxxxxxxxxtxxxtxxxxxxxxxxxxxxxxxxxxxtxxFxxrxxxxxzxxrxttx 

cxxxxxxxtxxxxxxxxxxxxxxxxxxxxxzxxxxxxxxxxixxxxxxtixxxxxxuxxxxxxxxxxxxxx 
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PROCIHW  STAITT' 


C -  INITIALIZE  UARIAIUS  AND  AmAVS 

DinCNSIOH  TSIN(Sa««),TS0UT(SM*>,rtnE(Z*M),l(2M> 

INTEGER  IINCUT(Scet) 

A(Din«)  •  S3;^T(ltt.  -4.  Vi*.)  <  .Mi  S  CM.)  S  CQIITCa.) 

ANPl  •  -42. 

ANPZ  •  •. 

PEAKl  •  A  (AAPl) 

PEAKa  •  A  (ANPa) 

KN  •  t. 

C -  INPUT  ANO  PRINT  T>C  WORKINO  UARIAtLU 

READ  I,  PS 

READ  *,PC1,  TC.  PATIO 

READ  (.SETA.  GA.iriA 

PRINT  »,  ■  TRANSIENT  RESPONSE  TEST  AT  ‘.PS.*  BPS* 

PRINT  t,  ■  WITH  TC  ■  •.TC."  AND  RATIO  •  *, RATIO 

PRINT  filter  PARANETERS  ARE.  BETA  •  *,BETA,*.  OANNA  •  ‘.OMSM 

0 - GENERATE  FILTER  COEFFICIENTS  AND  CUSD  SVSTEJ1  PARAAETER8 

CALL  FLTROENIBETA, GANNA. NP.B) 

CALL  unAXOPT(0«AX,uniH,FS.FCi,TC.RATIO> 

C -  INITIAL2E  PLOTTER 

CALL  FACTOR (.5) 

CALL  PLorca..  a.,  -3) 

C - GENERATE  INPUT  TINE  FUNCTION  SANPUS 

CALL  SICNAL2(TSIN.6SM.FS.feSS..S..PEAKl.PEAK2> 

C -  FILTER  THE  INPUT 

CALL  FILTERITSIN.SSM.NP.l) 

C -  PROCESS  THE  INPUT  TIRE  SERIES  THROUGH  Tl*  CVSO  SVSTER 

CALL  ENC0DEl<TSIN,BIN0UT.Sa88,FS.FCl,FC2.FC3.TC.UNAN,umN.I>C) 

CALL  DCCOO£l(31NOUT.TSOUT.5S9«,FS.FCl.FC2.FC3.TC.GNAX.uniN,I>C> 

C -  FILTER  THE  OUTPUT  OF  THE  DECOCER 

CALL  FILTER<TS0UT,46M,NP.t) 

C -  PLOT  THE  OUTPUT  SIGNAL 


yv  »  »  - 

TIRE(I)  •  1  /  FS 
S  CONTINUE 

DO  6  I  •  1.225 
IK  •  2SSS  *  I 
TSOUT(I)  •  TSOUT(IK) 

C  CONTINUE 

CALL  PLTINE(TINe.TS0UT.22S.227) 

CALL  PLOTS (N) 

END 

cnttsttts»<t((>(»x»sstssHntnssssssxuKtns»nsn»nsnusnxn 

CsssMXXxxxxxxxxxxxxxixxxxxxxxxxxxxxxxuxxsxtxxsxxxxxxnxsxxuMusMxn 

SUBROUTINE  PLTINEIX.  V.  N.  NX) 


t . . TIRE  US.  AflPLITUOE  PLOTTER - 

C  THIS  SUBROUTINE  RAKES  A  LINEAR  PLOT  OF  TIRE  US.  ARPLITUDE. 
csssssxxxxtxxstttxxtxxxsxH  VARIABLES 
C  X  •  THE  ARRAY  CONTAINING  THE  ORDINATE  VALUES 

C  V  •  THE  ARRAY  CONTAINING  THE  ABSCISSA  VALUES 

C  N  •  THE  NURtER  OF  VALUES  IN  THE  X  AND  V  ARRAYS 

C  NX  <  N  «  3 

Ctututtttttsututsts  SUBROUTINES  USED  mtmxmttsumxsttttxtxxxx 
C  SCALE.  AXIS.RECT.  PLOT,  PLOTC,  LINE  •  CALCORP  PLOTTING  ROUTINES 

cnssstsstsstsssnstts»stssttsssnsstntsstnsiasnHsssx»ssBsssmus 
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c - suBROuTiKt  rrmr - 

C -  IHITtALIZC  UAfilAILES  AND  ARRAYS 

DIHENSIOtl  X(HX).  V(NX) 

C - SCAU  THE  X  AND  V  ARRAYS 

CALL  SCALEtX.  IS..  N.  1) 

CALL  SCALE(Y.  6..  N.  1) 

C -  MX  IN  THE  PLOT 

CALL  RECT<S..  S..  IS..  •..  3) 

C - ORAU  THE  AXES 

CALL  AXIS(S..  S..  ISHTINC  ISEC).  -IS.  IS..  S..  X(N4t >.  XCN^B)) 

CALL  AXISCS..  S..  13MAnPLlTUCC  CU).  13.  S..  BS..  Y(H«1>.  V(N«8)) 

C -  PLOT  THE  POINTS 

CALL  L1NE(X.  V.  N.  1.  S.  t! 

RETURN 

END 

etttxtttttttttxttttttt*tmtnn*nrmTttnxmtrtstnttnma**ntstt 

CMxxttxxxxxtxxxzitxixxxxxxxixxixxxtxtitxxrstxtxxtxsxxttixttxtsxsxxattat 


iUlROUTINE  SlCNAL2(0urPUT.N.FS.PRE01.PRE03.AAPl.ANPa) 

ixxxttxtxxxtitxxxxxxxxxxxxxxxxxxxxxxxxxxxxxsxxxxxxxxxxxxxxxttxxxxxxxxxxt 

C  THIS  SUIROUTINE  GENERATES  A  SINGLE  FREOUENCY  SINUSOIDAL  SIGNAL 
C  THAT  ALTERNATELY  H*?  SSS  SANPLES  AT  ONE  ANPLITUDE  THEN  5SS  SAMPLES 
C  AT  A  SECOND  ANPLITUSE. 

cs*t«s<tt«»»sst<t<t  UARIABLES  tsstuiutsxsxmnststurMtsuntm 

C  OUTPUT  .  THE  ARRAY  CONTAININO  THE  OUTPUT  TINE  FUNCTION  SAMPLU 
C  H  •  THE  HURIER  Of  SANPUS  TO  K  PRODUCED 
C  F$  •  THE  SANPLE  RATE 

C  FREOl  •  THE  FREQUENCY  OF  THE  TEST  SIGNAL 
e  FREOZ  •  UNUSED 

C  ANPi  •  THE  AMPLITUDE  OF  THE  FIRST  SSS  SAITLES 

C  ANPZ  •  THE  ANPLITUCE  OF  THE  SECOND  SSS  SANPUS 

C  ANP  •  THE  SIGNAL  AMPLITUDE  CURRENTLY  KINO  USED 

cxssnnsKssxttssssxsxMxttxsistssnnstnxumzsusttsnnsxnsntns 

c - SUDROUTINC  START - 

C -  INITIALIZE  UARIABLES  AND  ARRAYS 

DIREHSION  OUTPUT(H) 

DATA  P1/3.141S93S53S/ 

K  >  S 

S  CONTINUE 

e -  KT  CURRENT  SIGNAL  AfPLITUDC 


ANP  •  ANPI 
AHPl  •  ANPa 
ANPa  «  ANP 
J  •  s 

IS  CONTINU 


C -  GENERATE  SSS  TINE  FUNCTION  SANPUS 

J  •  J  *  t 

IP  (J  .QT.  5SSI  M  TO  f 
K  •  R  A  1 


c - If  N  SANPUS  HAUt  SEEN  GCfCRATED.  STOP  PROGRAM 


sn 


IF  (K  .GT.  N I  CO  TO  MS 
OUTPUTtK)  •  AMP  X  SIN  (3. 
CO  TO  IS 

continue 

RETURN 

END 


X  PI  X  FREOl  /  FS  I  K) 
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APPhUnix  i> 


I  d  le  Qian  no  1  ’''oiso  Pror'i'an 


MKMmW  IDLEHOI<IMPUT.OinWT.TA««.OtrrPlfr,in.OT> 

: - IDL£  CHWtCt  HOlSt  MO«(W - 

i  this  P«0<1R/«  measures  the  idle  OlAHWt.  NOISE  Of  THE  CUSO  TRAHS- 
C  WSSIOM  SVSTEN  U  iN  THE  LNCCl 'R  ON3  t£CC?: »  ft't  CC»«ECTC5  lACX- 
C  TO-IrtCJC  AHD  The  INPUT  TO  TIiE  ViCCiVt  IS  CiOUNEtO. 

C  THE  SVSTEN  SAIN  IS  AOJUSTEO  SO  THAT  AH  PC®  H2  INPUT  SISMftt  AT 

c  -ae  CB»«  PRcxiCES  a  -ea  om  signal  at  TrC  out>>ut  of  t>€  de- 

C  CODER. 

c 

CMUutttttistmttssstxx  UARiABLES  ttmmttmoosauiitiTiiimm 

C  FREOl  •  THE  REFERENCE  FKEOUENCV  USES  TO  SET  T>C  SYSTEM  GAIN. 

C  TSIN  •  AN  ARRAY  CONTAINING  THE  It«>UT  TINE  FUNCTION  SAMPLES. 

C  TSOUT  •  AN  ARRAY  CONTAINING  FIRST  THE  rECOOER  OUTPUT  TlfS  FllNCTlOH 
C  SAMPLES.  THEN  THE  OUTPUT  TIKE  FUNCTION  SAMPLES  OF  -ftC  FIR  FIL- 

C  TER 

C  B  •  AM  ARRAY  CONTAINING  TtC  FILTER  COEFFICIENTS. 

C  SINOUT  •  AN  ARRAY  CONTAINING  THE  BINARY  OUTPUT  OF  THE  CUSS  ENCODER 
C  ANPl  •  THE  AflPLlTUOE  OF  THE  REFERENCE  SIGNAL  IN  DMW. 

C  FS  •  THE  SAMPLE  RATE. 

C  FCl.  Fca.  FC3  •  THE  ROLL-OFF  FREGUEHCIES  OF  THE  PRIMARY  INTEORA- 

C  TORS. 

«  TC  •  TIC  TIME  CONSTANT  OF  THE  SYLIAIIC  FILTERS. 

«  UMAX  1  UniN  •  THE  MAXIMUM  AMS  MINIMUM  INPUTS  TO  T>C  SYUASIC  FIL- 

«  TER. 

C  SETA  •  THE  NORMALIZED  3  Ci  FFEOUENCV  OF  THE  OUTPUT  FILTER.  TIC 
«  FREQUENCY  IS  NC.V.ALIZES  TO  THE  SAMPU  RATE. 

C  ORPIA  •  THE  NOTMALIZa)  WIDTH  OF  THE  ROLL-OFF  REOION  OF  TIC  OUTWT 
C  FILTER.  T>'  S  -  :;IC.N  IS  THE  FTSOUENCV  SAND  BETWEEN  TIC  65*  AI® 

C  5*  OUTPUT  ««»LITUD£8, 

C  PEAK!  •  THE  MAXINUM  AMPLITUDE  OF  THE  TEST  SIGNAL  IN  UOLTS. 

C  NP  •  TIC  UfOER  OF  FILTER  COEFFICIENTS. 

C  DC  •  nc  DUTY  CYCLE  OF  THE  SLOPE  OUERIOAD  DETECTOR, 
e  PIN  •  TIC  POWER  OF  TIC  tlffUT  SIGNAL  IN  DBMS. 

C  POUT  •  THE  POWER  OF  TIC  OUTPUT  SIGNAL  IN  DBMS. 

C  ICN  •  TIC  CALCVAUITCD  idle  CHANNEL  NOISE  IN  DSNS. 

?  GAIN  •  TIC  UOLTAQE  AMPLIFICATION  OF  TTC  SYSTEM. 


f  tTTT’^TT^TT  ■»*»***"«***«”»*«”»*”*  ***«Fm*«qg 

SUSaOUTIICS  USED 


FLTRGEN  •  THE  SUBROUTINE  THAT  GENERATES  TIC  OUTPUT  FILTER  COEFFI¬ 
CIENTS. 


C  SIGNAL  •  THE  TEST  SIGNAL  QEICRATOR 
C  ENCOOEl  •  THE  CU9D  ENCODER 

C  DECODE!  •  THE  CUSD  DECODER 

C  FILTER  •  THE  SUSROUTIIC  THAT  FILTERS  THE  IlWr  TIME  FltCTION  SAM- 
C  PUS  USING  TIC  FILTER  CCEFFICIENTS  GENERATED  lY  FLTRGEN. 

e  POWER  •  A  ROUTIIC  TO  CALCLXATt  THE  POWER  IN  A  SAMPUD  TIME  Fl»«- 
e  TION  WITH  IRPEOCNCE  ■  600  OlflS. 


etUMStttstxsxuti^*t2iiiau'itxfi*tfii******ui*i****i**n*nf****f**^ 
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39 

C— 


—  INITIALIZE  UAAIMLE3  AND  ARMVS 

DINENSION  TSINCSMO.  TS0UT(6M«>,  XSM) 

REAL  ICN 

INTECCR  •IN0UT(5«««) 

ACDtiNA)  •  SeRT(l«.  tXtCDSNR  -4.  VI*. )  t  .**1  t  *•*. )  «  9QRT(3.) 

—  INPUT  AND  PRINT  UQRXINO  UMIAtLES 

READ  t.  PREOl.  PflPl,  PS 
READ  t.PCl,  TC.  r^TIO 
READ  t.SETA.  C.<^.TA 

PRINT  ».  •  ISLE  awf*CL  NOISE  TEST  AT  *.«.•  »PS' 

PRINT  X,  •  UITH  TC  •  •.TC,*  ‘’C  •  <TI0  •  '.Pr^TIO 
PRINT  t,  •  OUTPUT  FILTER  PARAWTEflS  ASEt  SETA  •  •.SETA 
WINT  t,  •  6iV«A  •  •,6««A 

OCHERATC  nC  FILTER  COEFFICIENTS  AND  CVSD  SYSTEN  P«N«TER* 

CALL  FlT^KEN<t^TA.CA!<f^A.^P,i» 

CALL  WtAXCPT««tAJ<,WIIN.fS,FCl.TC.RATIO) 

—  CCMERATE  IlWr  TINE  FUNCTION  SMTLES 
PEAKS  •  A  (A,-*?!) 

CALL  SICNAL<TS1N.S««*.FS.FRE01.*.,P£AK1,*.) 

—  PROCESS  TW  irsvr  tiw  function  tnrooo*  t>«  cuso  svsteh 

CALL  CNC00Cl<T5IH.»INi>fr.5AA*.F*.FCl,FCS,FC3.TC,lPV>X.ljniN.DC) 
CALL  OECODCl«IINOUT.TSOUT.E6««.F*.FCl.Fca,FC3,TC.UOAX,UNW,DC> 

"  FILTER  the  OUTPUT  OF  THE  DECODER 

•'ALL  FILTCRaSOUT.S«**,NP,B) 

—  DEUV  THE  IWVr  SIGNAL  START  TO  CORRESPOND  TO  T»C 
FILTERREO  OUTPUT. 

DO  3*  ID  •  S,4«0S 
KD  <  2C«  4  ID 
TSIN(IO)  •  TSINKKSS 
CONTINUE 

~  CALCULATE  TW  REFERENCE  fVSTEN  GAIN 

CALL  POLCR(TSIN,40««.FS.PIN) 

CALL  POL'ea<TSOUT,4;:3,FS,POUT) 

MIN  •  SWT  (PIH/POUT) 

—  GENERATE  A  ZERO  INPUT  SIGNAL  ARRAY 

00  4S  I  •  1.5*** 

TSIN(I)  '  *. 

CONTINUE 

—  PROCESS  T»C  ZERO  SIGNAL  THROUGH  T)C  SYSTEN 

CALL  CNC05El<TSIH.BINCUT,S*A*.FS.FCl.FCa.FC3.TC.VRW«,<JNlN.X) 
CALL  0£COOei(BINOUT.rSOUT.S«4*.FS.FCl.FC3,FC3.TC,LnAX.CNIN,K> 

—  FILTER  Tie  OUTPUT  SIGNAL 
OAU  FILTER(TSOUr,SM«.W.B> 

—  ADJUST  THE  OUTPUT  SIGNAL  AMLITUSC  TO  THE  REFERENCE  WPUC 
X  4S  I  •  1,4B96 

rsourci)  •  rsourdi  t  gain 

CONTINUE 

->  CALCULATE  THE  IDLE  CHAMCL  NOISE 

OAU  POUERITSOUT,4*9«.FS.POUT} 

ICN  ■  IG.  I  ALOGt*(POUT) 

PRINT  OUT  THE  RESULTS 

URITC<«,6*«)  ICN 

FORNATdXf'THE  IDLE  OWCL  NOIK  •  *,Pf.l) 

END 


115 


r 


OO  AOO  A  OA  A  A  A  A  AA  A  AAA  AAA 


APPWiDIX  M 


Total  liarnonic  Pi!'tortion  Pro'Tan 


MMOIST(l»»UT.OUTHJT.TAPtS.l(»UT.TMU»OUTPUT) 

. . tot^l  MflsnoNtc  oiSToarioH  p*»«(w - 


CALCUtATtS  TV^  TOTAL  HAPf^OHIC  DT^TCiTTlCN  TN  TWf  AfTPilT 
UMCN  A  single  Ff-!-.  ;ncv  TV'T  St-xAL  13  r4cc|f'"DT/-^~  ■-V  A^?n^ 
w  ENCODER  AND  CECCC.  ;<  c :cT£o  v»<:-or: 

'j  *^1^***  ONLY  TuCrC  t  -£CT  iL  Cf  ?C'-.NTS  Cf  TV'S  CLri>^T  THAT  L'£ 
t  1^0  m2  iViO  4l'v3  H2a  TrC-A  PT  CX«:%CTLV  m2  Cil  K2  m5t 

C  INCLutEO  IN  TNE  CALCULATION. 

cnttttttttMtttxsttxM  UAfiiAaLC«  tnxntxssxsnxsstxmxmssxstssam 


c 

c 

e 

e 


c 

c 

c 

c 

c 

c 


*'*^1*1!.*'"'“"  COMTAININO  T*C  IlNARV  OUTPUT  OP  T>€  CUSB 

^*■*^1,*..?  CONTAINIHO  Tt»  TITS  FUNCTION  OUTPUT  CP  TVC 

OUTPUT*OF''TH£^£Mlifl.^“  PaOCESSIW,  THE  TirZ  FUNCTlS 

’^’‘oF  ePTCTtfAL  PO-ER  CC«W«WT* 

PWCESSINC  i-V  THE  FAST  FCUJIER 

luc.  u(.  cuK  •  uomciHO  arrays  usa  tv  the  fft  sujroutdc. 

PRCOUE  •  A  REAL  ARRAY  CCNTAINIKO  TH?  Fr^GUENCia  AT  1M104  ft*  FRY 
HAS  CALCULATO  the  SPECTRAL  COflPCNEHTS.  ^  ^ 

PREO  •  T«  FREOUEHCY  OF  T>«  TEST  SISNAL  IN  HZ. 

AW  •  THE  ANPUITUDC  OP  THE  TEST  SIGNAL  IN  D8N*. 

PS  •  THE  SANPU  RATE  IN  IPS. 

R  •  the  pear  UALUE  of  THE  TEST  SIGNAL. 

^‘eScO^  IWtGRATOR 

TC  •  THE  COHPANBINO  SPEEB  OP  THE  SYLIAIIC  ItfTEGUATOR  IN  SEC. 

^  •*  ^SEOL'ENC"  UARIASLE  USED  TO  tETSPNlFt  THE  TEST  STOHV  SPEC- 
SsPoSwTS  *1-S0  TO  CETEWlire  THE  K»«NONIC  SPECTRAL 

**  FREOUC^.^^*^  ^  ^  output  SPECTRAL  CONPONENT  AT  T*€  TEST 

■EP  •  THE  OUTPUT  SPECTRAL  CONPONENT  POUER  IN  »R. 

«fl  •  THE  RUNNINO  SUf  OP  T>«  POWER  OP  THE  HARNONIC  COWOCNTS. 

THD  •  THE  TOTAL  HARNONIC  OlfTCRTICN  IN  N. 

•CTR  •  NORTMLIZEB  3  01  FRCOUENCY  OP  TIC  OUTPUT  FILTER 

•RlflR  ■  T>C  NORHALIZEP  ROLL-OPP  lAHBWieTH  OP  T>C  OUTPUT  FILTER 

RATIO  •  T>C  RATIO  OF  THE  N.SXINI/1  ETr*  SITF  TO  THE  NINIIUI  RTTP 
513  IN  THE  CVSO  CNCOOCR  ANS  C£»..ER;  QIvIn  W  "l. 


iittttittitttytt  susroutines  uso  nsnnsxrasssnn 

CNCOOEl  •  T»C  CWSD  ENCODER 
KCOKl  >  THE  CUSO  DECODER 

'*“2^Jp^®«**TCS  WHAN  AW  WHIN  USES  IN  THE  CWSD  ENCODER  M« 

PLTERGCN  •  TIC  COEPPICEKT  QEICRATOR  POR  TIC  OUTPUT  FILTER 
FILTER  •  FILTERS  THE  OUTPUT  SIGNAL. 

POICIER  TMNSFORN  SUMOUTIIC  PROP  THE  INBL 


C<ftt«(ttitytrt«ttti.«;'^trrTtrr<-:;rTX7«'vTYnrTTfTTrTT7?tTTttmxsrTZtim 
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DO  IS  J  •  l.KF 

IF  (  F  .NE.  FRE3UC(J>)  «0  TO  IS 
SUH  •  sun  ♦  fSX(J) 

IS  CONTINUE 
M  CONTINUE 

0 -  CAICUIATE  TOTAL  HARNONIC  DltTORTION 

THD  •  SORT  (SIN)  /  n  «  IM. 

URlTECS.COe)  TKD 

•M  FORflAT(lX,*T»C  TOTAL  HARnONIC  DISTORTION  IS  *.Ft.l,*K.*> 
DB 


I 


t 

I 


r 


OU  UUO  UU  OOtf  W  OO  V  O  O  O  <>  ’i  o  o  u  oo  o  o  oo«>  oo 


APPENDIX  rt 


Total  >!ar''or!ic  ''irtortion  v?; .  Tn-v't  ^o  rT^l  Povr^r 

PWCRAH  tTHD( INPUT, 0UTPUT,TA«5«lNPUT.TAPCS»0trrPUT,Pt0T) 

C - TWO  OS.  INPUT  POUER - 

C  THIS  PROGRAN  INUESTICATES  THE  UARlftTICN  IN  HAANONIC  DISTORTION 
C  AS  THE  SIGNAL  INi^UT  PC.-*?  IS  U.iRlCD.  THE  In.^uT  PC'JER  IS  C-J-NCCD 

C  IN  .4  03  STEFS  FKC.I  D’Ta  TO  i  0T*a.  TH-  HARHONIC  DISTORTION 

C  IS  THE  OUTPUT  IS  THEN  PEASUtB  AND  PLOTTED. 

C  THE  PROGRAN  IS  REPEATED  THREE  TINES.  STEPPING  THE  STEP  SIZE 
C  RATIO  FRCn  33  03  TO  33  C8.  THE  THREE  SETS  OF  DATA  ARE  THEN 

C  PLOTTED  ON  THE  SA,£  GRAPH. 

cnt»s»sssx<»xs»«t  VARIABLES  ttxtttntttnstntstxtttxtxxxtxxttsxxs 


INPUT  -  AN  INTEGER  ARRAV  CONTAINING  THE  BINARY  OUTPUT  OF  THE  CUSD 
ENCODER. 

OUTPUT  •  A  REAL  ARRAY  CONTAINING  THE  TIPE  FLHCTICN  OUTPUT  OF  THE 
TEST  SIGNAL  CENiSATCS  AND  AFTER  PROCESSING,  THE  TINE  FUNCTION 
OUTPUT  OF  THE  DECODER. 

POUER  ■  A  REAL  ARRAY  CONTAINING  THE  POUER  THAT  EACH  SANPLE  IS 
TAEEN. 

PSX  •  A  REAL  ARRAY  CONTAINING  THE  OUTPUT  SPECTRAL  POLtR  COrtPONEHTS 
OF  THE  CECCCER  CJTPUT  AFTER  PROCESSING  BY  THE  FAST  FOURIER 
TRANSFORN  SUBROUTINE. 

IWC,  UK,  CUK  •  WORKING  ARRAYS  USED  BY  THE  FFT  SUBROUTINE. 

FREOUC  •  A  REAL  ARRAY  CONTAINING  THE  FREQUENCIES  AT  UHICH  THE  FFT 
»M$  CALCULATED  THE  SPECTRAL  CO-TPCNENTS. 

H  •  THE  NURBER  OF  TINE  SARPLES  TO  BE  TAKEN. 

mCO  «  THE  FREQUENCY  OF  THE  TEST  SIGNAL  IN  HZ. 

ARP  •  THE  ARPLITUOE  OF  THE  TEST  SIGNAL  IN  DBR«. 

FS  ■  THE  SARPLE  RATE  IN  BPS. 

A  •  TtC  PEAK  V»ALUE  OF  THE  TEST  SIGNAL. 

lUL.  UK.  CUK  •  WORKING  ARRAYS  USED  BY  THE  FFT  SUBROUTINE. 

FRE'Hie  ■  A  REAL  ARRAY  CONTAINING  THE  FREOUENCIES  AT  UHICH  THE  FFT 
HAS  CALCULATED  THE  SPECTRAL  CORPONENTS. 

N  •  THE  NURBER  OF  TIRE  SAHPLES  TO  BE  TAKEN. 

FREQ  >  THE  FREQUENCY  OF  THE  TEST  SIGNAL  IN  HZ. 

ARP  •  iHk  ARPLITUDE  OF  THE  TEST  SIGNAL  IN  DBRO. 

FS  •  THE  SARPLE  RATE  IN  BPS. 

A  •  T>C  PEAK  UALUE  OF  THE  TEST  SIGNAL. 

FCl.  FCa,  FC1  •  ROLL-OfF  FREOUENCIES  FOR  THE  PRINCIPLE  INTEGRATOR 
IN  THE  CUSD  ENCODER  AND  DECODER. 

TC  •  THE  CORPARDING  SPEED  OF  THE  SYLWBIC  INTEGRATOR  IN  SEC. 

I,  J.  K  ■  COUNTING  INDICES  FOR  THE  VARIOUS  ’DO*  LOOPS. 

F  •  A  FRECUENCY  VARIABLE  USED  TO  DETERRIHE  THE  TEST  SIGNAL  SPEC¬ 
TRAL  COP.PONENT  and  also  TO  DETERRINE  THE  HARRONIC  SPECTRAL 
CORPOHEHTS. 

CS  ■  THE  RRS  POWER  OF  THE  OUTPUT  SPECTRAL  CORPONENT  AT  THE  TEST 
FREQUENCY. 

C  REF  •  THE  OUTPUT  SPECTRAL  CORPONENT  POUER  IN  DBR. 
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C  SW1  •  THE  RUHHIHO  SUH  OF  THE  POUER  OF  THE  HWJflOHlC  COMPONENTS. 

TKO  •  AN  ARPAV  CONTAININQ  THE  IMLU£  OF  HARMONIC  DISTORTION  AT  EACH 
LEVEL  OF  INPUT  POUtR. 

•  •  AN  ARRAV  CONTAINING  THE  OUTPUT  FILTER  COEFFICIENTS. 

HP  >  THE  NUMBER  OF  OUTPUT  FILTER  COEFFICIENTS. 

BETA  •  THE  NORNALIEED  CENTER  OF  THE  TRANSITION  BANS  FOR  THE  OUTPUT 
FILTER. 

C  CAHMA  •  THE  NORMALIZED  UDTH  OF  THE  OUTPUT  FILTER  TRANSITION  BAND. 

C*ttXtXtXtXttXtXtttXtTttiXttXttttt1ttXJtttXXStttXttttXXXTtXXS»tXXXSXmS 

cxxxxxxtxxxxttxxxxxxxxxxxxxxsxxxxxxxxxxxxxxxxixxxxxxtxxzxxxxxxxtxtxxxxxx 


c- 


■PROGRAM  START' 


C-~—  INITIALIZE  VARIABLES  M4D  ARRAVS 


DIMENSION  INPUT< ) ,  0UTPUT<5«««  > , POtiOfCBeZ ). PSN(  ISO ) 
l4lUK(E0>,U<(r>a),Ff.£Ct'£(15«>,  THDiEOB).  B(E20) 


COMPLEX  CUK«3Ja) 
MtOSnOJ  •  SCAT (10. 


Xt  ((DBMO  -  4.)  /  10. )  t  .001  *  BOO.)  t  SORT 


1 1 

KHAR  •  -1 


INPUT  UORKING  VARIABLES 

MiNT*x, ’^Dynamic  range  test  at  *,fs,’  bps  and  '.freoi,*  hz* 

READ  >.  FCl.  TO 

READ  *,  XLEN,  VLEN.  XMIH,  XMAX,  VMIH,  VMAX 

xstep  •  (xr.AX  -  xniN)  /  xlf.n 

VSTEP  •  (VMAX  -  VNIN)  /  VLEN 
PRINT  f,*  TO  •  •,TC,’,  FCl  •  SFCl 
READ  B.BETA,  GArVIA 

PRINT  t,*  CETA  •  ',SSTA,',  GAMMA  •  SOAMMA 
CALL  FLTROEN(B£TA,GAfWA,NP,B) 


START  LOOP 

DO  1000  NR  ■  2. 6, a 
RATIO  *30.  ♦NR 

CALL  VMAXOPTC VMAX, VMIN.FS, FCl, TC.RATIO) 

DO  500  IS  •  1,109 
POUER(IS)  •  -49.  ♦  .4  *  IS 

determine  peak  value  OF  TEST  SIGNAL 
AMPl  •  A(POUER(IS)) 

GENERATE  INPUT  TIME  FUNCTION 

CALL  SICr«L(0UTPUT,50«9,FS,FRE01,0.,AMPl,0. » 

PROCESS  THE  TIME  FUNCTION  THROUGH  THE  CVSD  SYSTEM 

CALL  EHCOCEl (OUTPUT,  1*'PUT,SC0»,FS, FCl, FC3,FC3,TC,UNAX,VMIH, DC) 

MLL  DECOCEl  ( INPUT,  OU  ’PUT,  5000,  ,  FC2, FC3,  7C,  IWAX,  V« IN, DC  ) 


C -  FILTER  ne  OUTPUT 

CAU  FIlTER(0UTPUT,5O90,NP,f» 

C -  DETERMINE  THE  SPECTRAL  COMPONENTS  OF  THE  OUTPUT 

CALL  FTrpS(0UTPUT,DUM,409S,2S0,0,PSX,DUM,DUM,l«C,WC,C«C,IER) 

C— -  DETERMINE  THE  COMPONENT  AT  THE  TEST  FREGUCNCY  AND 
C  CALCULATE  REFERENCE  VALUES. 
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»;f  •  0 

00  3  K  •  3.189,2 
KF  •  KF  1 
F  •  »(IC-1  j/ass. )  <  FS 
FREOUE(KF)  .  f 

IF  (P5X(K>  .LE.  fi.E-19)  PSX(K)  ■  6.E-13 
PSX(ICF)  ■  PSXCK) 

IF  (F  ,ME.  FRtQl)  GO  TO  8 
Ee  •  SQRT(  PSX(K)  ) 
t  CONTINUE 

C - CAUULATE  POUER  AT  HARTIONIC  FREOUENCIES 

sun  -  8. 

DO  ae  1  •  2,10 
F  •  FREOl  /  I 

IF  ((F  .LE.  leO.)  .OR.  (F  .GT.  4000.))  CO  TO  20 
DO  IS  J  •  l.KF 

IF  (F  .NE.  FREOU£(J))  CO  TO  IS 
sun  •  SUN  ♦  PSX(J) 

IS  CONTINUE 

20  CONTINUE 

DO  30  I  •  2,30 
F  •  FREOl  t  I 

IF  ((F  .LE.  ICO.)  .OR.  (F  .GE.  4000. >>  GO  TO  30 
DO  25  J  •  I.ICF 

IF  <  F  .NE.  FREQUEfJ))  GO  TO  2S 

SUN  •  sun  ♦  psx(j) 

2S  CONTINUE 

30  CONTINUE 

C -  CALCULATE  TOTAL  HARflONlC  DISTORTION 

THD<IS)  •  SORT  (SUN)  /  EO  *  100. 

IF  (THDdS)  .GT.  100.)  THD(IS)  •  100. 

500  CONTINUE 

C -  PLOT  RESULTS 

ICHAR  *  JCMAR  4  1 

CALL  PLTRANG  ( POUER,  THO,  100,  ICHAR, XHIN.XLEN,  XSTEP,  VTIIN,  VLEN,  VSTEP ) 
1000  CONTINUE 

CALL  PLOTE(N) 

END 


5U1R0UTIHE  PLTRAN0(X,Y,N,ICHA«,XNIN,XL£H, XSTEP, STUN, VLEN, VSTEP) 
C - DVNANIC  RANGE  PLOT  SUSROUTI* - 


C  THIS  SUBROUTINE  CREATES  A  PLOT  OF  THE  Y  ARRAY  VERSUS  THE  X  ARRAY. 
iuumutttttttttwxxs  VARIABLES  t(sssxutn»xxnsn»s*usn»xnM 
C  4  •  AN  ARRAY  CONTAINING  THE  ORDINATE  VALUES 

C  V  •  AN  ARRAY  CONTAINING  THE  VALUES  TO  BE  PLOTTED. 

C  N  •  THE  NUNBER  OF  VALUES  IN  THE  ARRAYS 

Ctxxxxxxxxxsxxxtxxxtxxttxsxxtxtxtxxtxxtxxxxxtxxxtxttnxttxxxxxixxttxxxxt 

CttXtXXXXXtXXXXXXXXSXXXXXXXXXXXXXXXXXXXXXXXXXXXXXiXXXXXXttXXXilXXXXXXXXX 

c - SUBROUTINE  START - 

C -  INITIALIZE  ARRAYS 

DIRENSION  X(4e90),  Y(4096) 

X(N*1)  •  XNIN 

X(N>2)  •  XSTEP 

V(N*l)  •  VniN 

Y(N42)  •  YSTEP 

IF  (ICHAR  .CT.  0)  GO  TO  500 

C -  ESTABLISH  NEU  PAGE  ORIGIN 

CALL  FACT0R(.5) 

CALL  P10T(2.,  2.,  -3) 
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C -  WX  IN  THE  GRAPH 

CALL  RECT(e..  0..  VLCN.  XLEN.  3) 

C -  DPAU  THE  AXES 

CALL  AXIS<0.0,0.«.18HINPUT  PO'x-ER  (I!B«0).-lB,XLE»<.e  a.XNIN.XSTEP) 
CALL  AXIS(Q.e,e.e,l4H0IST0RT10N  (ll),t'<<yL£N,»d.e.yflIN,VST£P) 

C -  PLOT  UALUES 

6M  CONTINL'E 

CALL  LZN£(X.y*N.l.l«.ICHAA) 

RETURN 

END 


12;? 


APPSriDIX  0 


"  1  r: t;r  li'^  d  d'otM]  Pi  r. ‘orl: ;  on  vo  .  Inn^;'.  Pj  'nr;  1  Po  •.  ;o  r 


PROCRAfl  nnTHD(INPUT,0(jrPirr,T<tf»ES-INPUT.TW>M»0OTPUT,PL0T) 

C - niSflATCHCO  THD  US.  INPUT  POWER - 

C  THIS  PROCRPP  INC'PSTIGATES  THE  IW9IATI0N  IN  HAIfflONIC  DISTCRTICN 
C  PS  THE  SICN^tL  IdrUT  PC'.'ER  IS  UPSIEE.  THE  Ii..’UT  =0..SR  IS  CHANGED 

C  IN  .4  DB  STEPS  FRC.I  -4)  DS"3  TO  e  E:r3.  THE  W^JNONIC  DISTCRTION 

C  IS  THE  OUTPUT  IS  THEN  rEASURED  AND  PLOTTED. 

C  THE  PROGRAN  IS  REPEATED  THREE  TlfES.  WHILE  THE  ENCODER  PARAPETERS 
C  ARE  HELD  CC.'iSTANT,  THE  DECODER  STEP  SIDE  RATIO  IS  ALLCOED  TO  UATV 
C  FRON  3a  D3  TO  E3  D3.  THE  THREE  SETS  OF  DATA  ARE  THEN  PLOTTED  ON 
C  THE  SANE  GRAPH. 

ctxmttt«t»»t»x»t«  UARIABLES  tttxnstxxxtxttxtttsxxttttttwtmxxs 

C  1»#»UT  •  AN  INTEGER  ARRAY  CONTAINING  THE  BINARY  OUTPUT  OF  THE  CUSD 
C  ENCODER. 

C  OUTPUT  •  A  REAL  ARRAY  CONTAINING  THE  TINE  F'JNCTICN  OUTPUT  OF  THE 
C  TEST  SIGNAL  GEtERATCR  AND  AFTER  PROCESSING,  THE  TIPE  FUNCTION 

C  OUTPUT  OF  THE  DECODER. 

C  POWER  •  A  REAL  ARRAY  CONTAINING  THE  POUER  THAT  EACH  SAFIPLE  IS 
C  TAICEN. 

C  PSX  •  A  REAL  ARRAY  CONTAINING  TT'E  CUT?!JT  SPECTIAL  Pd.ER  COrRONEHTS 
C  OF  THE  DECODER  OUTPUT  AFTER  PROCESSING  BY  THE  FAST  FOtSIER 

C  TRANSFORN  SUBROUTINE. 

C  IWC,  UK,  CUK  •  WORKING  ARRAYS  USED  BY  THE  FFT  SUBROUTINE. 

C  FREQUE  •  A  REAL  ARRAY  CONTAINING  THE  F’^EGUEHCIES  AT  WHICH  THE  FFT 

C  HAS  CALCULATED  THE  SPECTRAL  CCfiPONENTS. 

C  H  •  THE  NUriBER  OF  TINE  SAITLES  TO  BE  TAKEN. 

C  FREG  •  THE  FREQUENCY  OF  THE  TEST  SIGNAL  IN  HZ. 

C  AFIP  •  THE  AFIPLITUOC  OF  T>«  TEST  SIGNAL  IN  DBNG. 

C  FS  •  THE  SARPLC  RATE  IN  BPS. 

C  A  •  THE  PEAK  UALUE  OF  THE  TEST  SIGNAL. 

C  FCl,  FC2,  FC3  •  ROLL-OFF  FRtOUENCIES  FOR  THE  PRINCIPLE  INTEGRATOR 
C  IN  THE  CUSD  ENCODER  AND  DECODER. 

C  TC  •  THE  COUPANDING  speed  of  THE  SYLLABIC  INTEGRATOR  IN  SEC. 

C  I,  J,  K  -  COUNTING  INDICES  FCR  THE  UAfilOUS  ’DO'  LOOPS. 

C  F  •  A  FREQUENCY  UARIABLE  USED  TO  CETEFNINE  THE  TEST  SIGNAL  SPEC- 

C  TRAL  COr.PONENT  AND  ALSO  TO  DETERNINE  THE  HARNONIC  SPECTRAL 

C  CONPONENTS. 

C  EB  •  THE  RNS  POUER  OF  THE  OUTPUT  SPECTRAL  CONPOfCNT  AT  THE  TEST 
C  FREQUENCY. 

C  REF  •  THE  OUTPUT  SPECTRAL  CCflPONENT  POWER  IN  DBN. 

C  SUH  •  THE  RUNNING  SUN  OF  THE  POUER  OF  THE  HARMONIC  COMPONENTS. 

C  THD  •  AN  ARRAY  CONTAINING  THE  UALUE  OF  HARNONIC  DISTORTION  AT  EACH 
C  LEVEL  OF  INPUT  POWER. 


CtXtXXXXXXXXXXXXXliXXXXXXXXXXtXXXXXZtXXXXXXXXXXXXXXXXtXXXXXXXXXXXXXXXXXX 

CXXtXSXXXXXXXXXXXXlXXXXXXXXXXXXXXXiXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 
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■PR05RW1  STORT- 


C' 

C - IM1TIRLI2E  MID  ARRAYS 

DiriEHSICN  INPUT* BC'’9). 0UTPUT(ra»9  ).POt;HB(88a).f>SX(  150  » 
l,ll!S(t0),W<<lSt!).P«t(ilE{150).  THDREFCcea),  TVID(20a),  BJ2W) 
COnPLEX  CL'XC  J?) 

AOBW)  •  SCAT*  10.  Xt  ((OBHB  -  4.  )  /  10. )  Z  .001  t  600.)  Z  S0fiT(2. 
1) 

ICHAR  •  -1 


C -  INPUT  yORtCINO  UARIASLES 

READ  Z.  FREOl,  FS 

PRINT  Z.*  DYN.V1IC  RPNGE  TEST  AT  '.FS,*  BPS  AND  ’.FREai,*  HZ’ 
READ  Z,  FCl.  TC 

READ  Z.  XLEN.  VLEH.  XNIN.  XPAX,  yfllH,  VNAX 

XSTEP  •  (X'^X  -  X.-'IN)  /  X'.EN 

VSTEP  •  (VllAX  -  VNIN)  /  VLFH 

PPIHT  I,*  TC  •  *,TC,*,  FCl  •  ’.FCl 

PEhD  I.BETa,  GAWIA 

PRINT  »,•  BETA  •  '..ETA,*,  GANNA  •  GANNA 
ChLL  FLTRCENIBETA.CA.VW.NP.B) 


C— -  STmRT  loop 

DO  1000  NR  •  a. 6. 2 
RATIO  •  33.  +  NR 

CALL  WVAXOPT(UnAX.ONIN,FS,FCi,TC,RATIO> 
IF  <ICHAR  .GE.  0)  GO  TO  2 
EVNX  •  UNAX 
EUNN  •  UPilN 
2  CONTINUE 

DO  500  IS  •  1.100 

POUERdS)  •  -40.  .4  Z  IS 

C -  DETERNINE  PEAK  UALUE  OF  TEST  SIGNAL 


ABPi  •  AlPOUtRaSn 

c -  GENERATE  INPUT  Tlf€  FUNCTION 

CALL  SlGNAL(OLrrPUT,5Oe0,FS,FREQl,0..ANPl,0. ) 

C -  PROCESS  THE  TINE  FUNCTION  THROUGH  THE  CVSD  SYSTEM 


CALL  ENCODPl < OUTPUT , INPUT , 5030, FS, FCl ,FCa, FC3, TC, EUNX, EUNN, DC ) 
CALL  DECODEH INPUT, OUTPUT, 5ee0.FS,FCl,FC2,FC3,TC,VrtAX,UNIN,DC) 


C -  FILTER  the  OUTPUT 

CALL  FILTER(OUTPUT,5C00,NP,B) 

C -  DETERNINE  THE  SPECTRAL  CCrPONENTS  OF  1>C  OUTPUT 

CALL  FTFPS(OUTPUT,OUN,4Og6,25S,0.PSX,DUN,DUN,IUK,UK,CUK,IER) 

C -  ELININATE  THE  CO^PO’ENTS  AT  ODD  NULTIPLES  OF  THE  SAMPLE  RATE. 

C  DETERNINE  THE  CO'  ^OT'ENT  AT  THE  TEST  FREQUENCY  AND 

C  CALCULATE  REFERENCE  VALUES. 

KF  •  0 

DO  8  K  ■  3,128,2 
KF  •  KF  ♦  J 
F  •  ((K-1 1/256. )  Z  FS 
FREOUE(KF)  •  F  ^ 

IF  (PSX<K)  .LE.  6.E-10)  PSX(K)  •  S.E-10 
PSXIKF)  •  PSX(K) 

IF  (F  .NE.  FPEOl)  GO  TO  8 
E0  •  SORT!  PSX(K)  ) 

■  CONTINUE 

C -  CALCULATE  POIdCll  AT  HARMONIC  FREQUENCIES 
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IS 

2» 


2S 


sun  •  e. 

w  20  I  ■  a,i0 
F  •  FREOl  /  I 

DO  Jr®’’  <»  TO  28 

InMz^uu  CO  TO  IS 

^htIhT  * 

CONTIN'JE 
DO  33  I  .  a. 30 
F  •  FREOl  S  I 

M  CO  TO  38 

IF  <  F  .NE.'fkEOUE<J))  CO  to  2S 

sun  •  sun  ♦  psx(J) 

CONTINUE 

continue 


c— —  calculate  TOTAL  HARWONIC  DISTORTION 


TMD(IS)  •  SORT  (Sun)  /  E0  t  100. 

^  IF  (THO(IS)  .CT.  100.)  THD(IS)  •  100. 
SOe  CONTINUE 


C— -  PLOT  RESULTS 


ICKAR  •  ICHAR  >  1 

S5J;V.Ef-P*''^‘^<^'^^‘*'‘'^0«I®®*ICHAR.xniN.>aEN,XST£P,vniH,VL£N,VSTEP) 

lVv9  COnTInut 

CALL  PLOTEIN) 

END 


CttxtxtttzttittttxittxtttttxxxtxtwxtTtttntitttxtxttttttxtxtTizrixx  m 
CttXXXtttXXXXXXtXXXXXXitXXZXZSXXXiXXXlZiXXXXXXXZXX'iliXXXXZXlXSXZZiXXT  izx 

SUBROUTINE  PLTRANC(X,Y.N.ICHAR.XniN,XL£N.XSTEP,vniN,YLEN,YSTtP) 

c - DYNAniC  RANGE  PLOT  SUBROUTINE - 

C  THIS  SUBROUTINE  CREATES  A  SEnl-LOC  PLOT  OF  THE  X  AND  V  ARRAYS. 
cttxttttsxxtttxttxzxxtxxt  UARiABLES  xnsnms»msmzn>su»mzBm 
C  X  •  AN  ARRAY  CONTAINING  THE  ORDINATE  UALUES 

C  V  •  AN  ARRAY  CONTAINING  THE  VALUES  TO  BE  PLOTTED. 

C  N  •  THE  NUflBER  OF  VALLES  IN  THE  ARRAYS 


Cmsi;snt»»tztzxzxttrrz»xzxxxzs:sz*xxz:z»zx»»tnznnzTnztzzsttt 

CtXSXXXtStXXSXZXXXXXXZXXXXXXXXZJXXXZXllZXXXXXXXXXXXiXXXlXXiXXXXXZXXXZtXt 

C - SUBROUTINE  START - - - - 

C -  INITIALIZE  ARRAYS 

DIHENSION  XC^ese).  y(4096) 

xiH+i)  .  xniH 

X(H+a)  •  XJTEP 

V<N+l)  •  VNIH 

V(H»a)  •  YSTEP 

IF  (ICHAR  .GT.  0)  CO  TO  S88 

C -  ESTABLISH  HEU  PAGE  ORIGIN 

CALL  FACTOR (.5) 

CALL  PL0T(2..  2..  >3) 

C -  BOX  IN  THE  GRAPH 

CALL  RECT(0..  0.,  YLEN.  XLEN,  9.,  3) 

C -  DRAW  THE  AXES 

CALL  AXIS(0.0,0.0,18HINPUT  POL’ER  <DBn0).-18,XLEN,0.0,XI1IN,XSTEP) 
CALL  AXIS(e.0,0.8.1'4HOISTOR7ICN  (»),10.YUN,6«.8.yniN.ysTEP) 
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PLOT  UPLUES 


CALL  L1Nc(X,V.H,1,1«.ICHAR) 

RETURN 

END 


APPKMDIX  P 


In torriodulation  ■’'i r;ljOri; ion  i’ni  ; 


lKTEimO<lMPUT.OUTFUT,T*«S«I»^T,T««*.etfmn> 

C . --INTERnOOULATlOM  DI8T0RTI0M  - 

C  this  program  CALCUU»TeS  Tvc  iKTPC-OOULATtOH  DISTOOTTM  OF  A 
c  cwsD  svsTtfl  I  '  r  r-'u  e*-cc“:r  ►  o  :  A-t  ccr.-ECTEO  »ack-to- 

C  »ACK.  OISTCr^ITH  IS  E3  SV  I'.’UTTUG  A  TT.':’’  CE  -PCSO) 

C  OF  TUO  EGJA'.  t  -LITL  '  S;-'  .D!:3  AT  i;,.l  v-z  pfj  7"1  K.;.  Th  ■ 

C  AJIPIITL'EE  Of  Z~z  tIEF-.r-*'--  r.  --XT  is  Zy  S  t  aS'.  -;0  f.-fo  COIVAJSB 

C  TO  TXE  IhPUT  Sli.SAL  TO  ti7d>.SlF;£  TViC  FEkCiKT  DliTG-iTIOfl. 

CttXUMttttcnsXMSU  UARIAIl£S  mxmnntKttnmwtttttunmtnn 


C  imn  •  AN  IKTC^  MUMV  CONTAININQ  T»C  tlNAAV  OUTPUT  OP  THE  CU6I) 
C  CNCOOCR. 

C  OUTPUT  •  A  F.'^AL  AT’AV  CONTAINING  THE  TIFE  PJNCTICN  OUTPUT  OF  Tv€ 

C  TEST  SlCnAL  ClfETflTCT  ASO  AFTER  FflOCESSINQ.  TFC  TINE  FUHCTICM 

C  OUTPUT  OF  THE  CECCikJJ. 

C  PSX  •  A  REA'..  ATPAV  CCNTaININO  TwE  CLfTTUT  SPECTRAL  POtGJ  CONPOfCMTS 
C  OF  TFJE  rZCCiER  OUT'-.rr  FFTER  PFXCESSIW  by  the  FAjT  fourier 

C  TRANSFOWI  (FFT)  SLZXUriME. 

0  tUK,  UK.  CUK  •  UORKINQ  ARRAYS  USES  IV  THE  FFT  SUBROUTINE. 

C  FREOUE  •  A  CEAL  ARRAY  CONTAINIHO  Tvc  FPE'-JEHCIES  AT  UHICM  THE  FFT 

C  HAS  CALCULATES  Th£  SPECTRAL  CC.'.PCNEHTS. 

C  N  •  T>C  NUnSER  OF  TITZ  SAMPLES  TO  IE  TAKEN. 

C  AHPl  •  TFC  AflPLITUDCS  Of  TFC  TEST  SIGNALS  IN  DIM. 

C  FS  •  THE  SAMPLE  RATE  IN  IPS. 

C  PEAK  •  nC  PEAK  VALUES  OF  THE  TEST  SIGNAL  COMPOtCNTS. 

C  A  'D*  OR  'E*  PFSFIX  ON  THE  NEXT  THPEE  SETS  OF  VARIABLE  IHSICATES 
C  THE  UARIAILE  IS  U3E0  lY  EITICR  TK£  UCCSER  OR  EHCC:£R.R£SPEO> 

C  TIVELV 

C  FCl,  FCa.  FC3  •  P-IL-CFF  FPEOUENCIES  FOR  THE  PRINCIPLE  INTEGRATOR 
C  IN  THE  CVSD  ENOOOER  AND  CECODCR. 

C  TO  •  THE  COMPANDINO  6PCEB  OF  T»C  SYLLABIC  INTEGRATOR  IN  SEC. 

C  RATIO  •  TFC  FMXIMUN  STEP  SIZE  TO  MINIMUM  STEP  SIZE  RAHO  IN  Cl 

C  F  •  A  FREOUENCY  UARIA.3LE  USED  TO  rSTERMIHE  THE  TEST  SIGNAL  SPIC- 

0  TRAL  C\  .-CrCNT  AND  ALSO  TO  CETEkMINE  THE  FMRMONIC  SPECTRAL 

C  COMPONENTS. 

C  REF  ■  TFC  OUTPUT  SPECTRAL  COMPONENT  POUER  IN  DIR. 


:  iUM  .  THE  SUM  OF  THE  POUER  IN  THE  TEST  SIGNAL  COMPONENTS. 

•  NP  •  TFC  NUMSER  OF  FILTER  COEFFICIENTS 

C  •  •  AN  ARRAY  CONTAINING  TFC  INPUT  AND  OUTPUT  FILTER  COEFFICIENTS 
•ttsixttttxttxittus  SUBROUTINES  USED  ttxaxxsxxxxxaaxxxxuuxxtxxx 
C  ENCODCl  •  THE  CUSO  ENCODER 

C  DECODE 1  •  THE  CVSD  DECODER 

C  W1AXOPT  •  DETERMIFCS  THE  VALUES  OF  VMAX  AND  VMIN  USED  IN  TFC 
C  ENCODER  AND  DECODER 

C  FLTRGEM  •  GENERATOR  OF  TFC  OUTPUT  FILTER  COEFFICIENTS 

C  FILTER  •  FILTERS  TFC  OUTPUT  USING  TFC  COEFFICIENTS  PROCOCED  SV 
C  FLTRGEN 

C  SIGNAL  ■  TFC  TEST  SIGNAL  OENERATOR 

C  FTFPS  •  THE  FAST  FOURIER  TRANSFORM  FRCP  THE  IMSL  LIBRARY 

ct»«tssms<xxttuftx»*nzxnuunuxusnmmnnt>xsx8zaxtuMB 
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c -  'Hwaww  rr««T - - - 

C -  INITIALIZE  VMIAILES  AND  ARRAYS 

DINENSION  INPUT(50^«>.0LT1>UT(Se«*>.PSXtl5S) 
i,iuK(ea),uK(is«).Fi^ouc(2M).  »(a«e) 

REAL  ir«D 
COMPLEX  CUX(30«) 

ACDSnS)  •  SORT  (iS.  t<  ((DIN*  -  4.  VIS.)  t  .Ml  S  CM.)  t  SORTIS.  > 

C -  iwvr  AND  PRINT  UORKINQ  UARIADLES 

READ  t.  ANPi.  rs 
READ  »,  EFCl.  ETC,  rUTIO 
READ  ».  tFCl,  DTC,  C-iATlO 
READ  t,  seta.  r-."<A 

PRINT  t,'  lhT£-^0  -rST  FOR  INPUT  SIONALS  Of  7F#  <WD  ISM  HZ  AT  *, 
lANPi,*  CBfre,  A.M)  S.V."LE  KATE  .  •,FS.*  tivS* 

PRINT  ».•  ETC  •  •.rrc,*,  E.-Cl  .  EPATIO  •  '.DJATIO 

PRINT  t,*  DTC  •  •.DTC.*,  tFCl  •  •.CPCl.*,  C-iATlO  •  '.DRATIO 
PRINT  t,<  SETA  ■  ".BETA.*,  CFVtflA  •  ‘.CANMA 

e—  DCTERNINE  PEAK  UALUE  Of  TEST  StCNAL 

PEAR  •  A  (ANP1> 

©— —  OENERATE  INPUT  TINE  FUNCTION 

CALL  SIONAL(OUTPUr.  5«M,  fS,  75*.,  ISM..  PEAK,  PEAK) 

C -  OEICRATE  FILTER  COEFTICIENTS  AND  CUSD  SVSTEN  PAAANCTERS 

CALL  OnA>»PT(FW*X.nw,FS,Frci,rrC,EPATIO> 

CALL  «1AXOPT(C.-X,':v  'N.F9,:fCl,OTC,C3ATIO> 

«.aLL  FLTRCENI  beta, GANNA, NP, I) 

—  PROCESS  the  tine  FUNCTION  T)«OUGH  THE  CVSO  SYSTIN 

Call  ENCOOEl (OUTPUT. INPUT, WS.FS. EFCl, FC3,FC3, etc, EWK.ELNN, DC) 
CAU  DECODE! <  INPUT, 0UTI>UT,6iM,fS.0FCl,FCS,FC3,0TC,DUR<.DtW,DC) 

c - FILTER  the  output 

CAU  riLTeR(ourpur.sM«.NP,a) 

C— —  SUITRaCT  the  AUERAQE  UALUE  PROA  THE  SANPU  STRINC 

SUN  •  S. 

DO  «  I  •  1,4696 
SIR)  •  SUN  ♦  OUTPUTCI) 

C  CONTINUE 

AUER  •  sun  /  4996 
DO  7  I  •  1,46CS 
OUTPUT(I)  •  OUTPUTd)  -  AUER 
7  CONTINUE 

e -  CALCULATE  THE  SPECTRAL  CONPONEHTS  OF  THE  OUTPUT 

CAU  FTFPS(OUTPUT,  Dun.4696.RSS.  S.PSX,  OUH,  DUN,  lUK .UC.CUC,  lER) 

C  DETERniNE  the  COnPOTCNT  AT  THE  TEST  FREQUENCY  MO  CALCULATE 

C  THE  RNS  UOLTACC. 

sun  •  6. 

DO  I  K  •  2.139 
r  •  <(K-l)/cS6.)  X  FS 
fREOUE<K-l)  •  F 
PSX(K-l)  •  PSX(X) 

IF  ((F  .HE.  TJO.)  .AND.  <F  .HE.  ISM.))  00  TO  S 

sun  •  sun  ♦  SGRT  (  psx(K)  > 

REF  •  IS.  S  AL0Q16(PSX(K>  /  6M.  <'.M1) 

I  COffTlNUC 

C -  CALCULATE  INTERnODOLATlOn  DISTORTION 

BIf  «  2SS. 

DO  IS  I  •  1.12S 

If  (FREOUE(t)  .NC.  OIF)  M  TO  IS 
OIS  •  PSX(I) 

»  CONTITKJE 

inOD  •  SORT  (DIS)  /  Sin  t  IM. 

If  ItflOO  .or.  1(T«.)  INOO  •  ISO. 

UR1TE(6,6M)  ITWD 

MS  FORHATdX.'THE  INTERnOOUlATlON  OISTORTION  It  '.fS.I.**.') 

M 


1 
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Ai’PKnriix  Q 

In  I  0':P:  1 . 1 '  i  (.’1  "  ■  -  i  o;-  ‘  i  or'  Tn'  -  .1 


ntocMM  oiNoo(imT.oum;r.TMc«*otrrvur,nAr) 
c- - tHTtgnoa  w*.  imir  wun - 

C  this  WJOOWW  W*Sl!«i  IKTtWf^CUtATlOM  BISTOTTIOM  OS  •*  FDHCTICH  Of 
C  INPUT  sif'fPL  PC  I  IN  p  C-.  o  ■."7  i  T,';  r.:,^c::,s  pno  rrccrai 

C  PPRANC-rtlW  Pi<E  Pi^TECTLV  P^TO'O.  C  ’  CC;- -A  T I  i  ha  fc.,'  ^^-EpTEP  T>*«* 

C  TINES  AS  TME  STEP  SIZE  RATIO  la  M^lEJ  F«.JM  3J  DS  TO  .3  OS. 

C  nc  TMtCZ  SETS  OS  DATA  A«  THEM  PlOTTtO  ON  THE  SANK  OMAN. 

ctmutntxstttttxtttts  uaaiaius  tttn-ntTTTWtinimmmtmtms 

C  X  «  AN  ARRAY  CONTAININO  THE  TEST  SZ6NAL  INPUT  POiCII  AT  UHICH  HEA- 
e  SURENENTS  ARE  NACE. 

e  V  ■  AN  ARRAY  CONTAININO  THE  1NTE38IOOUIATION  OXSTOmON  HEASURS- 
e  rCNTS. 

e  ■  •  AN  ARRAY  CONTAININO  THE  OUTPUT  PILTER  COeFTlCIEMTS. 
e  INOO  •  INTERNOOULATION  DISTORTION  AT  TIC  PRESENT  TEST  SIONAL  PMER 
C  PS  •  THE  SANPU  RATE 

CO  THE  FOLLOgiMa  UARIAJUS  USE  AN  •£•  AND  A  'O*  P«^tX  TO  INDICATE 
C  use  tY  CITKER  the  ENCCCER  CR  CECOSCR.  SESPECTIUEUY. 

C  PCI  •  ROU-OPP  FREOUENCY  OP  THE  PRINftRY  IKTECRATOR 

C  TC  •  THE  TINE  CONSTANT  OP  THE  SYtLASIC  FILTER 

e  «MN  ■  THE  NAXINUn  INPUT  TO  T>X  SVILASIC  FILTER 

«  (RM  •  THE  NIHIHUN  imiT  TO  THE  SVLLARIC  FILTER 

e  KTA  •  THE  NCRNAL12C0  CENTER  FREQUENCY  OP  THE  OUTPUT  HLTER  TRANt- 
C  ITION  BAND 

e  OAnW  ■  1>C  NORNALIZEO  WIDTH  IP  THE  OUTPUT  PtLTER  TYANSITtON  BAND 

coBimuBtimnotw  susroutiiks  used  moommrionioioiim 

C  FLTROeN  •  T>e  OUTPUT  FILTER  CCCFPICIENT  OENERATOR 

»*  FILTER  •  the  SuraoUTIK  THAT  PILTtAS  THE  «C0Da  OUTPUT  SIOHM. 
USING  the  CC€Fr:CliNTS  CALCULATED  BY  fLT,<8IN. 

:  UNAXOPT  •  CAlC'JLATtS  UNAX  AND  UNIH  FOR  THE  CWSD  ENCODER 

C  AND  DCCOCER  SUBROUTINES. 

e  INTCRNO  •  THE  SirtROUTINE  THAT  CALCULATES  THE  INTERNODUUITION  DIS- 
C  TORTION  IN  T>£  CUSD  SYSTCN  OUTPUT  SIONAL. 

e  FACTOR,  PLOT.  SCALE,  AXIS,  RECT,  LINE  •  CALCCNP  PLOTTER  ROUTIICB 

C— - PROQMM  STMT - 

C-~-  INITIALIZE  VARAIBLU  AND  ARRAYS 

DINENSION  XdM),  Y(tM»,  KMO) 

KAL  INOO 
ICHAR  •  -1 

C - READ  AND  PRINT  TIC  HORKINQ  UMIRSLES 

READ  t.  FS 

READ  t,  CFCl,  ETC 

READ  t,  CFCl,  DTC 

READ  (,  BETA,  <«4NNA 

PRINT  t,*  FS  •  ■,FS 

PRINT  •,•  EFCt  •  ',CFCt,*,  ETC  •  •,ETC 

PRINT  t,*  CFCl  •  M'Cl,*,  CTC  •  VDTC 

PRINT  I,*  BETA  •  META,*.  OANNA  •  *,aANNA 

GENERATE  LOU-PASS  FILTER  COEFFICXENTf 

CALL  FLTR6eN(BETA,aAmA,M,B> 

e -  START  CALCULATION  LOOP 

90  MM  >  f.S,I 


IPO 


CMUuum  cvf»  tvmn 


ICHM  •  ICHM  *  t 
DMTIO  •  S*.  ♦  J 

CALL  v»wxcirT(cwn.cunM.rf.crci.rK;.OMrioi 

IM  cociTinuz 

MATIO  •  r4.  *  J 

CMX  MMxoprioum.oum.pt.Brct.o'‘e,MMTie> 

CMXUum  tNTcsntuutTioH  BirronrioH  ut.  tmrr  Mun 


M 


I  I  • 
-4*. 


I.IM 
*  .4  S 


1 


Kl  I  I  •  AfWt 

CALL  tNTcmo(r«,«#»i.crei.cTC,cum,cun«.Brci.B)Tc,nOTc,DMm,NP,B,in 
tool 

vitl  •  INOD 
COMTtIUI 


e —  nor  nc  nwLrt 

If  (tCHMI  .or.  •>  00  TO  BM 

MU  rit:TC4(.s) 

MU  ^toT<a.,a.,-3» 

MU  SMU(X.10..1C}.t> 

MU  SC4lt<V.6..1.l.l) 

MU  .i^TfirTioN  («t.i4.o.,M..vctei).y(iea>i 

MU  <l;.A0>,-U,10.,0..X(l(;t>,X(t»2>) 

MU 

MO  CONtINL-? 

MU  LI.' £(X.V,lM.l.tO,ICK4M) 

MM  C0«<T1^<l£ 

MLL  PLOTCIN) 

CKO  _ 

•uwouTiME  iNTEfmD(rt.Anpt,r:rci,rrc.Evn(,cufM.orci.OTc.O(m,OMt<.N 
10,0,  imo  I 


e- . — . IWTPWOBUUITIOH  OIStCIITlOM  (UMOUTIIC - 

e  TWH  MOMXfl  MLCULxiT*  TR:  iirrT’.'^nuMTiOM  BiSTrTric**  of  » 

e  CWSO  SVSTZd  I  •:■£  Iri  F-CC-"')  F  -.J  I  X.  'A  «f.Z  Cf.  ,  *Crt3)  S4W-T0- 
e  MM.  BioTC4Ti;d  IS  '•-'  L..!)  sv  I'-crr;.-.!  x  rjf  .;l  c;-..’»oo£» 
C  Of  Tuo  E.'j.u  f  .-LIT>, : ;  at  i;.)  ri  f  j  ?-•  fX.  : 

e  <m.lTUL£  Of  T>--  oif-:.5r<;-  f  xvxr  lo  t;- "d  f.  -a 
e  TO  THC  lAfur  siOMAL  TO  LEYiiMiht  T>i£  FencafT 


IMur  •  m  INTCOER  MMV  C0MT02N1N0  T>«  IIIIMV  OUTPI/T  Of  TtC  CVSO 
OCOOCR. 

owTRvr  •  •  ffxi  p»“»xy  coMTiiNjMo  Twf  TIP?  ruNCTxoM  ournrr  tvc 
Test  OIC^tAL  CL  .  ATC4  WU)  •fTt*  fROCMfXM,  TX  TXX  fUlCTIOH 
OUTPUT  Of  Tx 

MX  •  *  «AL  #’««y  ccwtaxhim  tx  output  fPCCTPxi  po*."^  ccppwprro 
Of  TVi  C  ;ec  «  CU  -  jT  e-T*-*  PHOCSOilNO  iV  TX  FAOT  fOLMlOl 
TMHOfOM  tffT) 

IX.  X.  CX  •  WOKIX  AMMYt  UMO  OV  TX  fPT  fUMOUTlX. 

fKOUC  •  A  WAL  AMAV  COPTAINIdO  TX  fBTQUDtCXCS  AT  UHICH  TX  fPT 
HAO  MLCULATCO  TX  frffCTXAL  COAPOXATO. 

<  AMPt  •  TX  APPLITUMf  Of  TX  TCOT  tlQHAU  IN  DM0. 

C  ff  •  TX  OANPU  NATC  IN  DM. 

C  PCMC  •  TX  PCX  IMUXf  Of  TX  TCDT  OXOMM.  COHPOXOTO. 

c  •  *0*  OP  Pftfix  ON  TX  kfjit  Tvm  irro  cf  waniaslc  inoicatti 
C  tx  UARIADU  U  UUI  DT  CITXI  tx  CCCOBCP  op  CNCOECR.RCSPeC* 

C  TlXlV 

«  ret.  fci,  fC3  •  Pou-ofT  fo'oufNoxu  fo«  tx  poxxxpu  intccmtoo 
C  In  TX  CUSD  CnCCCU  and  kccur. 

C  TO  •  TX  CONPAratX  tPCSD  Of  TX  OVUAIIC  XNTCOPNTOII  XN  OCC. 

C  MTIO  ■  TX  NNKINUN  OTCP  iXZt  TO  NXNlAtPI  BTCP  fIS  NATIO  IN  Cl 

e  f  •  A  fXOUCNCV  UAPIADU  USCD  TO  ffTTWlX  TX  TMT  fICNAL  CPtC- 
0  TML  CC.'..>0XI(T  AND  ALSO  TO  DCTCRNXX  TX  HARMMIC  fPCCTAAL 

0  CONPOXNTf. 

e  Rif  •  TX  OUTPUT  SPCCTKAI.  COWPOXNT  POUCR  IN  DM. 

C  BUN  •  TX  DM  Of  TX  POXR  IN  TX  TffT  SXQNAl  COXOXNTS. 
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e  IV  ■  Tw  NURKR  or  riLTDi  coerricitMT* 

C  t  •  AM  MMV  COHTAININQ  TIC  IMTUT  MO  OUTPUT  PltTEP  COCPF1CIO<TS 

aMMuiHttniiutro  nnwouriict  uses  iminnnnii«»mnn.«t««i 
e  DICOOCI  ■  TMC  evto  CNCOOOl 

e  Kcocci  -  nc  cun  oecoooi 

e  rLTmn  ■  cdcmtoii  or  tic  output  filtci  cocpFXCiDrr* 

riLTCR  •  riLTCRS  TIC  OUTPUT  UOXNO  TIC  COCFPXCXOfrO  ppoeuco  tv 
riTMCN 

e  tzoNM.  •  nc  rerr  txoiMi  oocmtoc 
e  rrm  •  tic  pmt  pourici  tumcporm  pmm  tic  inti  lxmmtv 


t - RIMOUTXIC  TMIT'--  . . .  .  -  - . 

C - IMITIALXZE  UMIMUS  MIB  MHMVt 

DtfCNSION  tNPUTtS9M),0UTPUT<SMB),P«X(lM> 
t.IUK(£ai.UK(l!:4).F'>^aX(ciMI. 

R^L  IKOO 

COMPLEX  evKOet)  _ 

AlOBflR)  •  SC;iT  (IB.  n  <(CtfN  -  4.1/lR.)  t  .Ml  B  Mt.l  t  MRTd.) 
c-»-  orrtRMXic  pcmc  uruc  or  Tcrr  txoMRL 


PCRK  ■  A  (AMPl) 

C -  OOCRATC  XIVUT  TIIC  FUCTXON 

CAU  •XCNALC OUTPUT.  64M.  Ft.  7SA..  IMt..  PEAK,  PEAK) 
e~—  PROCESS  nc  TXIC  FUNCTXOH  TMIOUOM  tic  cuss  tVSTER 


CALI  tHCOCEl<OUTPUT.IIWT.C-’.M.FS.EFC».FC3,FC3.ETC,EUre<.EUnH.DC) 
CAU  BCCOAEX < X«PUT. OUTPUT, M4,F0, DFCl,FCi,FC3,BrC, OOfU, CWtfl.BCI 

e —  rxLTER  nc  output 

CAU  rXLTERCOUTPUT.SMR.MP.f 

c —  turrwcT  nc  aucragc  uauc  from  nc  sampu  tTRXw 

MR  «  t. 

M  S  I  >  1.4»3« 

SUM  •  SUM  ♦  OUTPUTCI) 

•  CONTINUE 

AUCR  •  SUM  / 

DO  7  X  •  1.4 _ 

CUTPUT(I)  •  OUTPUTCXI  -  AUER 
T  CONTINUE 

e— —  CALCULATE  THE  tPCCTIIAl  COMPOMEMTt  OF  TIC  OUTPUT 

CALL  rTFPt<OUTPUr.OUR.4«M.S«.R,PtX,OUM.SUM.XUK.WC.CUC,ZER> 

C  BETERMXIC  THI  COMPONENT  AT  TIC  TEST  FREOUEMCV  AMB  CALCULATE 

C  TIC  RMS  UOLTACE. 


B 


C.— 


It 


Ml  ■ 

BO  t  K  •  t.l39 

F  •  <<K-l>A.Jt.)  t  Ft 

FRE«JE(K-1)  •  F 

PSX(K-t)  •  PSXIKI  ■  _ 

IF  tiF  .NC.  7:1.)  .MO.  IF  .HE. 

SUM  <  SUM  *  Si.'T  I  P-XIKI  I 

REF  •  lA.  S  ALOOIBIPSXIKI  ^  SM 

CONTXNUE 


IMR.I)  00  TO  t 
/.Ml) 


CAUULATE  XNTERHOBUIATXON  BXtTORTIOM 


tXF 

to  If  X  «  1,133  „ 

IF  (FREOUEiXi  .1C.  OXFI  00  TO  It 
tXt  ■  PSXIXI 

eONTXHL*  „ 

XMOO  «  SORT  (tXtl  ^  tin  t  IM. 

IF  IIMOB  .OT.  IM.  I  inOP  •  IM. 

RETURN 

END 
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APPEflDlX  R 


Sip,nal-to-’'’oi  Ii'ilio  Fro'  r  ii'. 

MogMA  $»•«•«(  iNPuT.oumiT,Tw>tt.oumrr> 

- - —fUTumcB  $m  wiogwH  .  , 

C  THIS  PROCftW  WASURtS  T>€  SlC^♦At-TO-HOTSE  PtSFORHANCC  OT  A  COSD 
C  tNCOOER  AWD  D£C0:LJ>  CC^‘^£Cr;B  '•,^-TO-T«C>:.  A  5V  ;LI 
C  site  UAME  IS  I.'.  UT  TO  Trt  SVSTt;^  AM  THE  DIFFEREHCE  KTUCEN  T>C 
C  OUTPUT  MiO  INPUT  CCrVUTED. 

e 

C  A  JIAXMALLV  FUtr  LI^•5«  PHASE  riR  PIITCT  IS  PLACE3  CN  T»C  OUTPUT 
C  Of  the  DECODER  TO  KENOUE  SIONAL  COTIPOfCNTS  ASOUE  3600  tQ. 

C 

e<m(ttt*SSU«tttXXXtSt  UARIAtLCS  . . . 

C  fREOl  •  THE  TEST  SIGNAL  FREOUENCV 
C  SMt  ■  TIC  rCASURED  SHt  UALUE 

e  TSIN  •  AN  ARRAV  CONTAINIHO  THE  INPUT  TINE  SERIES  SMVtES. 

TSOUT  •  «»AV  COMTAIHTHO  FIEST  T^r  C^COTER  OUTPUT  T:.»E  SERIES 
SANPLES,  THEN  THE  OUTPUT  TlfS  SERIES  SAS^'LES  CP  THE  SIR  PiLTtU. 

cm  •  AN  ARRAV  CONTAINING  THE  OlFTERENCE  BETUEEN  TIC  OUTPUT  SAft- 
PUS  AND  THE  INPUT  SA.  J>t£S. 

•  •  /)N  ARRAV  COKTAININQ  THE  FILTER  CCEFf  tCIENTS. 

■INOUT  •  AN  ARRAV  CONTAININQ  TrC  IlNARV  OUTPUT  Of  TIC  CUSS  ENCODER 

AHPl  •  THE  ANPLirUDC  Of  THE  TEST  SIGNAL  IN  DMO. 

rs  «  nc  SAMPLE  RATE. 

FCl.  fC2.  FC3  •  THE  ROLL-Off  FREOICNCIES  Of  THE  PRINARV  IHTCGRA* 
TORS. 

TC  •  THE  TIME  CONSTANT  Of  THE  SVLLAIIC  FILTERS. 

UNAM  E  <JNIN  •  THE  MAXIMUM  AND  MINIMUM  INPUTS  TO  THE  SVUMIC  fll> 
TER. 

tCTA  •  THE  HORMALlreS  3  TS  FREOUEHCV  OF  T»€  OUTPUT  HlTtR.  TIC 
FREOieiCV  IS  KCP..VILIZED  TO  THE  SAMPU  RATE. 

OAMMA  •  TIC  HCWTALirrS  WIDTH  C*"  TV^  POLL-Ofp  RSOICN  Of  THE  OUTPUT 
FILTER.  V-X  r  r.l  -.'i  IS  TIC  PREOUENCV  SAND  ICTT.CEH  TIC  BS«  AM 
R  OUTPUT  ANPLITUTES. 

PEAKS  •  TIC  MAXITVM  AMPLITUDE  Of  THE  TEST  tIONAL  IN  UOLTt. 

HP  •  THE  NUnSER  Of  FILTER  COEffICICNTt. 

OC  •  The  DUTY  CYCLE  OF  THE  SLOPE  O^RLCAS  CCTECTOR. 


ittitit  <«»!««(  t«(t((t((Tttzt«tnxmn>nst3ntnxnnnttsnssmxtzss 

FLTRGCN  •  THE  SUROUTINC  THAT  GENERATES  TIC  OUTPUT  FILTER  COEFFI¬ 
CIENTS. 

signal  •  TIC  TEST  SIGNAL  OEICRATCR.  PRODUCES  SAMLCS  OF  tlrXJSOI- 
OAL  WAUES  WITH  AT  POST  TWO  Pf.'cOUENCV  COT^EJfTS. 

ENCODES  •  TIC  CUSD  CNCOCSR  SLTOOUTTIC  WITH  A  SINGLE  ROU-OFF 
FREOUCNCY  IN  nc  PRl.TARV  ItirECRNTCR. 

OCCODCS  •  THE  CUSS  CEC03ING  SL'SrOUTXHE  WITH  A  SINGLE  ROU-OFF 
FREQUENCY  IN  THE  PRIIWRY  IHTECRATCR. 

FILTER  •  THE  SUBROUTINE  THAT  FILTERS  DC  INPUT  TIRE  SERIES  SAFPLCS 
USING  THE  FILTER  COEFFICIENTS  GENERATED  »V  FLTRGEN. 

POUCR  •  A  ROUTINE  TO  CALCULATE  THE  POWER  IN  A  SAMPLED  TUC  SERIES 
WITH  IMPEDCNCC  «  600  OMIS. 


CTOxmxttxtwxxiaxiniwxxtxnxwmRMMamxxxmxxnnxmmm 
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c - rtttn . . .  .  . .  ■  - - 

C -  IHITIMLIZE  UMIAILES  AM)  AAMV* 

D1REN5I0N  TS:h(S4M),  TSOUT(M««>. 
i.tiaco) 

tHTtCOi  IlNOirDStM) 

«(Din«)  •  MRTda.  »<(Din«  -a.vi*.  )  t  .mi  *  m«.)  t  soirrd.) 

c - imrr  vomiNa  vMftmss 

iKAO  s.  rs 

ICAO  t.rrci,  rrc,  eoatio 

RCAO  t,L'CJ,tTc.::.T:o 

CALL  W/'<C^T(f'.  -x, fS.tFCl. ETC, riATlO) 

CALL  Wit.'.Xo.»T<:-.  x,./i,>«.fS.L#'Cl,DTC,UlAT10> 

AEAO  l,l£TA,  C  ..)A 
KAXl  •  A(-<  a. ) 

MIKT  «.  •  TEST  AT  ET*  KZ  AM  SAaVtJt  RATE  •  •,rf 

MIWT  *,  •  EECl  •  •.'■-Cl.*.  K-t  •  C'^TIO  •  '.'.^TIO 

MIKT  ».  •  ffci  •  *.: -cs. *,  r-c  •  ‘..Tc.*.  :  .^rio  •  *,: '^no 

MINT  f.  •  flLTE.II  A>*«AICTtS«  AAC.  »£TA  •  ‘.tTrA.".  QAHnA  -  '.QAftl 


©~—  OOCRAIT  Ol/mff  flLTEA  COCTflClEWTa 
CALL  rLTIKeA(ttTA.aAAnA.»#.t) 

0 - <iCAC*ATE  infur  Tine  series  SArsfLES 


call  ;iCNAL<TSIA.SM«.fS.SM..R.  .fZAEl.R.  > 

f— —  AROCCSS  TM  IHPUT  TtflE  SERIES  T>«OOQM  T>C  CWSS  SV5TEII 

CAU  CHCOMl(TS!H.iT*<luT.t'x^.fs,ErCl.rca.rir!.E“C.tjrt.r>*<.DC) 
CALL  CECOC£l(llN(Xrr.TSOi.T.i,w«),f<.tfCl.Ki,f£3.5:c..V<)»«.OW^.DC) 


c - riLTER  ne  oimvr  cr  nc  oeccbcr 

CALL  flLTERlTSOUT.SMR.KR.i) 

C-~-  OELAV  THE  INAUT  SICAAL  START  TO  COR5JCSROMO  TO  T« 
C  flUTERREO  OUTIVT. 


M 


DO  M  10  •  1.4R9S 
«  •  2TS  ♦  10 
TSiNdO)  .  TSIH(lCO) 
CONTIME 


0—  ADJUST  OUTPUT  ARPLITUCC  SO  INPUT  POWER  •  OUTPUT  POWER 


CALL  POUERtTSiK. fS.PIN) 
CALL  P«>ER(TSC'..T.-*.  .-i.fS.POUT) 
OAIN  •  SwlT  IPIM/?0UT) 

00  4S  1  •  1,<  .i 
TSOUTin  •  TSCUTd)  S  CAIN 
M  CONTIMUE 

e -  CALCULATE  TIC  NOISE  POWER 


M 


M  M  I  •  l.ATSS 

ERR(H  •  TSOUT(I)  -  TSINdl 

CONTINJE 

CALL  P0UCR(ERN,4MS.fS,CRRP> 


CALCULATE  TIC  S/n 


Ml  «  It.  t  AL041S  (PIN  /  ERRP) 


PRINT  TIC  RESULTS 

PRINT  S,*  THE  SIONAL  TO  NOISE  RATIO  • 
CONTinS 

09 
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APPEIIDIX  E 


*'*1  ‘  1  i  ,  .‘Ht  'it.  *''i.  'r/'i  I  :‘r:‘0’i«  'ic;/  -ro '' i'"; v 


e 

c 

c 

c 


W>OOWOW  *«1H«<lNWT.0UTWT,T««i»0UTWT,^VOT> 
- FILTWOCO  SMU  WOOHOW— . 


THIS  PPOC»»n  TV€  STOftPL-TO-^^OI'E  PCTPO»F«WC€  Of  •  CU«0 

CHcoocP  «N0  cc;-  -ctej  r.v:ic-TO-'x<. _ n  jlz  f-EouENCV 

SIfC  WOWE  IS  I*-  UT  "C  T  :  tvSTErt  ivto  T>it  DlFrCKEMCE  *£TW££H  THE 
OOTPCT  p«0  INPUT  COrlPUTEO. 

«  TWWIfWUV  FEPT  LIKtP  PHASE  PI»  FIETfP  19  PLPCEO  ON  TVC  OUTPUT 
Of  TVC  DECOCEP  TO  ^iT^DUE  SIO.-<AL  CClVOKENTS  A90UC  36##  M2. 

THE  PPOOPiMI  H  PEPEA-^O  TWCE  T1FE9  AS  TVS  UAUJC  CP  TVC  STEP  flZC 
RATIO  IS  CVMNOCO  PrtOn  32  TO  36  29  IH  STEPS  OF  8  C3. 


uARipauf  tttmttxtxxxxatxtasntxfttmtsm 

r*C01  •  ***  APPAV  CONTAINIHQ  TVS  FnSOL'EKIES  THAT  TVC  SNR  HA9  lEEH 
RCASURCO  AT.  ThE  RANCC  IS  J3#  h2  TO  330#  H2. 

SNR  •  AN  ARRAY  CONTAININO  TVC  REASUREO  ST«  UAUCS. 

TSIH  •  tm  <«RAV  CONTAININO  TVC  INPUT  TIT£  FUNCTION  SAfPLES. 

TSOUT  •  AN  «N>AV  COtfTAlNTNO  F1?ST  TVE  tSCCtER  OUTPUT  TIRE  FUNCTION 
SANPLCS,  TVC  TVE  fl'.TER  OUTPUT  TITS  S"..pUS. 

ERR  •  AN  A?7RAV  CONTAININO  TVS  DIFFERENCE  lETVlEEH  TVC  OUTPUT  SAN- 
PLES  TVC  INPUT  SA/VUS. 

I  .  Mi  ARRAY  CONTAININO  TVC  FILTER  COEFFICIENTS. 

SlHOUr  •  Mi  <«»RAY  CONTAININO  THE  SIHARY  OUTPUT  OF  TVC  CVSB  ENCODER 
ARPl  •  the  AfTLITUDE  Of  THE  TEST  SIGNAL  IH  OSRS. 

FS  •  TVC  SAflPU  RATE. 

rcit  rca.  rc3  •  the  Roii-ofr  freouckies  of  tvc  prirary  Integra- 

rOiWI# 

t  TC  •  TVC  TIRE  CONSTANT  OF  TVC  SVLLASIC  FILTERS. 

C  UWW  6  URIN  •  TVC  RAXinUN  VMD  RINIRUN  INPUTS  TO  TVC  SVILA9IC  FIl- 

C  TER. 

C  STTA  •  TVS  HOTIWLITED  3  Cl  fPEOLCNCV  OF  TVS  OUTPUT  FILTER.  TVC 
€  FRiaUENCY  IS  NJ..nALI2E0  TO  TVS  &AAPLE  RATE. 

e  nnnnn  •  tvc  NC"fviLiztD  uicth  of  tvs  poll-off  psoion  of  tvc  output 

C  FILTER.  Tvl  >  :.S(.N  IS  TvC  FPEOtCNCV  lAHO  SLTVEEN  TVC  9S»  ANS 

C  S«  OUTPUT  AHPLJTU2C9. 

c  PEAK!  •  TVC  fWIAUN  AV^ITUT*  OF  TVC  TEST  tIONAL  IH  UOETS. 
e  vr  •  TVC  PtJHKR  OF  filter  coefficions. 

e  KM  •  TVC  NUNSER  OF  TEST  FRCQUCNCIES. 

e  K  >  TVC  BUTY  CYCU  OF  TVC  SLOFl  OUCRLOAO  DETECTOR. 


e 

c 

c 

c 

c 

e 

e 

c 

c 

c 

e 

c 

c 

c 


SuiNOUTlrCS  UViO  SS(tU*UU9UI(SU>Sm<t«»XXXn 

FLTRGEN  •  TVC  SUSROUTIVC  THAT  OOCRATES  THE  OUTPUT  FILTER  COEFFI- 
••  CIENTS. 

C  PLOT,  SCALE.  AXIS.  RCCT,  LINE.  PLOTf.  •  CALCOW  PLOTTINQ  ROUTIVCS. 

StCNAL  •  TVC  TEST  SIGNAL  CINtRATOP.  PROOUCSS  SANPLEl  OF  flHUSOI* 
DAL  UAUCS  UlTM  AT  ROST  TWO  FRiOUENCV  CCNPOfCNTS. 

CllCOOCl  •  TvC  CU9D  ENCOOSR  S'.'3RCUTIN£  WITH  A  SIN6U  ROLL-OFF 
FREOUENCY  IN  THE  PHirARV  lhT£i»»TOR. 

KCODCI  •  TVC  CVSD  DCCOSINO  SUTPOUTINC  UITH  A  SIHOU  ROLL-OFF 
FREQUENCY  IN  THE  PRIRFRT  INItw^iVTOR. 

FILTER  •  TVC  SUIROUTINE  TVVRP  FILTTPS  TVE  INPUT  TIRE  FUNTION  SAH- 
PLCS  USING  THE  FILTER  CCEFFICILNTS  GYRATED  SY  FLTKCCN. 

e  FOUCR  •  A  routine  TO  CALCULATE  TvC  POWER  IN  A  SAfVlCO  TIRE  FUNC- 
e  TION  UITH  IRNCDENCC  •  Gw#  OHRS. 


1  :iA 


c - PROORM  rrpwT - 

C -  INITIALIZE  UAAIAIUS  AND  ARRAVS 

DINCNSION  nrOKiet),  SMKIM).  T$IN(6«M>.  TSOUTIEMt).  DOXMM) 
1,*(R?9).  C.A:h(£’9) 

IHTEr.£R 

A(Pfcn«)  •  soMTdA.  ttuttm  s  .Ml  a  6M.)  a  soirro.) 

lOM)  •  -1 

C -  iwvr  AND  MINT  UORXINO  UAfflABLEt 

REM)  a.  AMI,  FS 
HEAD  a.FCl,  TC 
R£A0  a, BETA,  CAMA 
PCAKl  •  A(A."F>n 

MINT  a.  •  TEST  AT  •,AM1.*  DBNa  AHJ  ‘.Fi,*  IP»* 

MINT  a,  •  WITH  TC  •  '.TC 

PRINT  a.  •  aCTA  •  META.*.  OANNA  •  •.OANNA 

C -  OOCRATE  FILTER  COEFFICIENTS 

CAU  FLTR6eN(irrA.0A«*M,NP.a) 

C -  INITIAL2E  PLOTTER 

CAtX  FACTORT.S) 

CALL  PLOTta..  a..  -3) 

e -  START  OF  tIONAL-TO-NOltC  LOOP 

DO  ISM  NTiNEs  •  a.a.a 
KN  •  0 

ICV^)R  •  ICHAR  «  1 

RATIO  •  3a.  ♦  HTirSS 

.'mLL  WV».'<OAT<unAX,yniN.FS.FCl.TC.RATI01 

DO  3M  F  >  3a«.36M.lM 

FN  •  FN  *  1 

C-~-  OCfCRATE  TEST  SIQNAL  FREOUCNCV 
FRCOllFN)  •  F 

C- - C£»€RATC  INPUT  TIRE  FUNCTION  8ARPLCS 

CALL  SI<iNAL<TSIN,S0M,Ft.FVE01(FN).S.,PEAKl,S. ) 

C -  PROCESS  the  input  TIFK  FUNCTION  TKROOCH  TIC  CVSD  CVSTER 

CALL  ENC0^61(TSIM,aIN0UT,5■>^».FS.FCl,FC2.FC3,TC,WV'X,yHIN,DC) 

CAU  D£COC£l(aiNCUT,TSOUT,SOM,FS,FCl,FCa,FC3,TC.UaAX,W1IH,CC> 

C -  FILTER  THE  OUTPUT  OF  THE  DECODER 

CAU  F1LTER(TS0UT.EM«,NP,I> 

C -  DEIAV  THE  INPUT  SIGNAL  START  TO  CORRESPOND  TO  TIC 

C  FILTERRED  OUTPUT. 

DO  3R  ID  •  1.4«M 
KD  •  ZM  «  ID 
TSIN(ID)  •  TSIN(KD) 

M  CONTINJE 

C -  ADJUST  OUTPUT  AMPLITUDE  SO  INPUT  POVER  •  OUTPUT  POUER 

CALL  P0LCR(TSIN,4a06.FS.PINI 
CALL  POLES <TSXrr,4c,s,FS, POUT) 

OAiNicN)  •  scir  (PiN/i>oun 

DO  40  I  •  t,4,>cs 

TsouT(i)  •  Tsourai  a  gain(kn) 

Ot  CONTINUE 

e——  CALCUATE  THE  NOISE  POWER 

DO  M  I  •  I.40M 

ERRCI)  •  TSOUTdl  -  rSIN(t) 

M  CONTINUE 

CAU  P0UER(ERR.40DS,PS.ERRP> 

C -  CALCULATE  THE  S^ 

SMKKMI  ■  IS.  t  ALOOlO  (PIN  /  ERRPl 
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PLOT  THC  RCSULT* 


M* 


CCmiNUE 

IF  (10  .-9  .QT.  •>  00  TO  CM 
MLL  PL0T(ia..«.,-3) 

CPU  i':;cr.L<F-'3i,n..«x» 

CPU  I 

CPU  U.MIS(«.,«..14^KCdUC7CV  (KS>, 
|N*2 1 1 

CPU  p.yif(«.,«..SH$m  (Ci).t.6.,M., 
CpU  P£cr(0..«..6..1«..»..i) 

COMTIHO? 


-14.1*.,0..F1«ai(DHl).FREQl(K 

MdnPD.SMRCOfft)) 


CpU  U<.INe(FREai,S*«.KN.l«,lCHPA,-l) 

continue 

CPLL  PLOTE(N) 

CNO 


i 
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PipTatchort  Si  I  n.-i  l-io-i'loi  c-o  ,  I,, r  , , >  r-  j 

fWOCHAH  OUTPUT, T<«**eUTfUT.PtOT> 


C- . .  — . mSfWTCWO  6HH  PKOOHm - - 

C  THIS  PROCRftH  flEAfuWS  Tt-S  SICNftU-TO-HOISS  PESTORflPHCt  OF  *  CUSO 
0  tNCOCCR  AMD  d;c3:;r  cc.'-'-:cted  j-'Csc-To-t^x.  «  kciouencv 

C  SINE  UPUE  IS  I.-r'-JT  TO  :  SVSTEn  AMB  THE  DIFFESEHCE  EETliEEH  THE 

0  OUTPUT  MIB  input  CCPUnjTEJ). 

<  THE  CPLCULPTIONS  APE  PEOrOftnEO  PS  THE  UPUJE  OF  THE  ENCOCCH  STEP 
C  SIZE  RATIO  13  »:i:?T  ca-  ;T.->fT  pho  r,-t  cecc;.r  step  si;r  p^rio  ts 

C  CHANCED.  T>:£  EnC5D2:.<  Ri^TIO  IS  P2  £».  UHILE  THE  EtCUj^  IS  STEPPED 

C  FROn  32  TO  32  Dt  IH  STEPS  OF  2  D3. 

C 

C  P  nAXinALLV  FLAT  LIKEPR  PHASE  FIR  FIITER  IS  PtPCED  SH  THE  OUTPUT 
C  OF  THE  DECOtSR  TO  KE-WE  6ICNAL  CC.TP0NEHT8  AJCVE  3GC«  K2. 

C 

cttatxttttttttxsxtixxxtt  UARIAILES  txxtxxxtxxxxxtttxxatxtatxxxxitttjx 

C  FREOl  •  AM  P-'^AV  CCMTAIHtNO  THE  ntOUENCIES  THAT  THE  SMR  HAS  »EEM 
C  MEASURED  AT.  THE  RANCE  IS  ECO  KZ  TO  3.'>3»  HZ. 

C  SMt  •  AN  ARRAV  COMTAININO  THE  MEASURED  SNR  UAUJCS. 

C  TSIH  •  AN  ARRAY  COMTAININO  THE  INPUT  TIFZ  SERIES  SANPIES. 

C  TSOUT  ■  AN  <«RAV  COMTAININO  FIRST  THE  DECCIHR  OUTPUT  TirE  SERIES 
C  SAMPUS.  TV3*  THE  OUTPUT  TINE  SERIES  SAiPLES  CF  THE  FIR  FILTER. 

e  DW  •  AN  ARRAV  CONTAINING  THE  DIFFERENCE  BETWEEN  THE  OUTPUT  8AN- 
C  PES  AND  TIC  INPUT  SA.NPLES. 

C  •  •  AN  ARRAY  COMTAININO  T>C  FILTER  CCEFFICIENTS. 

C  TIF6  •  AN  ARRAY  CONTAININO  THE  TIPE  THAT  TVE  FIRST  2M  SAHM.ES 
C  ARE  TAICEN  SO  THAT  THEY  MAY  BE  PLOTTED. 

€  BtNOUT  •  AN  MiRAV  COMTAININO  THE  BINARY  OUTPUT  OF  THE  CUSD  ENCODER 

«  MPl  •  T>C  AMPLITUDE  OF  THE  TEST  SIGNAL  IN  OBNB. 

€  Ff  ■  THE  SANPU  RATE. 

C  rci,  Fca.  FC3  •  THE  ROLL-OFF  FREQUENCIES  OF  THE  PRIMARY  IMTEORA- 
C  TORS. 

C  TC  •  THE  TIME  CONSTANT  OF  THE  SYLLABIC  FILTERS. 

C  UMAX  %  VMIN  •  T»C  mxiMUM  AND  MINIHUH  INPUTS  TO  T>C  SYLLABIC  FIL- 

€  TER. 

€  BETA  •  THE  NORMALIZED  3  D»  FREQUENCY  OF  THE  OUTPUT  FIlTEJt.  TVC 
e  FREOUCNCY  IS  NOHPALIZED  TO  THE  SAMPLE  RATE. 

C  OAHMA  •  THE  NORMALIZED  WIDTH  CF  TTS  ROLL-CFP  REGION  CF  T»«  OUTPUT 
C  FILTER.  THE  .  It  THE  FREQUENCY  BAND  BETWEEN  TVC  SS«  AND 

e  »  OUTPUT  AfiPHTUEES. 


C  f#  •  THE  NUMBER  OF  FILTER  COCFFICltNTB. 

O  EN  •  THE  NUMBER  OF  TEST  FRCOUENCICS. 

C  OC  •  TvC  DUTY  CYCLE  OF  THE  SLOPE  OUERIOAD  DETECTOR. 


CBttBtttKttKttttttMdtxtBtsxsxnsztnxnsstuxssBxxnTttnttnxxTfvrt 
ixtxtutxtxtttttxtxxx  suiPcuTiNES  USED  txxixixmtxxxxxziitxzxxxtxxis.xxt 

C  FLTRCEN  •  THE  SUBROUTINE  THAT  GENERATES  TVC  OUTPUT  FILTER  COCFFI- 
e  CIENTS. 

C  PLOT,  SCALE,  AXIS,  RECT,  LINE,  PLOTE.  •  CALCOMP  PIOTTINO  ROUTINES. 

C  SIGNAL  •  THE  TEST  SIGNAL  CEHERATOR.  PROOtCES  SMTPIEB  OF  SIHUSOI- 

C  DAL  MAUES  WITH  AT  MOST  TWO  FREQUENCY  CC»*>0N1NTS. 

C  EMCODCl  .  THE  CUSO  ENCCOER  SUSROUTITC  WITH  A  SINQU  ROLL-OFF 
C  FRCOUENCY  IN  TVC  PVilFiARV  INTEGRATOR. 


C  KCOKl  •  THE  CUSO  CECOOINO  SU5R0UTIVC  WITH  A  BINQU  ROLL-OFF 
C  FREQUENCY  IN  TVC  PRIMARY  IKTCCRATOR. 
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riLTCT  ■  THt  SU!■:?^lCTr^T  “HOT  r’lTTP?  Tvr  I»TVr  Tirs  SERIES  S<«PLES 
USINd  THE  FILTER  C-CCFFISltjns  iElERRTED  #Y  FLTREEN. 

POUER  •  A  ROUTlrr:  to  CTLCL'I-ATE  the  PCJER  1H  a  S^WX®  Tir.E  SESIES 
WITH  IKPECENCC  •  e««  ChflS. 

PlTlne  •  A  LIhEA«  PLOTTIKQ  ROUTINE  TO  PLOT  SIGNAC  (WPCITUDE  VS. 
TIRE. 

CXtxtxttttttxtrtTtttTttTtttTrrnstrtTtTTxtrttrttmxxsttrtfKtxitttittftt 

CttxxsxxtxtxxtxxxtxxxxxxxttxxxxxxxxxxxxxxttxtxxtTrxtifttfrfTrtf**xxx-fj^t 

C - PROCRAN  START— - 

C -  INITIALIZE  OARIAtCES  AND  ARRAYS 

OirCNSICN  FPEQIMM).  TSIN(5«et).  TSOUTtSAO*).  ERR>&«M) 

l.tiEZi),  TI'Kc'J),  CAINIHSa) 

IKTECER  BII.'C'JT(t;.‘ 

A(CS;t«>  •  SCRTCIA.  UKCSn*  -4.  Vl«.>  t  .Ml  t  6M. )  <  SORTIS. ) 
ICHAR  •  -1 


INPUT  UOtaCINO  UARIASLES 


READ  t.  AflPl.  FS 
REA®  t.FCl.  TC 
READ  *.Enn.  C/TKrtA 
PE«1  •  A(F.??n 

PRINT  t,  •  SM  TEST  AT  •.AAPl.*  ETTA  ANO  '.FS,* 
PRINT  t.  •  L'ITH  TC  •  '.TC.*  F.O  T-ITIO  -  ‘/^TTIO 
PRINT  1,  •  lETA  •  •,J£TA,*.  GASftA  •  •^QA.VM 


IPS* 


e -  GENERATE  OUTPUT  FILTER  COEFFICIENTS 


CALL  FlTRO£N<trrA.GANnA,W>.|> 
e -  INITIAL2I  PLOTTER 


CALL  FACTOR (.S» 

Call  plotiz..  a..  -31 
c— -  start  op  sicnal-to-noisc  loop 
00  too*  NTiNES  •  a.s.a 

KN  •  C 

Khar  •  ICHAR  ♦  1 

RATIO  •  39.  ♦  NTirrS 

CALL  U?W<CPT(Uf,AX,L^IN,F8,rCl,TC,RATW> 

IF  (ICKAS  .GT.  #5  CO  TO  3 
Eimx  •  WAX 

CUHN  •  WIN 
3  CONTINUE 

00  3M  K  •  3««.36M.1M 
KN  •  KN  A  1 

C -  GENERATE  TEST  8IOWL  FREQUENCY 

FREOKKN)  •  K 

C -  GENERATE  irW  TINE  SERIES  SANPUS 

CALL  SICNAL(TSIN.SMS.FS,FRE01(KN),t..PEAKl,S.) 

C -  PROCESS  THE  IfSVr  TINE  SCRIES  THROUGH  THE  CUSD  SYSTEN 

CALL  CNCOECl«TSIN,IINOUT,5^9».PS,FCt,FC?,FC3.TC.EUNX,EVNH,CC) 
CALL  D£COKl<IINOUT,TSOLrT.SCO*,FS,FCl,FC«,FC3,TC.,UflAX,UnlN,PC) 

C -  FILTER  THE  OUTPUT  OP  T«  DCCOCER 

CAU  FIlTeR(TS0UT.6MS,r«>.l> 

C -  OeiAV  THE  INPUT  SIGNAL  START  TO  COBRESPOtC  TO  T>C 

C  FILTERRED  OUTPUT. 

DO  3S  ID  •  1,4«9S 

KO  •  SOS  A  ID 

TSlNdD)  •  TSINIKB) 

30  CONTINUE 

C -  ADJUST  OUTPUT  AMPLITUDE  SO  INPUT  POOCR  •  OUTPUT  POUER 

CAU  FOWFRITSIH.ArJR.FS.PIHI 
CAU  POV.ER<TrCUT,4?56,FS,POUT) 

CAIHIFN)  •  S'  tr  (PINXf»CUT> 

DO  4*  I  •  1,- 

Tsourd)  •  rsouTd)  t  sainckni 
40  COKTlftJC 


CALCULATC  THE  NOISE  POUEX 


n 


DO  S*  I  •  1,4«« 

ERR(I)  •  TSOtfTd)  -  TSIN«1) 
CONTIKT 

CALL  FXlU£:il(E)n.4»£«,n,ER») 


0. - CALCIiLATE  THE  *>T< 


tNRiiCNI  •  IS.  t  AUMtS  (PIN  /  CPPP) 
0 -  PBINT  AND  PLOT  THE  RESULTS 


««« 


1*M 


URITE(6,€0«)  rREOKtHI.  V(R<rH>.  OAIH(ICH> 
F0BnAT(?0X,Fl*.I,i,X,fld.«,bX,Fl*,8) 

CONTIh’.>i 

IF  (KHAR  .CT.  8)  CO  TO  S-T4 
CALL  LC£CAL(Fr:C1.10..C'(» 

CALL  SCftL£(S^^»,6.,CN,t ) 

•  <KZ),-lA.l*.,*..FR€ai(IO4+l).F1JE01(K 

In^c  )  * 

fill!!  <E»>.«.«*,9«..6»«(IOffl).SI«(DH4)) 

CALL  FrCf(8,,S,,6.,Ii*.,i.,3) 

CONTIh-.t 

CALL  LCLIHt(FR£ai.SM?.ICH,l*,ICHA«.-l) 

CONTIKUe 
CALL  FLOTE(N) 

END 
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APPF.riDIX  ’J 


Sir  nn  1 -to-iV)i y,-. ,  t.,.-  Po-’;-?;'  '’ro.-r,  !■' 

POOGIMN  5H*BMNPUT,OU'Wr,T*Pei*eOTPUT.PW) 

r-— . - . SNP  WS.  INPUT  P0U»  P«OQ<Mn - 

C  THIS  PPO(a»«1  t*?:A5t'ReS  TH?  STC.NAL-TO-NO!$E  PtaFCRN<i*iCt  OP  P  cwo 
CNCOOEP  f>ND  CiCC':/?  COf<.-f:C7£Q  i'AClC*TO-S'<ciC.  A  'U  P;?ECU£.SCY 
C  5lr«  UflUC  IS  I.---UT  TO  T-  •.  SvSTtfl  PhD  Th£  DlPPtKEi«£  lETWcEN  THE 
0  OUTPUT  PhD  if^ur  COPHUTED. 

C  A  fWXIIMLLV  PLAT  LIht;«  PWPSC  PIR  PILTER  IS  PIACED  CH  TNC  OUTPUT 
C  OP  the  DECODEJl  TO  REhOWC  Sl&hAt  CCKrt»CNTS  P»OU£  HZ. 

C  THE  INPUT  PO'.'ER  IS  W*niED  PHON  -4#  D?h#  TO  •  rt*9.  TVE  P«OC»PP1 

C  IS  REPEATED  AS  THE  STEP  SIZE  RATIO  IS  VARIED  PRCA  32  Ct  TO  3o  D». 

C  r>C  RESULTS  ARE  T>CN  PLOTTEO  ON  A  SINOU  CRAPH. 


custttttrttxstttwtxtxi  uwiAtus  xstxxaxxxttsxmxxtxmxtxttxxxxxxx 

C  PREOl  •  an  A/IRAV  CCNTAtHlha  THE  PREOUEHCIES  THAT  THE  SNR  HAS  KEN 
C  NEASUREO  AT.  THE  kfinU  IS  3s4  hZ  TO  HZ. 

C  sm  •  m  ARRAV  CONTAININO  THE  HEASURED  SNR  UAUiES. 

C  TSIH  •  AN  ARRAY  CONTAININO  THE  INPUT  TINE  SERIES  SANPIES. 

C  TSOUT  •  ARRAV  CONTAININO  PIRST  TVT  DECOCER  OUTPUT  TINE  SERIES 
C  SAflPtES,  THEN  the  OUTPUT  TlPiE  SERIES  S.ViPLES  Or  THE  FIR  FILTER. 

C  ERR  •  ARRAY  CONTAININO  THE  DIFFERENCE  lETVEEN  TIC  OUTPUT  SAN- 
C  PUS  AND  the  INPUT  SA/Y=L£S. 

C  •  •  AN  ARRAV  CONTAININO  THE  FILTER  COEFFICIENTS, 

C  TINE  •  Ml  ARRAY  CONTAININO  THE  TINS  THAT  THE  FIRST  2M  SANPLE8 
C  ARE  TAKEN  SO  THAT  THEY  NAY  K  PLOTTED. 

e  IINOUT  •  AN  ARRAY  CONTAININO  THE  UNARY  OUTPUT  OF  THE  CUSD  ENCODER 

C  AHPI  •  THE  AflPLITUOE  Of  THE  TEST  SIGNAL  IN  DJ!». 

C  rs  •  THE  SANPU  RATE. 

C  PCI.  FC2,  FC3  •  THE  ROU-OFF  FREQUENCIES  OF  TW  PRIHARV  INTEORA- 
C  TORS. 

C  TO  •  T«  TINE  CONSTANT  OF  THE  SYILASIC  FILTERS. 

C  VNAX  S  OWN  •  THE  NAXINUN  AND  NIHINUN  INPUTS  TO  THE  SVLIAJIC  FIL- 

C  TER. 

C  seta  •  TIC  NIDPOINT  OF  THE  OUTPUT  FILTER  TRANSITION  IANS. 

e  OAAWA  •  the  NORIWLirSO  WIDTH  OF  THE  ROLL-OFF  REOION  OF  THE  OUTPUT 

e  FILTER.  ■!>£  R£3XCN  IS  THE  FREOUENCY  lAHO  lETUEEN  THE  C5»  AND 

C  S»  OUTPUT  ANPLITLDES. 

c  PEAKi  •  TIC  PVMinun  anplitudc  or  the  test  signal  in  uolts. 

C  NP  •  tic  nunkr  of  filter  cccfficicnts. 

c  oc  •  the  duty  cycu  of  tic  slope  ouerload  dctector. 

<:tltlXllllXXXXXXXXXXX  SUSROUTItSES  Us£0  XXXXXXXXiXXJXXXXXJlXXXXXXXllXXXXXS 

C  FITROEN  •  THE  SUIROUTINC  THAT  GENERATES  THE  OUTPUT  FILTER  CCEFFI- 
C  CIENTS. 

V  PLOT.  SCAU.  AXIS,  RECT,  LINE.  PLOTE,  •  CAICOAP  PtOTTINO  ROUTINES. 

C  SIGNAL  •  nC  TEST  SIGNAL  CDfSRATOR.  PRODUCES  SANPLES  OF  Sm*50I- 

C  DAL  UAUES  WITH  AT  NCST  rjO  FREQUENCY  CONPONEJfTS. 

C  ENCODEl  •  THE  CVSO  ENCODER  SUIAOUTINC 

0  KCOKl  ■  THE  CUSD  DECODER  SU3ROUT' 

C  FILTER  •  THE  SU9ROUTINE  TWIT  Flir  THE  INPUT  TTTE  SERIES  SANPUS 
C  USING  THE  FILTER  CCEFFICIDTTS  CE/CRATED  lY  FLIVOth. 

C  POWER  •  A  routine  TO  CALCULATE  THE  POWER  IN  A  SAPPUD  TINE  SCRIES 
C  WITH  INPCCCMCC  •  6GS  CMtS. 

CtXXtXXXXiXXXXXXXXXXXXXtXXXXSXtXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXlXXX 
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c - pfioolwfl  *T<«rr - 

C -  INITIALIZE  UAfilASLES  A«0  ARRAY* 


DINENSION  AnPl(»*).  SNA(2««),  TSIH(E«M>.  TS0UT(5«M).  0!R(SM«> 
IHTTCER  IINOUT<S(?(»*) 

A(Dlri»>  •  S^lRTCie.  XXCDSN*  -4.  )/!•.)  t  .MI  S  fiM. )  S  S0RT(2. ) 
ICHAR  •  -1 

C -  INPUT  AND  PRINT  THE  UORKINO  VARIAILES 

READ  *.  FREOl.  FS 
READ  *,FC1.  TC 
READ  t.lETA.  C,Vr^ 

POINT  t.  •  SK-»  test  at  •.FSEOl,*  D5JM  AND  ‘.FS.’  IPS* 

PRINT  »,  *  WITH  TC  •  ‘.TC.*  r  0  R;UIO  •  ‘.O-^TIO 

PRINT  *.*  FILTER  PAOAAtTER*  ARE,  iCTA  •  '.I£TA. '.O^WA  •  '.aAmA 

0 -  OENEfMTE  OUTPUT  FILTER  CCEFFICIENT* 

CALL  FXTRC£H<irrA.CA«IA.NP,*) 

e -  INITIALZE  PLOTTER 

CALL  FACTOR (.S> 

CALL  Pl.OT(a..  a..  -3) 

0 -  START  LOOP 

M  IftM  NR  •  8. A. a 
RATIO  •  TJ.  ♦ 

CALL  UrMXOPT(UNAX.UNIN.FS,FCI.TC,RATIO) 

- ^TART  of  SI«NAL-TO-NOISt  LOOP 

DO  MO  F  •  l.tM 
ANPItE)  •  •4*.  4  ,4  <  K 
PCAKl  •  AIANPKK)) 


C— -  CCNCBATE  INPUT  TINE  SCRIM  SAHPU* 

CALL  SIONAL(TSlN.SMA,F*.FRCai.*..PEMCl,l.) 

C— ~  PROCESS  THE  INPUT  TINE  SERIES  THROUGH  THE  CV8B  SVSTEH 

CPU  CHC0DEl<TSIN,BINC«JT,ST**.FS,FCI,FCa.FC3,TC,l*WX,UniN,0C) 

CALL  CCCOI«l<>lhOUT,TSOUT,50va,FS,FCI,Fca,FC3,TC,W!AX,YNIN,eC) 

C - FILTER  THE  OUTPUT  OF  THE  DECOK* 

CALL  FILTCR(TSOUT.SM«.NP,l) 

C -  DELAY  the  input  SIGNAL  START  TO  CORRESPOND  TO  THE 

C  FILTERREO  OUTPUT, 

DO  3*  ID  •  1,4«9« 

KD  •  aa«  4  10 
TSIN(ID)  •  TSIN(KD) 

M  CONTINUE 

C -  ADJUST  THE  OUTPUT  SIQNAL  ANPLITUDC  SO.  OUTPUT  POWER  •  INPUT  POWER 

CALL  P0UCR<TS1N. 40^5*. FS, PIN) 

CALL  POL^R<TSOUT,-<.;^,FS,POUT) 

OAIN  •  SCOT  (PIN/POUT) 

DO  4*  I  •  i.4j;3 
TSOUT(I)  •  TSOUT(I)  t  GAIN 

4*  CONTINUE 

C -  CALCULATE  THE  NOISE  POWER 

DO  6«  I  •  1.4SCS 

ERRU)  •  TSOUT(I)  -  TStml) 

M  CONTINUE 

CALL  PCWER(ERR,4M*,FS.ERRP) 

C— —  CALCULATE  THE  S/N 

SHR(K)  •  1*.  S  ALOOt*  (PIN  /  CRRP) 
ir  (ICHAA  .LT.  S>  CO  TO  «»« 

IF  (SNR(K>  .LT.  SM)<l*i>)  SWCE)  ■  SNIdSt) 
tM  CONTINUE 

DM  CONTINUE 
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MIKT  MO  nor  T»C  RESULT* 


ICHM  •  tCMfA  ♦  1 
IF  (1CH!«  .CT.  •)  CO  TO  9«* 

CALL  SC«t.C<i: 'P1.10..U-4.1* 

CALL  scAU'(£.'^,s.,::«. n 

^CAU.  (O(28(«..*..l$KAAHLtTV0C  (09n*>.-l*,l*..*..M1Rl(l*l 

CALL  AyiS(«.,*..8HSNR  IDS ),*,«. ,8*., SMR(  1*1  ),*NR(1M)) 
CALL  RECrve..4..6..1«..«.,3) 

M*  CONTIHUE 

;aLL  LlNE(AM>t,SMR.lM.i.t*,lCHAR) 

IM*  CONTINUE 

CALL  PLOrE(N) 

END 
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APPKMDIX  V 

Mlsnatcl-.c:  -  -it  'o  v'~, ,  ln:  i 

Si  n;  i  ''ov.'er  ira  r  .-. 


mSI«0<INPUT.0irmiT.T<»»>«»0UTFUT.PWT) 

c - mSfWTCMEO  SN*  w*.  IMPUT  POUCa  PBOaiMH - 

r  THIS  PROCRAfl  ^A'SL.'PeS  TW5  SICK»U-TO-+«irt  P«F0‘?W«t  CP  P  cvso 

C  CKOD'R  PMO  t’.cc:  a  cck;  icted  rw:K-To-t;  ■;<.  a  ,l£  FE-iLVfCv 

0  SINE  UFiOE  IS  If-'UT  TO  T-:  SvSTEH  and  1>£  0IFTES£NC£  irrUEEN  TFC 
C  OyTfVT  AMO  INrUT  CO^?UTEO. 

C  A  JWX1NAI.LV  FLAT  LlNT'tR  PJ-A«F  pm  FILTTR  IS  PUCEO  CN  TVE  OUTPUT 
C  OF  THE  OECOCER  TO  SIGNAL  CCNACf^NTS  A20LE  3b0«  KZ. 

C  the  INPUT  POVER  IS  UFRIEO  FRCH  -4«  DTPS  TO  •  tT'9.  TVT  PROORPH 

e  IS  REPEATED  PS  TF'E  STEP  S1"E  PATIO  IS  UfRIEJ  F  .;,1  22  El  TO  ES  Of 

C  IN  the  decoder.  CHILE  THE  EhOODER  STEP  SIZE  RATIO  IS  FELD  CON- 
C  STAHT. 

C  THE  RESULTS  ARE  T>CN  PLOTTED  Ott  A  6IHSU  QRAFH. 

C 

cmmssts»tsntx»ttM  uariailes  ttMxxxxxtansttttaxxxaxntsxxts* 

C  FREOl  •  AN  ARRAV  COHTAIHINO  Tl-5  FriFOUEPOlES  T»'AT  TIC  SNR  HAS  SEEN 
C  nEASUHEO  AT.  THE  RAHCE  IS  3C«  HZ  TO  3.C»  HZ. 

C  SNR  •  AN  ARRAV  COtfTAlNINO  THE  fCASUREO  SNR  HALLES. 

C  TSIH  •  AN  ARRAV  COKTAIHIMO  THE  INPUT  TINE  SERIES  SANPLES. 

C  TSOUT  •  AH  ARRAV  CCNTAIHIN3  FIRST  TVE  CECCDCR  OUTPUT  TIPS  SERIES 
C  SAflPLES,  THEN  TV*  OUTPUT  TI.'E  SERIES  SF.-PLES  OF  THE  FIR  FILTER. 

C  ERR  •  AN  ARRAV  CONTAIHIMO  THE  DIFFEROXE  SETUEEN  TIC  OUTPUT  SMt- 
C  PUS  AND  THE  IWVr  SA-tPlES. 

C  t  •  AN  ARRAV  CONTAINING  THE  FILTER  CCEFFICIENTS. 

C  TINE  •  AN  ARRAV  CONTRIHINO  THE  TIPS  TVPT  THE  FIRST  RS*  SANPUS 
C  ARE  TAKEN  SO  THAT  THEY  NAV  SC  PLOTTED. 

C  tINOUT  «  AN  ARRAY  CONTAINING  TIC  DINARV  OUTPUT  OF  TIC  CUSD  ENCODER 

C  ANPI  •  THE  ANPLITUBE  OF  THE  TEST  SIGNAL  IN  DSNS. 

C  FS  •  THE  8AHPU  RATE. 

cn  NOTtl  A  ‘D*  OR  •€•  PTEFIX  ON  TFC  FOLICUII*!  UPRIA8U8  INDICATES 
«  THAT  THE  W'-RIPIU  IS  U.'TO  IV  THE  TtCorCR  CR  ENCC'ER.RTIPSCTIUE- 

<  LV.  UNAX  AND  LAIN  RAY  IE  CONTRACTED  TO  UNX  AND  U>V1,  RESPCCTIUE- 

«  LY. 

C  FCl.  FC2.  FC3  •  THE  ROLL-OFF  FREQUENCIES  OF  THE  PRINAAV  INTCGAA- 
C  TORS. 

C  TO  •  THE  TINE  CONSTANT  OF  THE  SVlLAStC  FILTERS. 

C  VHAK  S  UnXN  ■  THE  NAXXflUN  AND  HININUN  INPUTS  TO  THE  SVllASIC  FXL- 

e  TER. 

e  SETA  •  nc  NORNALIZEO  HIDPOINT  OF  THE  TRANSITION  SAND  OF  THE  OUT- 
e  PUT  LOy  PASS  FILTER. 

C  GANNA  •  THE  NOPNALIZEO  WIDTH  OF  THE  PCLl-OFF  RraiCN  OF  THE  OUTPUT 
FILTER.  The  FEGICH  IS  THE  FREOUENCV  SAND  lETCCEN  THE  8S«  AND 

C  5%  OUTPUT  AHPLITUOES. 

0  PCAKl  •  THE  NAXINUN  ANPLITUBE  OF  THE  TEST  SIGNAL  IN  UCLTS. 

C  NP  •  TIC  NUnSER  OF  FILTER  COEFFICIENTS. 

C  KN  •  nC  NUNSER  OF  TEST  FRCOUENCICS. 

C  DC  •  THE  DUTY  CVCU  OF  TIC  SLOPE  OVERLOAD  CCTECTOR. 

Cttttittt«t*t(ssnts<XK«tu»ussx»n»xsnn»nst3»nsxnsttttsx»t 

CUSSJtSIIIISISSSISXt  SLZ.'^OUTIJCS  USED  SSII»XU»XXI»UXZIUX»»SUXX 

C  FLTRCEN  «  THE  SUBROUTINE  THAT  GENERATES  THE  OUTPUT  FILTER  COCFFI- 
C  CIENTS. 

C  PLOT.  SCAU,  AXIS.  RECT.  LINE,  PLOTE.  •  CALCC!»  PLOTTING  ROUTIKLS. 

C  SIGNAL  •  TIC  TEST  SIGNAL  GCNERATOR.  PRODUCES  SANPUS  OF  SINUSOI- 

C  DAL  UAUCS  WITH  AT  HOST  TWO  FRCGUCNCV  CONPOICNTS. 
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C  OlCOOCt  •  T>C  CVSO  OCOQCI  sutiKvri>c 

e  eCCOKl  •  THE  CUSO  KCOCCX  SUMOUTINC 

riUTER  •  THE  SL'SPOUTINC  THAT  FILTfPS  THE  INPUT  TT*E  ECTIES  SAflnES 
USINO  THE  FILTEjI  CCEPFICIENTS  CCNESATEO  tV  FXTRSEH. 

POUEX  •  A  ROUTINE  TO  CALCtItATE  TX  POIIER  IN  A  SAflPUD  TIW  SERIES 
UlTH  IRPECCNCE  •  C««  CHTS. 

CtmtttmtttttrxtxmtttntirtxtttxmtttrmttxtsxrttttxttTtttttstttx 

Ctnxxxx3.xxtxxixxxxxxtxxtxiixxxxxxxxxxxxxxxxxTxxxxxxx:xiixxxxxxxxxixxxxxx 

C - PROCfi«(  START - - - 

C -  INITIALIZE  VARIAILES  AND  ARRAYS 

DtnptSION  AnPKZSS).  SNR(2M),  TSINtS***),  TSOUT<S«M),  ERRCBMS) 
t  *1  ) 

INTECER  lINOUTCStMt 

A(DS«*)  .  6i*T(lS.  XSHBXm  S  .Ml  t  SM.)  t  SOrTCE.) 

ICHAR  «  -1 

C -  IXUT  AND  PRINT  UORKINO  UHRIAtUS 

READ  S.  PREOl.  PS 
READ  R.PCl,  TC 
READ  t.lETA,  CV-A 

PRINT  *,  •  t;3T  at  •,PRE31.»  dim  and  '.PS,*  IPS* 

PRINT  »,*  UITH  TC  •  •.•'C 

PRINT  ».•  FILTER  PARANETERS  ARC.  SETA  •  •.SETA,  •.MAM  •  •.MfW 
e-~  OEJCRATf  OUTPUT  PILTER  CC4TPICIENTS 
OAU  FlTR«EN(SrTA.QAmA.NP.Sl 

C - -  INITIALZC  PLOTTER 

■ALL  PACTOR( .51 
caul  PLOTiE..  2..  -3$ 

C - START  LOOP 

00  i«M  tm  •  E.s.a 

RATIO  •  3».  ♦  fp 

Call  v«axopt(una)(,uaih.ps,pci,tc,ratio> 

IP  «  ICH’R  .CC.  •)  00  TO  6 
EYRX  •  ONAX 

com  •  WHIN 
•  CONTINUE 

c -  start  of  SIONAl-TO-NOISf  LOOP 

DO  3M  K  •  t.lM 
ANPl(K)  •  -49.  ♦.ASK 
PCAKl  •  A(A.'V>t(K>) 


C -  GOCHATC  INPUT  TINE  SCRIES  SAWPLES 

CPU  SIGNAL(T8IH,5«M.F8.PRC0l.R..PeM(l,R. ) 

C -  PROCESS  THE  INPUT  TIHC  SERIES  THROUGH  TX  CVSO  8VSTER 

CALL  CNCOD*:t<TSIN,iIKXrT,5BOS.PS.PCl.FC3,PC3.TC.EWX.EUnN.DC> 
CALL  OCCOOEl(IINOUT,TSCiUT,5e««.FS,FCl.FC3.FC3,TC,U3AX.WHIN,OC) 

C -  FILTER  TX  OUTPUT  OP  TX  OECODER 

CALL  FILTERCTSOUT.SAM.X.I) 

C -  DELAY  TX  INPUT  SIGNAL  START  TO  CORRESPOND  TO  TX 

C  FILTERRED  OUTPUT. 

DO  3«  ID  •  l.AtDS 
KD  •  EM  ♦  10 
TSIN(IO)  •  TSIN(KO) 

M  CONTINUE 

C—  ADJUST  OUTPUT  SIQNAl  AMPLITUDE  SO.  OUTPUT  POUCR  ■  imUT  POUCR 

CALL  POUCRCTSIN.MM.PS.PIN) 

OAU  POUCR(TSOUT,a:36,FS.POUT) 

RAIN  •  SC«T  <PIN/POUT) 

DO  40  I  •  t.4:.s 
rsouTdi  •  Tsourci)  s  rain 
4#  CONTINUE 
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C -  CMXUtATC  7«  NOIK  KMOI 

DO  M  t  •  1.4?;8 

•  TSOOTd)  -  rSINd) 
M  COffTIffl.>E 

CMX  P0ue)l(ER)t.4M«,Ft.CltV> 

0 -  CMXULATC  TW 


$MiKl  •  1«.  I  ALOQt*  (PTN  /  rU9) 

IF  dCHftR  .LT.  •)  CO  TO  fO 

SN«d»l)>  S.NRdC)  •  MldDl) 

••  ^0HTIf<u€ 

3M  COMTIhUC 

C - PRIMT  AND  PtXT  TMC  RESULT* 

lOMR  •  ICHAR  *  t 
IF  dCHAR  .GT.  a )  CO  TO  CM 
CALL  SCALE<^.C?l.l9..SC«,n 
CALL  SCAL£<S.'-?  £•  10,1) 

^MLL  AXIS«e.,#,,ii;>v4';pLITU0E  (DSM),-l(,l*..*..ArRl(l*l ),MPld*2) 
CALL 

N*  CO^INLE 

CAU  L1V£(A,TP1.»«.1**.1,1*.ICHIW) 

IM*  CO.'fTD^.'E 

CALL  PLOTC(n) 
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APPKMDIX  W 


V’-i  '■'(■'0  'Voqii';  f  7  I  r-o'  i' m 


PROCRrtrt  0IFGAIN(1HPUT,0UTPUT,TAPEG*CHJTPUT,P10T) 

^ . . — OIFFERNTIAL  OAIH  PROGRAM - 

V  THH  program  measures  the  SVSTEM  CAIN  US,  FREOUENCV  RESPONSE  FOR 
•;  M  C'.'SO  DICITAL/ANALOG  svstem  connected  BACIC-TO-SACIC. 

C  H  NhXIMALLV  flat  linear  phase  fir  filter  is  placed  on  the  OUTPUT 
C  OF  THE  DECODER  TO  REMOUE  SIGNAL  CC?.PON£NTS  A30UE  3308  HZ, 

UARIABLES  ttttxtxsxsttxxtxxxxxtittttxxxsxxsxtt 

C  FREOl  •  AN  ARRAV  CCNTAININQ  THE  FREOUENCIES  THAT  THE  SHR  HAS  BEEN 
C  MEASURED  AT.  THE  RANGE  IS  380  HZ  TO  3308  HZ. 

C  TSIH  •  AN  ARRAV  CONTAINlNQ  THE  INPUT  TIME  FUNCTION  SAMPLES. 

C  TSOUT  •  AN  ARRAV  CONTAINING  FIRST  THE  DECODER  OLTPUT  TIME  FL'NCTICN 
C  SAMPLES,  THEN  THE  OUTPUT  TIMS  FUNCTION  SAMPLES  OF  THE  FIR  FIL- 

C  TER. 

C  •  •  AN  ARRAV  CONTAINING  THE  FILTER  COEFFICIENTS. 

C  BINOUT  •  AN  ARRAV  CONTAINING  THE  BINARY  OUTPUT  OP  THE  CUSD  ENCODER 

C  AMPl  •  THE  AMPLITUDE  OF  THE  TEST  SIGNAL  IN  DBMS. 


C  FS  •  THE  SAMPLE  RATE. 


C  FCl,  FCa,  FC3  •  THE  ROLL-OFF  FREQUENCIES  OF  THE  PRIMARY  INTEGRA- 
C  TORS. 

C  TO  •  THE  TIME  CONSTANT  OF  THE  SYLLABIC  FILTERS. 

C  UMAX  L  UMIN  •  THE  MAXIMUM  AND  MINIMUM  INPUTS  TO  THE  SYLLABIC  FIL- 

C  TER. 

C  BETA  •  THE  NORMALIZED  CENTER  FREQUENCY  OF  THE  OUTPUT  FILTER  TRANS- 
6  ITION  BAND. 

C  OAMMA  •  THE  NORMALIZED  UIDTH  OF  THE  ROLL-OFF  REGION  OF  THE  OUTP'JT 
C  FILTER.  TrE  KEOICN  IS  THE  FREQUENCY  BAND  BETWEEN  THE  95*  AND 

C  5«  OUTPUT  AMPLITUDES. 

C  PEAKl  •  THE  MAXIMUM  AWLITUDE  OF  THE  TEST  SIGNAL  IN  UOLTS. 

C  MP  •  THE  NUMBER  OF  FILTER  COEFFICIENTS. 

C  KN  •  THE  NUMBER  OF  TEST  FREQUENCIES. 


C  DC  •  THE  DUTY  CYCLE  OF  THE  SLOPE  OVERLOAD  DETECTOR. 


Ctsxxxtxxxxxxxxtxmxxxxzxxxxxxxzxxxxxtxxxxxxxxxxxxxxxtixxxxxxxxxxxxxtxx 

eXXXXXSXXXtXXXXXXXXXX  SL3R0UTINES  USED  XiXXSXXXXXXUXXXXXXXXXXXXXXXXiXXS 


c 

■% 

k 

K 


FLTRGEN  •  THE  SUBROUTINE  THAT  GENERATES  THE  OUTPUT  FILTER  COEFFI¬ 
CIENTS. 


»LiT,  SCALE,  AXIS,  RECT,  LINE.  PLOTS,  •  CALCCMP  PLOTTING  ROUTINES. 

'IGNAL  •  the  TEST  SIGNAL  GENERATOR.  PRODUCES  SaMPLES  OF  SINUSOI¬ 
DAL  WAVES  WITH  AT  MOST  TWO  FREQUENCY  COMPONENTS. 


ENCODE!  '  the  CVSD  ENCODER  SUBROUTINE 


V 

C 


DECOOEl  •  the  CVSD  DECODER  SUBROUTINE 

FILTER  •  the  SUBROUTINE  THAT  FILTERS  THE  INPUT  TIME  FUNCTION  SAM¬ 
PLES  USING  THE  FILTER  COEFFICIENTS  GENERATED  BY  FLTRGEN. 


C  POWER  ■  A  ROUTINE  TO  CALCULATE  THE  POWER  IN  A  SAMPLED  TIME  FUNC- 
C  TION  WITH  IMPEDENCE  ■  680  OHMS. 


estxtstixsxxsstsxxxsxsssxxtxxxsxxxsiststxxxxxxsstxtxxxatsnxxxxtxxxtxxst 

CXXXXXXXttXitXXXtXXiXXXXXSXXXXSXXXXXSXXXXXXXStXXttSXXXXXXXXXXXXXXXXXXXXX 
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c - PfiOGRMI  START - 

C -  INITIALIZE  yAfilABLES  AND  ARRAYS 

DinCHSION  rREQKiea),  TSIN<5««e),  TSOOT<SMe) 

I,S(CM),  CAIN(c‘;0) 

INTEGER  BINOUT(5C30) 

ACDBTW)  •  SGRT(18.  «(<OB(W  t  .Ml  S  6M. )  S  SORT(8.] 

ICHAR  •  -1 

C -  INPUT  AND  PRINT  THE  WORKING  VARIABLES 

READ  t,  ANPl.  FS 
READ  «,FC1.  TC 
READ  B.BETA,  G»rf»A 

PRINT  »,  ■  S^•R  TEST  AT  '.ANPl,*  DBM  AND  ‘.FS,’  BPS* 

PRINT  »,  •  WITH  TC  •  '.TC 

PRINT  »,*  OUTPUT  FILTER  PARANETERS  ARE,  BETA  •  •,BETA.*  GANNA  * 
INNA 

C -  GENERATE  OUTPUT  FILTER  COEFFICIENTS 

CALL  FLTRG£N(BETA.CANnA,NP,B) 

C -  INITIAL2E  PLOTTER 

CALL  FACTOR (.5) 

CALL  PL0T(2.,  2.,  -3) 

C -  START  LOOP 


DO  aOM  NTINES  •  2.6,2 
ICHAR  ■  ICHAR  *  1 
RATIO  •  30.  ♦  NTirSS 
hLL  ' RIAXOPT ( UfWX .uniN.FS.FCl.TC, RATIO) 

i.H  •  J 

.* —  start  of  signal-to-noise  loop 

00  300  K  •  300,3600,100 
KN  <  KN  4  1 

C -  GENERATE  TEST  SIGNAL 

FREOKKN)  •  K 
PEAKl  •  A(A,"?H 

CALL  S1GNAL(TSIN,5000,FS,FRE01(KN),0.,P£AK1,0. ) 

C -  PROCESS  THE  INPUT  TINE  FUNCTION  THROUGH  THE  CVSD  SYSTEN 

CALL  EHCODE1<TSIH,BIHOUT,S3-30.FS.FC1,FCH,FC3,TC.VNAX.VNIH  DC) 
CALL  DEC00E1(BIN0UT,TSCUT,5G30.FS,FC1,FC2,FC3,TC,U;".AX,W1IH,DC) 

C -  FILTER  THE  OUTPUT  OF  THE  DECODER 

CALL  FILTER<TSOUT.5C0O,rtf»,B) 

C -  DELAY  THE  INPUT  SIGNAL  START  TO  CORRESPOND  TO  THE 

C  FILTERREO  OUTPUT. 

DO  30  ID  -  1.4096 
ICO  ■  200  4  ID 
TSIN(ID)  «  TSIN(KD) 

30  CONTINUE 

C -  FIND  SYSTEN  GAIN 

CALL  P0UER(TSIN.40S6.FS.PIH) 

CALL  POUER(TSOUT,40SC,FS.POUT) 

CAIN(KN)  >  10.  S  ALCGiecPOUT/PlN) 

3M  CONTINUE 


.CA 
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^  C -  ADJUST  MIN  yftLUES  70  8M  H2  REFtRERENCE 

DO  6  X  •  I.ICN 

IF  (FREOi(i)  .E9.  8M.)  REFGAIN  ■  MIN(I) 

6  CONTIUJE 

DO  8  I  •  1,KN 

CAINCI)  •  MIN(I)  -  REFMIN 
8  CONTINUE 

C -  PLOT  THE  RESULTS 

IF  (XCNAR  .07.  8)  GO  TO  S8« 

CALL  SCflLE(CAIN.6.,ICN.n 
CALL  LG':C«L(fAril.ia..'C.S) 

CALL  LCA)(ISia.,0.,i<r^k.;GCENCV  (KZ).-i4, l«. ,8. .FREOldOHD.FRECl (K 
lNt2)) 

CALL  AXIS(e.,e..22KDIFFERENTIAL  CAIN  (D> ),22,6.,M.,GAIN(ICN41 
1GAIN(»CN'^2)> 

CALL  R£CT(e.,9..6..18..e«>3) 

8«e  CONTINUE 

CALL  LGLlN£(n?£01,CAIN,ILN,l«,ICNAR,-l) 

ZDM  CONTINUE 

CALL  PLOTE(N) 

END 
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APrt'.MDiy  X 


**.'  t. ’  ’lit  y  pr*  :  I '  *K‘ !  rc)  — 

PROCRAn  WIDCAINdHPUT, OUTPUT. TAPE6-0UTt>UT,PL0T) 

— . . niSHATCHEO  SYSTEW  CAIN  RESPONSE - 

THIS  PROCRAN  PEASt'RES  THE  SVSTEW  GAIN  US.  EPEQUENCV  RESPONSE  FOR 

A  CVSD  DIGITAUAfiALOC  SYSTEN  CONNECTED  BACIC-TO-BACIC. 

A  NAXirWLLV  FLAT  LlNEfR  PHASE  FIR  FILTER  IS  PLACED  ON  THE  OUTPUT 

OF  THE  DECODER  TO  REHCUE  SIC.NAL  CC«>CN£NTS  ABOUE  3S0«  HZ. 

UARiABLES  tttttxxximtttttttttxttxxtxtxtttxtx 

FREOl  ■  AN  ARRAY  CCNTAININO  THE  FREQUENCIES  THAT  THE  SNR  HAS  BEEN 
NEASURED  AT.  THE  RANGE  IS  3&a  HZ  TO  2338  HZ. 

TSIH  •  AH  ARRAY  CONTAINING  THE  INPUT  TIME  FUNCTION  SAMPLES. 

TSOUT  •  AN  ARRAY  CONTAINING  FIRST  THE  DECODER  OUTPUT  TIME  FUNCTION 
SAMPLES.  Th€N  THE  OUTPUT  TIME  FUNCTION  SAMPLES  OF  THE  FIR  FIL¬ 
TER. 

•  •  AN  ARRAY  CONTAINING  THE  FILTER  COEFFICIENTS. 

•INOUT  •  AN  ARRAY  CONTAINING  THE  BINARY  OUTPUT  OF  THE  CUSD  ENCODER 

AHPl  •  THE  AMPLITUDE  OF  THE  TEST  SIGNAL  IN  DBMO. 

FS  •  THE  SAMPLE  RATE. 

FCl.  FC2.  FC3  •  THE  ROLL-OFF  FREQUENCIES  OF  THE  PRIMARY  INTEGRA¬ 
TORS, 

TC  •  THE  TIME  CONSTANT  OF  TW  SYLLABIC  FILTERS. 

UMAX  1  UMIN  •  THE  MAXIMUM  AND  MINIMUM  INPUTS  TO  THE  SYLLABIC  FIL¬ 
TER. 

BETA  -  THE  NORMALIZED  CENTER  FREQUENCY  OF  THE  OUTPUT  FILTER  TRANS¬ 
ITION  BAND. 

CAFV1A  •  THE  NORMALIZED  WIDTH  OF  THE  ROLL-OFF  REGION  OF  THE  OUTPUT 
FILTER.  Th-  KiGICN  IS  THE  FRcQU£.NCY  BAND  BETWEEN  THE  95«  AND 
SR  OUTPUT  AFPLITUDcS. 

PEAKl  •  THE  MAXIMUM  AMPLITUDE  OF  THE  TEST  SIGNAL  IN  UOLTS. 

»•»  .  THE  NUMBER  OF  FILTER  COEFFICIENTS. 

f.;i  •  THE  HUMBER  OF  TEST  FREQUENCIES. 

DE  •  THE  DUTY  CYCLE  OF  THE  SLOPE  OVERLOAD  DETECTOR. 

t«itttttttt«tt»x(tiztxtxtt»xxttxtsxxszssisxnHsxtx»txttx<iitti>stxs 
txxtxtf IIXXXIXXXXZXX  SUBHCUTIHES  USED  XXX32ZZ2XXZXXSXZXXXXXZZXXXXXZXSXS 

FLTRCEN  •  THE  SUBROUTINE  THAT  GENERATES  THE  OUTPUT  FILTER  COEFFI¬ 
CIENTS. 

PLOT,  SCALE,  AXIS,  RECT,  LINE,  PLOTE,  •  CALCOMP  PLOTTING  ROUTINES. 

SIGNAL  •  the  test  SIGNAL  GENERATOR.  PRODUCES  S«MPL£S  OF  SINUSOI¬ 
DAL  UAL€S  WITH  AT  MOST  TWO  FREt  ENCY  CCrPONENTS. 

ENCOOCl  •  THE  CVSD  ENCODER  SUBROUTINE 

DCCODCl  •  THE  CVSD  DECODER  SUBROUTINE 

FILTER  •  THE  SUBROUTINE  THAT  FILTERS  THE  INPUT  TIME  FL'NCTION  SAM- 
PUS  USING  THE  FILTER  COEFFICIENTS  GENERATED  BY  FLTRGEN. 

POUER  •  A  ROUTINE  TO  CALCULATE  THE  POWER  IN  A  SA«»LED  TIME  FUNC¬ 
TION  WITH  IMPEDENCE  •  6M  OHMS. 

txxxxxxxxxxxxxxxxxxxxtxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxttxxxxxttxxtxxxxxx 

ettXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXii:iXXXXXtXtXXXXXXXXXXiXZXXX 


- - PROCRftH  STPRT - - - 

C -  IHITIPLIZE  ypRIABLES  AND  ARRAYS 

fREOKica),  TsiNcsees),  Tsourcsew) 

C;^Ih(r:3) 

IKTECER  BIhOUTCb,)Oa) 

1^'^’  ;jSCST(10.  xt((DBR«  -4.  )/!•.)  S  .Ml  t  6(M.)  t  S«IT(2. ) 

5 -  INPUT  AMD  PRINT  THE  UORKINO  UAR1A8LES 

READ  t,  ANPl,  FS 
READ  t.FCl.  TC 
READ  X.BETA,  CF.-CiA 

wKl  i: :  ^  -.rs.-  »ps- 

j^MT  t,*  OUTPUT  FILTER  PARAMETERS  ARE,  BETA  •  •,BETA,*  CAMMA  ‘.CA 

; -  CEMERATE  OUTPUT  FILTER  COEFFICIENTS 

CALL  FLTR<iEN< BETA, GAMMA. NP,B) 

-  INITIALZE  PLOTTER 

rtLL  FACT0«(  .S) 

•hLL  PLOTia.,  a.,  -3) 

-  JTmRT  LOOP 


DO  aooe  NTIMES  •  2,6,2 

ICHAR  •  IChAR  ♦  1 

RATIO  •  30.  ♦  NTirSS 

CALL  unA>COPT(VnAX,umN,FS,FCl,TC,RATIO) 

KN  •  e 

IF  (ICHAR  .GT.  9)  CO  TO  IM 
EWW  •  UTAX 
EVMN  •  UMIN 
IM  CONTINUE 

C -  START  OF  SIONAL-TO-NOISE  LOOP 

DO  3M  K  •  3M,3&M,1M 

KM  ■  KM  T  1 

C -  GENERATE  TEST  SIGNAL 

FREQKKN)  •  K 
PEAKl  •  ACAfIPl) 

CALL  S1GNAL(TSIN, 5000, FS,FREQl(KM),e., PEAKl, e. ) 

C -  PROCESS  THE  INPUT  TIME  FUNCTION  THROUGH  THE  CVSD  SYSTEM 

CALL  ENCOOEl<TSIN,BlNOUT.5CO0.FS,FCl,FC2,FC3,TC,EVNX,ELrN,DC) 
CALL  DECODEl<BINOUT,TSOUT,500a,FS,FCl,FC2,FC3,TC,UiVlX,uniN,DC) 

C -  FILTER  THE  OUTPUT  OF  THE  DECODER 

CALL  FILTER(TSOUT,S0M,NP,i) 

C -  delay  the  INPUT  SIGNAL  START  TO  CORRESPOND  TO  THE 

C  FILTERREO  OUTPUT. 

DO  30  ID  •  1.4096 

KO  >  200  ♦  ID 

TSINdDI  •  TSIN(KD) 

30  CONTINUE 

C -  FIND  SYSTEM  GAIN 

CALL  PO«ER(TSIN. 4096. FS, PIN) 

CALL  POUERCTSOUT.^eJS.FS.POUT) 

QAIN(KN)  ■  10.  «  ALCC10(POUT/PIN) 

3M  CONTINUE 

C -  ADJUST  GAIN  UALUES  TO  BOO  KC  REFERERENCE 
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IF  (FREQl(i^^*EO*  80®. )  REP<1A1N  •  6A1N(X) 
6  CONtIHUE 

CAlN(i)  -^CAlHd)  -  REEfiAIH 
8  COMTINUE 


c -  plot  the  results 


0M 

bVW 


ArF'’:fti  TV 


Tent  Sj ''nal  r-'  ;  I';” 


5Ut«0UTirC  SICMAL(0irmjr,N,rt.FREQt,FR£0I.M91,MVS) 

c  1 1 1 « » • « t  *  1 1 1  ( 1 1 1 1 1 1 «  f  1 1 1 « tsnttusnsnsniuxsssnnnxxsxnttstunn 

C  this  >0»»0UTINE  CEMCRATES  R  TEST  SICNRL  COf'i'OSCB  OF  UP 
C  TO  TUO  SIK  UAues  OF  OIFFEREMT  FREGUt.'tCIES  WiD  PftPLXTUDCS. 

Cttttttutttttittttttxxtt  uAAiRltXS  tzixxxttxxtxnttxutsxtttxxtttmtxM 

C  OUTPUT  .  dti  PRRftV  COMTAIMINQ  THE  OUTPUT  TIPS  FUHCTIOM  SPFPLES 

C  H  •  TIC  NURBER  Of  SAPPUS  Of  THE  TIME  FUNCTIOM  DESIRO 

C  fS  •  nC  SAPPU  RATE  IN  K»yS 

C  fREOl  •  THE  FSEOUENCY  Of  TIC  FIRST  SIONPL  COPPOWENT 

C  FRE02  •  THE  FREOUENCV  CF  TIC  SECOfO  SlCHAl  COWONEKT 

C  MVi  •  THE  PE«  ANPtlTUOE  OF  THE  FIRST  SIONPL  CCfPOICMT 

C  MVZ  •  THE  PEAK  APPLITUtE  OF  THE  SEC0H9  SIGNAL  COPPi^CNT 

CMtxtxtstttxtttixzixtsxzxiXi3xtizxizxiaxxixixtzizx*x3Xi>  xxsxxxxxtsxxxxx 

C -  INITIALIZE  UARIABUS  ANB  ARRAVS 

DINEHSIOH  OUTPUT(N) 

DATA  Pl/Z.lStSi^SSaS/ 

C -  GENERATE  OUTPUT  SANPU8  Of  TEST  SIGNAL 

DO  56  I  •  t.N 

OUTPUT! I)  •  APPl  t  SIH<  a.  »  PI  1  FRE91  /  FS  «  I ) 

!♦  FNPa  t  SIN  <2.  f  PI  «  FHEGa  /  FS  »  I) 

6S  CONTINUE 

RETURN 
END 
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VO 


ArrKr'ii.i;:  / 


■'n-il 


1  ,'  i1  ;  I  ioi;  v.^roi  •  ! 


JUWOUTIK  POWtHfX.M.FS.^) 

. — <WCRA«  POUER  9UWKXni«' 


This  iusfwriNE  wt.cuu»n:9  tve  pot^n  or  the  i^cvt  ti?* 

fufiCTiON  s  F-^.-oss  A  i.  i  c  -i^  li -• -r'.-:c£.  TVS  itprrr  Tir>£  fiuc- 

TI0t>  SAftHUS  «*)£  1I-90T  TO  Th£  i'JlS<Xm^£  Ai  A  1  X  H  A-iJAV. 


itttttttttttttxttitxtttx  OAOlAtLE)  s»tt3n»xxtsxstxx«tttzn»stiMn<s 
C  X  •  A«  ARRAV  COKTAIMINO  THE  IKPUT  Tir<  rLHCTION  SAMPLES 
i  H  ■  THE  NUfOER  Of  SAWTiS  TO  It  PROCESSED 

c  FS  •  the  sakpu  rate 

C  P  •  THE  CALCUWTED  SICMAt  POWER 

cxtttttttxixtxxntttxxxxtxxmsxxtxxxzttxttttxatxxxtxttxtitxstxxtttsxn 

c - SUtaOUTIKE  START - 

C -  IKITIALIZE  UARIARUS  AMD  ARRAYS 

DIREHSIOH  X(N) 

sun  • 


C -  sun  TVC  SQUARES  OF  nc  INPUT  TINE  FUNCTIOW  6AI«>US 

DO  IS  I  •  l.N 

sun  •  $<.*1  ♦  xcn  »  2. 

IS  COffTINUC 

C -  CALCULATE  THE  AUCRAOE  POUER  ACROSS  A  SSS  OHH  inPEDCNCE 

p  •  sun  /  N  /  6SS. 

RETURN 

END 
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